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1 . 1  MBE  j rowth  and  Characterization* 

During  this  period,  a  considerable  effort  was  placed  on  the  preparation 
and  characterization  of  single  and  multiple  period  modulation  doped 
( A1 ,Ga)As/GaAs  heterostructures.  A  theoretical  study  of  these 

he taro structures  was  also  undertaken  and  a  good  agreement  with  experiments  was 
obtained.  Segregation  effects  of  Sn  in  GaAs  grown  by  MBE  was  also 
investigated  in  detail. 

Single  and  multiple  period  modulation  doped  structures  with  mobilities 
several  times  those  reported  by  other  laboratories  were  prepared.  Transport 

properties  of  these  structures  were  studied  at  low  and  medium  electric  fields 

and  at  room  and  cryogenic  temperatures.  The  maximum  room  temperature 

mobility,  8,400  cm* /Vs,  wa3  obtained  when  the  AlAs  mole  fraction  was  about 
3SH.  In  addition,  the  single  period  structures  so  far  appeared  to  exhibit 

*  This  work  was  supported  bv  the  Joint  Services  Electronics  Program 

(U.S.  Army,  U.S.  Vavy,  U.S.  Air  Force)  under  contract  N0001  4-79-0-0424 ,  and  by 
the  Air  Force  Office  of  Scientific  Research  under  contract  AFOSR  80-0084. 
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larger  room  temperature  mobilities  then  the  multiple  period  ones.  At  nn  '<  end 
10  '<,  however,  comparable  results  were  obtained.  77  '<  mobilities  of  100,000 
cm -  h]3  end  10  M  mobilities  of  about  000,000  -  210,000  cm* 'Vs  were  obtained 
both  with  single  and  multiple  period  structures.  If  the  AIA3  mole  fraction  is 
made  larger  than  about  ^3"o,  the  direc  t- indirec  t  transition,  the  electron 
mobilities  steadily  decline  up  to  the  AlAs  mole  fraction  of  lOOt. 

Secondary  ion  mass  spectroscopy  (SIMS)  has  been  used  to  investigate  the 
accumulation  of  3n  at  the  growing  film  surface  and  the  corresponding  3n 
depletion  near  the  film- substrate  interface  during  MBE  deposition  of  Sn-doped 
3aA3.  3n  surface  concentrations  3gn  more  than  two  orders  of  magnitude  larger 
than  the  3teady  state  3n  doping  level  (Cgn)  have  been  observed.  9gn  n-n  -ills 
grown  at  approximately  constant  As^Sa  surface  concentrations  was  found  to 
increase,  while  remained  constant  with  increasing  growth  temperatures  for 
?3  C,  530°C.  Both  and  0gn  decreased  at  higher  T3  values.  These  results, 

as  well  as  those  reported  by  previous  investigators,  have  been  explained  by  3n 
surface  segregation  during  deposition  in  which  the  segregation  is  driven  by 
the  relief  of  strain  energy  due  to  the  incorporation  of  over-sized  3n  atoms. 
Calculated  Sn  profiles  were  3hown  to  provide  a  good  fit  to  experimental  data. 

• . 2  Transport  ? roperties  of  Hetero junction  3 aye rs* 

Transport  properties  of  modulation  doped  structures  up  to  an  electric 
field  strength  of  2  bV-'cm  were  investigated.  The  differential  electron 
mobility  drops  suite  rapidly  up  to  200  7 ''em  and  then  stays  about  constant  up 
to  2  <7 'cm.  Higher  electric  fields  cause  negative  differential  resistance 

effects  -which  can  result  in  nonuniform  electric  fields.  Once  the  field  is 
nonuniform,  no  reliable  u(E)  data  can  be  obtained.  The  electric  fields  in  1  u 
gate  no rm ally-off  ?ETs  are  estimated  to  be  about  2  kV/cm;  therefore,  our 
results  should  be  applicable  to  predicting  device  performance.  Structures 
with  an  extremely  high  low-field  mobility  showed  substantial  decrease  with 
electric  field,  particularly  at  low  temperatures.  Generally  the  room 
temperature  mobility.  3,A00  cm- -'Vs,  was  maintained  up  to  2  kV'cm.  However,  at 
'?~r  mobility  at  2  <V/'cm  was  about  23-!0co  of  its  zero- field  value.  The  '0  '< 

*  This  work  wa3  sunoorted  by  the  Air  Force  Office  of  Scientific  Research  under 
contract  A?0SR  33-0031,  the  Office  of  Vaval  Research  under  contract 
70001 i-7o-0-0306,  and  by  the  Army  Research  Office  under  contract  DAAO 
29-30-0-0011. 


3 


1 .  MOLECULAR  BEAM  EPITAXY 

mobility  at  2  kV/cm  was  about  15-20«  of  its  zero- field  value.  These  results 
indicate  that  a  differential  mobility  of  about  25,000  cm^/yg  at  2  kV/cm  with 
an  electron  concentration  of  about  IQ^cm"-  can  be  obtained  at  77  Thi3 

represents  an  intrinsic  speed  improvement  of  5  and  a  speed-power  product 
improvement  of  25  if  normally-off  ?ETs  were  fabricated  and  operated  at  ~’’7 
Room  temperature  operation  can  result  in  an  improvement  of  50^  and  about  200* 
in  speed  and  speed-power  product,  respectively.  3y  "improvement"  we  mean  to 
compare  it  to  a  conventional  normally-off  GaAs  PET  using  a  bulk  channel  layer. 

We  have  done  considerable  work  during  this  period  on  the  real- space 
electron  transfer  effect,  in  which  electrons  in  GaAs-AlGaAs  heterostructures 
transfer  at  high  fields  from  the  GaAs  to  the  AlGaAs  layers.  The  resulting 
negative  resistance  as  been  simulated  by  Monte  Carlo  techniques  and  has  been 
observed  experimentally.  This  effect,  which  has  great  potential  for  device 
applications,  is  discussed  further  in  Unit  4. 

1  . 3  Theory  of  Surfaces  and  Heterojunction  Interfaces* 

The  deep  levels  associated  with  surface  defects  in  compound 

semiconductors  have  been  studied  theoretically;  the  results  indicate  that 
Fermi-level  pinning  and  Schottky  barrier  formation  is  produced  by  surface 
antisite  defects  in  most  cases,  rather  than  by  surface  vacancies  or  subsurface 
defects.  We  also  have  investigated  the  possibility  of  designing  Schottky 

barrier  heights,  and  barrier  heights  at  hetero junctions  and  in  superlattices, 
via  the  deposition  of  selected  impurities.  Predictions  of  intrinsic  surface 
states  and  surface  core  excitons  in  III-V  and  II -VI  semiconductors  were  found 
to  be  in  excellent  agreement  with  the  available  experimental  data.  Finally, 

it  was  found  that  impurities  that  are  shallow  donors  or  acceptors  in  the  bulk 

can  become  deep  traps  at  semiconductor  interfaces. 


*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
('J.S.  Army,  U.S.  Navy,  U.S.  Air  Force)  under  contract  N0001 4-79-3-Oi24. 
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3.  3.  Stillman 


3raduate  Students 
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2 . 1  Intro ducting 

The  primary  goal  of  this  research  i3  to  evaluate  in  a  potential  of 
high-purity  In?  and  InCaAsP  for  optical  and  microwave  device  applications. 
3esiies  investigating  the  various  growth  parameters  and  techniques  that  can 
influence  the  purity  of  the  grown  layers,  an  important  part  of  this  work  i3 
the  characterization  of  the  purity  and  crystalline  quality  of  the  samples 
grown.  The  total  impurity  concentration  is  studied  through  Hail  coefficient 
measurements  at  3GO°!<  and  at  77°X.  The  origin  3.od  influence  of  growth 
conditions  on  acceptor  concentrations  are  studied  through  low  temperature 
photoluminescence  measurements.  Deep  levels,  and  the  influence  of  different 
growth  parameters  on  the  concentrations  and  energy  levels  of  these  centers, 
will  be  studied  through  DLTS  measurements.  The  crystalline  quality  and  degree 
of  lattice  match  are  studied  through  x-ray  diffraction  measurements.  The 
residual  shallow  donor  impurities  are  studied  using  far  infrared  Fourier 
transform  measurements  of  photothermal  ionization  spectroscopy.  This 
technique  is  being  used  to  study  the  residual  donor  impurities  in  high-purity 
3aAs  and  InP.  Far  infrared  photoconductivity  measurements  on  In 3a As P  and 
In3aA3  have  been  used  to  determine  the  effective  mass  in  these  materials,  and 
this  technique  may  be  useful  for  studying  alloy  disorder  in  these  materials. 


"*•*'  ^  -  -  - 
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2 . 2  Hydride  Vapor  Phase  Growth  of  High  Purity  InP* 

The  vapor  phase  system,  constructed  with  funds  from  3NR,  is  now  fully 
automated  and  under  computer  control  for  added  reproducibility.  Previously 
high  purity  GaAs  (^77  ^  *  38,000  cm^/V-sec,  1177  ^  =  1.7  x  lO1^  c®""^ )  had  been 
grown  in  this  system  but  comparable  InP  had  not.  Difficulties  encountered  in 
growing  InGaAs  lattice  matched  to  InP  led  to  a  more  detailed  and  systematic 
analysis  of  the  system's  growth  parameters  for  InP. 

Prior  to  growth,  the  substrates  are  etched  in  situ  with  a  D.131?  HC1  in  H? 
mixture  to  provide  a  clean  surface  upon  which  to  initiate  growth.  In  addition 
the  Arsine/Phosphine  line  has  been  extended  further  down  into  the  gas  stream. 
The  growth  rate  has  increased  to  the  point  where  it  is  possible  to 
substantially  decrease  the  total  gas  flow  (by  nearly  90^).  These  lower  flows 
result  in  much  less  downtime  for  cleaning  purposes  and  a  substantial  saving  in 
the  amount  of  source  material  used.  The  InP  grown  in  this  system  is  now 
comparable  to  that  grown  by  any  technique,  with  the  best  material  to  date 
having  a  liquid  nitrogen  mobility  U77  K  =  70, 300  cm-/V-sec  and  carrier 
concentration  n^y  K  =  5  x  I  01  4  cm-"5 . 

2 . 3  LPE  G rowth  of  InGaAsP** 

High  purity  InP  and  InGaAsP  epitaxial  layers  are  essential  for  wide 
depletion  width  and  low  capacitance  PIN  and  avalanche  photodiodes  and  for 
materials  characterization  in  general.  Also,  for  small  band  gap  InGaAsP 
alloys,  high  purity  material  is  required  to  avoid  Zener  breakdown  rather  than 
avalanche  breakdown  in  p-n  junctions.  Methods  for  obtaining  high  purity  InP 
and  InGaAsP  epitaxial  layers  on  (100)  InP  substrates  using  liquid  phase 
epitaxy  have  been  developed  and  this  high  purity  material  has  been 
characterized  through  Hall-effect  measurements  and  far  infrared  photothermal 
ionization  measurements.  Samples  with  net  carrier  concentrations  as  low  as  1 
x  I  01  7  cnT^  for  InP  and  2  x  1 01  ^  cm-^  for  In> ^Ga 4 47A3 ,  with  corresponding 
liquid  nitrogen  temperature  electron  mobilities  of  67,500  and  67,900  cm^/V-s, 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
(’J.S.  Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  Contracts  DAAG-2 9-73-C -001  6 
and  NQ001 4-79-C-0424  and  by  DARPA  Contract  N000014 -77 -C -00036 . 

**  This  work  was  supported  by  the  Naval  Research  Laboratory  under  Contracts 
N001  77-73-C-D1  29  and  N(X31 77-79-C-01  34 ,  and  the  Joint  Services  Electronics 
Program  (U.S.  Army,  U.S.  Navy  and  U.S.  Air  Force)  under  Contracts 
DAAG-2 9-73-0-001 6  and  N0001  1-79-C-0424. 
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2.  SEMICONDUCTOR  MATERIALS  AND  DEVILS 
respectively  have  been  grown  routinely. 

3oth  the  growth  solution  baking  technique  and  the  growth  temperature  have 
an  influence  on  the  epitaxial  layer  purity.  Baking  times  of  24  to  IS  hours 
are  required  to  reduce  the  carrier  concentration  to  these  levels. 
Incorporation  of  small  amounts  of  water  vapor  (<  .1  ppm)  in  the  hydrogen 
atmosphere  in  an  attempt  to  reduce  silicon  transport  from  the  quartz  tube  to 
the  growth  solution  during  baking  has  not  improved  the  purity,  and 
incorporation  of  larger  amounts  (>  I  ppm)  increases  the  carrier  concentration 
and  reduces  the  mobility.  An  increase  in  the  growth  temperature  has  improved 
the  purity,  especially  for  In0aA3,  presumably  due  to  a  decrease  in  the  silicon 
distribution  coefficient  as  the  growth  temperature  is  increased. 


2 . 4  Fourier  Transform  Spectroscopy  of  Shallow  Donor  Levels* 

Shallow  donor  energy  levels  have  been  measured  using  photo  thermal 
ionization  spectroscopy  in  high  purity  LPS,  Hydride  VPS,  and  bulk  InP,  grown 
in  our  laboratory  and  elsewhere.  This  .echnique,  which  involves  the 

measurement  of  p'notocond active  response  versus  incident  photon  energy  at  low 

temperature  (i  '<)  and  in  a  high  magnetic  field  (_<  5. 5  T)  has  also  been  applied 

to  high  purity  L?S  In Da As  and  InGaAsP  samples  lattice  matched  to  InP.  The 

minicomputer  based  Fourier  transform  spectrometer  facility  was  set  up  with  the 
support  of  the  Naval  Research  Laboratory. 

In  the  LPS  and  Hydride  '/PS  InP  grown  in  our  laboratory  and  the  LPS  In? 
grown  by  L.  F.  Eastman  at  Cornell,  the  sample  purity  wa3  sufficient  that  peaks 
corresponding  to  different  donor  species  could  be  resolved  in  the  phototherxal 
ionization  spectra.  Though  it  is  not  yet  possible  to  identify  these  donor 
species,  the  two  donors  present  in  the  Hydride  VPS  material  were  both  common 
to  the  Illinois  LPS  material,  and  one  of  these  was  also  the  dominant  donor  in 
the  Cornell  LPS  material.  A  third  donor  species  may  also  be  present  in  the 
Illinois  LPS  material  at  low  concentration.  It  is  reasonable  to  speculate 
that  a  donor  common  to  both  LPS  and  VPS  In?  may  be  silicon,  incorporated  from 
the  hot  quartz  growth  reactors,  especially  in  view  of  the  high  distribution 
coefficient  for  silicon  in  InP. 

*  This  work  wa3  supported  by  the  National  Science  Foundation  'under  Contract 
N3F  DHR-7"'-23329,  by  the  Office  of  Naval  Research  under  Contract 
NOOOt  4-77-C-C6cd  and  by  the  Joint  Services  Electronics  Program  ('J.3.  Army, 
•J.3.  Navy,  and  'J.3.  Air  Force)  'under  Contract  NOOQ1  4-79-C-042 4. 
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2.  SEMICONDUCTOR  MATERIALS  AND  DEVICES 

The  magnetic  field  dependence  of  the  peak  positions  in  these  spectra 

allows  accurate  determination  of  bottom  of  the  band  effective  mas3  (m*), 

* 

donor  ionization  energy  R  ,  and  low  temperature  static  dielectric  constant  e 
Values  of  these  important  materials  parameters  have  now  been  obtained  for  InP, 
In^Ga.^As  and  In.g^Ga  _  ^  -^Ga > ^qP  _  gQ.  For  these  materials,  the  respective 
values  determined  for  these  parameters  are: 

(m*)  =*  .00796  ^  .005,  .0437  ^  .001  ,  and  .0635  ^  .001 

R*  »  7.4  .17  meV,  4-30  +_  .17  meV,  and  5.29  .17  meV 

e0  »  12.06  +_  .13,  13.5  1  .5,  and  13.4  .5. 

The  peaks  in  the  spectra  of  the  ternary  and  quaternary  alloy  samples  are 
much  broader  than  would  be  seen  in  GaA3  or  InP  of  similar  liquid  nitrogen 
temperature  carrier  concentration  and  Hall  mobility  =  3  x  1 Q1  ^  cm~\  vw7 

3  53,000  cm^/V-s).  A  mechanism,  not  present  in  the  binary  alloys,  where 
compositional  alloy  disorder  results  in  a  different  local  alloy  environment 
for  different  donors,  could  account  for  the  excess  broadening  seen  in  the 
ternary  and  quaternary  spectra.  The  degree  of  this  broadening  might  then  be 
interpreted  as  a  measure  of  this  compositional  disorder. 
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3. ’  Excitation  and  3 tudy  of  Highly-Excited  V ibrational  States  of  Molecules 

' . 1 . '  Oirection  and  Significance 

Euantum  electronics  hai  it3  roots  in  the  microwave  3pe c tro.ecooy  of  low 
lying  energy  levels  of  molecular  gases.  Optically  p'umpai  far  I?,  lasers  hare 
employed  the  fundamental  vibrational-rotational  energy  levels  of  -,he 
molecules.  Mot  bands  and  excited  3tates  of  molecules  ao  to  the 
" quasi- continuum"  have  been  little  used  in  quantum  electronics  and  the 
spectroscopy  and  state  kinetics  are  relatively  unexplored. 

Ic  i 3  the  thrust  of  chis  work  to  study  the  kinetics  of  the  populations 
and  pi  lari cat ions  of  molecular  states  up  to  the  " quasi -continuum"  for 
potential  u3es  in  frequency  tunable  lasers,  nonlinear  optics  and  chemistry. 
To  implement  the  3tudy  requires  the  development  of  state-specific  excitation 
and  diagnostic  techniques,  the  initial  stage  of  the  wo~k. 

*  Ehis  wo  rk  wa3  suaported  by  the  Toint  Eervices  Electronics  Program 
.  H.E.  Army,  'J.E.  Mavy,  "M.E.  Mr  Eorced  under  contracts  EWE  ">  a-W-E -001  c  and 
H000 U-71-M1U,  and  by  contracts  MOOT  1-EO-EO? SO,  HAE\  HEE-1  Wa,  ME?  E?E 
EO-ObWE,  AEOEH  "a-oi  and  00E  Eandii  PE-10O0E1;. 
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3.1.2  Summary  o f  ?  rogress 

A  pulsed  multimode  DO.,  laser  has  been  used  to  excite  a  hot  band 
transition  in  room  temperature  DCS  gas  to  a  state  wi  th  three  quanta  of 
vibrational  excitation.  A  transient  microwave  gain  was  observed  and  is  a 
measure  of  the  excited  stats  population  difference  between  l  doublet  states. 
This  inherently  small  signal  has  been  observed,  without  signal  averaging,  by 
using  a  microwave  reflection  bridge  which  is  transiently  unbalanced  by  an 
infrared  pulse  exciting  the  gas  sample  in  part  of  the  bridge.  An  infrared 
pulse  truncation  shutter  and  lower  noise  microwave  detection  system  are  being 
added  for  direct  observation  of  the  excitation  relaxation. 

3 . 2  Energy  Transfer  ? recesses  in  Electronically  Excited  Atoms  and  Molecules 

3.2.1  TuOw-?ield  Ionization  Coefficients  of  Rare  Eases 

The  electron  ionization  and  excitation  coefficients  for  Ar,  Kr,  and  Te 
have  been  measured  in  the  low  E/p  region  (1-12  V / cm  torr')  using  a  drift-tube 
apparatus.  Discrepancies  between  previous  measurements  have  been  resolved  and 
sho  wn  to  be  due  to  photoelectric  emission  attributed  to  the  rare  gas 

metastables.  rJsing  this  information  a  revised  set  of  self-consistent  near 
threshold  inelastic  cross  sections  for  each  rare  gas  studied  was  determined. 

3.2.2  Kinetics  of  Interhalogen  Piatomic3 

The  spectra  of  the  S  ■+■  A^Hj  states  of  I?  and  Ill  were  obtained  and  for 
the  former  showed  clearly  resolvable  rotational  structure.  The  lifetimes  of 
the  S  states  of  I?  and  IC1  were  also  estimated  from  pulsed  discharge 

experiments.  The  vibrational  self  quenching  rate  of  the  IE  1  T,  state  was 

estimated  and  also  the  argon  quenching  rate  of  this  same  state  was  measured. 

3.2.3  Multiphoton  Interactions 

In  an  unrelated  study,  we  have  been  exploring  some  effects  associated 
with  the  multi-photon  nature  of  near  resonant  light-matter  interactions 
involving  more  than  one  field  and  a  sat  of  coupled  energy  levels.  Ye  have 
developed  both  algebraic  algorithms  and  a  graph-algebraic  approach  for  solving 
this  class  of  problems  in  the  adiabatic  limit.  These  techniques  were  then 
applied  to  the  specific  case  of  resonant  second  harmonic  generation  in 
three-level  systems,  an  example  of  which  occurs  in  magnetic  resonant  cases. 
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t.  tJa'ititi  sLscTjmns 

Thi3  study  sho ws  that  it  13  possible  to  have  1  vary  high  conversion,  great?? 
than  >■)  1,  of  pu.xp  energy  into  the  3ec0.nl  harmonic  under  selected,  bat  not  too 
restrictive,  coalitions.  The  results  of  this  theoretical  work  will  appear  in 
a  forthcoming  publication.  ^e  are  now  applying  the  finlings  of  this  3tudy  to 
specific  systems  in  orier  to  determine  nhe  feasibility  of  future  experiments. 

5 . 3  Txcited  Itate  Tnemistry  in  lases 

The  research  project  on  Txcited  State  Than is try  in  lases  has  grown  luring 
the  past  year  frox  the  initial  lischarge  worse  on  plasxa  etching  to  include 
pla3xa  annealing  ani  growth  'by  electron  impact"!  ani  xultiphoton  dissociation 
of  sexi  coni  uc  to  r  bearing  gases  30  as  to  grow  thin  filxs  of  Silicon  or 
lermanium  'll 


3.?.1  °lasma  Etching 

The  previous  report  xentionel  the  evidence  for  the  role  of  negative  ions 
in  the  etching  of  3i  or  3il?  with  discharges  in  11?^.  Ve  had  found  that  the 
etch  rate  of  these  xaterials  was  greatly  enhanced  -when  the  3axple  was  biased 
to  a  few  tens  of  volts  positive  with  respect  to  floating  potential.  /e  have 
since  correlated  this  behavior  to  Langmuir  probe  characteristics  which 
indicate  a  massive  negatively  charged  carrier.  7nis  behavior  is  only  evident 
for  discharges  in  '17^;  discharges  in  the  xore  conventional  mixtures  of 
'  ^4.:  5-j )  in  not  exhibit  this  behavior.  This  is  consistent  with  the  fact  that 
T?g  ha3  a  xuch  larger  att achxent  cross  section  than  the  other  gases  V  . 

r  „  n 

These  results  have  been  accepted  fo-  publication  c.o;.  'Jnfortunately,  the 
identity  of  the  ion  is  still  unknown  ani  await.s  the  installation  of  the  xass 

see ctroxeter  on  the  etching  system.  7ni3  is  scheduled  for  coxpletion  in  the 

forthcoxing  year. 

’.%2  °lasxa  Inneallng 

Ii  this  phase  of  she  work,  we  have  shown  that  one  can  use  a  ?ixp!  ? 

abnorxal  discharge  in  a  low  pressure  gas  'typically  heli’im  at  about  ~t  mtor~N 
to  obtain  a  large  area  electron  beam  which,  in  turn,  can  be  used  to  anneal 

large  areas  of  ion  implanted  semiconductor  materials  l'1.  It  has  been  ■■  ir 

experience  chat  the  samples  annealed  in  this  manner  for  times  less  than  ‘  ~ 

seconds  have  sheet  resistivities  consistently  1 es3  than  those  annealed  in  an 
oven  for  cl  minutes.  It  the  present  time  we  are  attempting  to  infer  the 
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temperature  of  the  simple  by  measuring  the  change  in  optical  reflectivity  as 
the  sample  is  irradiated  by  the  electrons.  This  is  a  variation  on  the 

technique  described  by  Olson,  Xikorowski,  Roth  and  Hess  ^R2]. 

"5.3.3  Laser  Thin  Film  Growth 

The  objective  of  this  research  program  is  to  investigate  the  chemistry 
and  optical  properties  of  excited  atoms  and  molecular  radicals  produced  by 

ultraviolet  photo  dissociation  and  the  interaction  of  these  excited  species 
with  a  surface  to  generate  a  thin  film.  Specifically,  an  ultraviolet  excimer 
laser  i3  being  used  to  photodissociate  inorganic  molecules  that  contain  a 
semiconductor  atom  (such  as  3iH^,  CeH^  or  SiCl^) .  The  spectroscopy  of  these 
laser-generated  plasmas  as  veil  as  the  electrical  and  chemical  properties  of 
the  resulting  Ce  or  Si  films  are  under  investigation. 

To  date  amorphous  polycry3talline  3i  and  -3e  thin  films,  around  1  cm*  in 
area,  have  been  grown  by  photolytic  laser  chemical  vapor  deposition  'LCVD)  at 
room  temperature.  Dilute  SiH^,  (CHgj^Si,  or  SeH^  mixtures  in  an  ine^t  carrier 
gas  were  pho todissociated  using  ArF  (193  nm)  or  XrF  (24-3  vm)  lasers.  Film 
thicknesses  of  up  to  0.5  um  were  obtained  on  SiO-,  substrates  at  deposition 
rates  ranging  from  i  to  10  A/sec  for  active  gas  partial  pressures  between  5 

and  30  torr  and  incident  laser  intensities  up  to  25  MW/ cm*.  Absorption 

coefficients  of  as-deposited  3e  and  3i  layers  were  2.3  x  10^  cm'*  and  c.'3  x 
10*  cm-1  ,  respectively,  and  sheet  resistivities  _>  !  0^  ft  /  □  have  been  measured. 
Incubation  times  and  initial  growth  rates  have  been  determined  _in  situ  as  a 
function  of  ’JV  photon  flux  from  He-Ne  laser  absorption  measurements.  A1 
concentrations  of  up  to  several  atomic  percent  have  been  incorporated  into  3e 
films  by  the  simultaneous  photodissociation  of  3a  H^  and  (CH^)^Al. 

Finally,  the  chemical  characteristics  of  these  films  are  being  studied  by 
infrared  absorption  Auger  spectroscopy,  and  T-ray  scattering. 

References 
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of  Laser  Induced  Solid-^hase  Crystal  3rovth  by  Time-Rasolved  Dptical 
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i. 1  Introduction 

This  research  involves  study  of  basic  properties  of  semiconductors, 
methods  of  device  orocessing,  and  new  device  concepts,  3oth  theoretical  and 
experimental  methods  ars  employed  in  each  of  these  categories.  We  are 
examining  a  variety  of  hot  electron  phenomena  and  their  effects  on  oresent  and 
future  device  performance,  especially  in  connection  *ith  modulation  doping. 
We  are  studying  ion  implantation  and  annealing  of  -Si  and  III -7  compounds, 
including  laser  and  electron  beam  processing.  These  experimental  studies 
include  examination  of  deep-level  impurities  and  defects  arising  ?-ot. 
implantation  ani  annealing.  Several  aspects  of  this  ■form  a~e  done  in 
collaboration  with  other  units,  particularly  the  studies  of  materials  gro  vn  by 
molecular  beam  epitaxy  'W33,  'Jnit  t  ) . 
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4 • 2  Laser  and  Electron  3eam  Annealing  of  Semiconductors* 

We  have  performed  both  theoretical  and  experimental  w ork  on  pulsed  laser 
annealing  of  amorphized  Si  and  swept  line  electron  beam  annealing  of  Si  and 
CaAs.  In  collaboration  -with  T.  T.  Verdeyen,  we  have  also  employed  gas 
discharge  annealing  to  process  implanted  Si. 

In  studying  pulsed  laser  melt  annealing  of  Si,  we  have  examined  heating 
by  band-to-band  excitation  (Ndrglass,  1.05  pm)  and  free  carrier  absorption 
(-9?,  10.5  ym)  .  The  technique  of  finite  difference  equations  was  used  to 
solve  the  one -dimensional  heat  conduction  equations  for  the  solid  and  liquid 
phases,  taking  into  account  the  temperature  dependence  of  all  the  parameters 
involved,  especially  the  substrate  parameters,  to  give  a  realistic  model.  In 
the  JQj  case,  we  have  calculated  the  temperature  variation  of  free  carrier 
absorption  in  n-3i.  'Ve  have  taken  into  account  acoustic  deformation  potential 
scattering,  optical  deformation  potential  scattering,  and  ionized  impurity 
scattering.  The  deformation  potentials  were  adjusted  to  fit  the 
experimentally  observed  values  at  300° '<•  Ye  have  investigated  the  "thermal 
runaway"  behavior  and  calculated  the  threshold  power  density  to  cause  thermal 
runaway  and  surface  melting  of  3i  during  C0o  laser  annealing.  For  the  pulsed 
Nd:glass  laser  annealing,  the  melting  model  is  supported  by  time- resolved 
reflectivity  measurements  using  a  He-We  and  an  Ar  laser.  We  used  3EH 
channeling  patterns  to  study  the  quality  and  orientation  of  crystalline 
regrowth.  To  investigate  the  electrical  activation  we  determined  the  change 
in  3heet  resistance  of  the  implanted  material  before  and  after  a  laser  pul3e. 
In  each  case  we  found  excellent  reordering  of  the  amorphized  3i. 

Successful  recrystallization  of  amorphous  BF,*  ion- implanted  3i  has  also 
been  demonstrated  by  the  swept  line  electron  beam  f SLE3)  technique.  High 
percentages  of  electrically  active  boron  were  observed  with  a  minimum  of 
spatial  redistribution.  Active  boron  wa3  also  found  in  the  amorphous- 
crystalline  transition  region,  an  effect  not  observed  in  low  temperature 
furnace  annealed  material.  We  are  currently  employing  the  BLEB  method  in 
recrystallization  of  deposited  3i  layers  over  3iD?. 

*  Thi3  work  was  supported  by  the  Joint  Services  Electronics  Program  'U.S. 
Army,  U.3.  Navy,  and  rJ.3.  Air  Force)  under  Contract  N00D1  d-7)-C-D424,  by  the 
Army  Research  office  under  Contract  DA  AC  29-BO-C-DO1 1 ,  and  by  the  Office  of 
Naval  Research  under  Conract  N0001 4-7S-C-D306. 
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1.3  Studies  of  Transport  ° roperties  in  3 emi conductors* 

Me  have  investigated  carrier  transport  in  semiconductor  hetero.j uno tion 
layers  and  developed  the  theory  of  mobility  in  modulation-doped  structures. 
The  results  of  this  theory  vers  compared  with  experiments  end  found  to  give 
excellent  agreement  i.n  certain  temperature  and  doping  ranges.  Me  have 
continued  to  work  on  the  concept  of  real  space  transfer  and  have  shown 
experimentally  the  validity  of  this  concept  'the  spill  out  of  electrons  from 
potential  wells) .  Me  believe  that  thi3  effect  has  high  device  potential  and 
also  many  interesting  basic  consequences  for  the  transport  in  layered 
structures  and  superlattices.  The  effect  can  also  be  generalised  as  a  real 
3?ace  analog  to  a  k-space  effect  (the  Ounn  effect^  and  thi3  generalization 
leads  to  several  other  phenomena  such  as  "resonant"  impact  ionization  because 
of  the  band  edge  discontinuity  in  hetero.yunction  layers. 

Me  have  expanded  our  Monte  Tarlo  program  'including  the  band  structure  as 
calculated  from  the  empirical  pseudo potential  method)  to  include  t~nnsport  in 
silicon  and  3a  A3.  Ve  were  able  to  calculate  ionization  rates  in  these 
materials  without  too  many  assumptions.  These  calculations  include  'for  the 
first  time)  the  population  of  more  than  one  conduction  band.  Me  have  shown 
that  che  higher  conduction  band  is  crucial  in  silicon  at  high  electric  fields 
a3  they  are  currently  reached  in  small  devices. 

Me  have  continued  to  develop  numerical  programs  to  model  transnort 
phenomena  in  various  semiconductor  devices  and  are  currently  including 
transient  phenomena  { overshoo  t  ^ undershoot  of  velocity''  in  the  operation  of 
charge  couolei  devices.  Part  of  our  theoretical  work  was  devoted  to  the 
possibilities  of  ballistic  transport  and  to  the  field  dependence  of 
heterolayers  at  low  electric  fields. 


*  fhi3  work  wa3  supported  by  the  Toint  lervices  Electronics  °rogram  U.3. 
Army,  '1.3.  Navy,  and"  'J.3.  Anr  Force)  under  Contract  '0331  1.  by  the 

3ffice  of  Naval  Research  under  Tontracts  N3031  V-7 A- 3-333A  and 
U333' l-33-I-3*b3,  and  by  the  Army  Research  3ffice  under  Tontract  TAXI 
a  a- 30-0 -001 1  . 
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4 • 4  Materials  Studies  of  InP* 

As  a  part  of  our  continuing  studies  of  ion  implantation  in  compound 
semiconductors,  we  have  studied  InP  carefully  during  this  period.  The 
redistribution  of  the  compensating  dopants,  iron  or  chromium,  in  semi- 
insulating  indium  phosphide  has  been  studied  using  secondary  ion  mass 
spectrometry  (3IMS).  Annealing  with  silicon  nitride  encapsulation  results  in 
impurity  accumulation  within  the  first  1000  A  of  the  surface  followed  by 
depletion  extending  to  a  depth  of  ~  t  ym.  Profiles  resulting  from  the 
implantation  of  "neutral"  elements  (He,  3)  exhibit  accumulation  at  the  surface 
and  also  accumulation  at  the  projected  range  peak.  The  profiles  can  be 
explained  in  terms  of  gettering  of  the  compensating  dopant  to  defect- rich 
regions . 

We  have  examined  polycrystalline  and  single  crystal,  liquid  encapsulated 
Csochralski  (L3C)  grown  InP  for  evidence  of  boron  contamination  using  3 IMS  and 
photol'uminescence .  Precipitates  of  boron  or  a  boron  compound  have  been  found 
in  In?  grown  by  the  LSC  method  'with  boric  oxide  encapsulation  and 

pyrolytic  boron  nitride  (pBN)  crucibles.  The  density  of  precipitates  appears 
to  be  solely  responsible  for  the  observed  phenomenon.  At  present,  the  source 
of  boron  contamination  most  consistent  -with  our  observations  is  the  result  of 
interaction  between  3,0^  and  pBN.  A  diffusion  coefficient  of  D|'750°C)  <  1  x 
1  0_1  ^  cm-'  s"'  for  boron  in  InP  has  been  determined  in  this  work. 

We  have  observed  cracking  in  chemical  vapor-deposited  3i0o  encapsulating 
layers  on  In?  when  these  samples  were  annealed  at  temperatures  above  c50c0. 
Using  optical  microscopy,  scanning  electron  microscopy  and  Auger  electron 
spectroscopy  we  have  studied  the  detailed  nature  of  this  form  of  film  failure. 
Thermally  induced  stress  between  th9  film  and  the  InP  does  not  fully  account 
for  the  phenomena  observed.  The  cracking  apparently  results  from  an 
interaction  between  mechanical  stress,  chemical  effects  3nd  possibly  defects 
in  the  SiO,  layers. 


*  Thi3  work  -was  supported  by  the  Joint  Services  Slectronics  Program  (7.3. 
Army.  7.3.  Wavy,  and  7.3.  Air  Force)  under  Contract  70001  4-70-C-O 424,  and  by 
the  Office  of  Waval  Research  under  Contract  70001 i-?6-C-030o. 
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For  ion  implantations  of  In?  that  are  highly  dominated  by  electronic 
stopping,  the  standard  Cau3sian  approximation  provides  a  poor  fit  to  the 

observe!  experimental  profile.  Using  four  moments  base!  on  theoretical 
considerations,  we  have  constructed  profiles  using  a  terminated  Edgeworth  and 
Pearson  I  distribution.  It  was  found  that  neither  the  Edgeworth  nor  Pearson  I 
provides  the  accuracy  that  would  be  expected  from  a  four  moment  distribution. 
The  moments  of  the  experimental  profiles  have  been  empirically  determined. 
These  moments  together  with  the  Pearson  IV  distribution  aere  shown  to  fit 

experimental  profiles  over  at  least  two  orders  of  magnitude. 

With  the  background  described  above,  we  are  currently  studying  the 
effects  of  3e,  Mg,  and  Si  implants  in  InP.  These  studies  include  impurity 
profiling,  electrical  activation,  and 
impurities . 


interaction  with  defects  and  background 
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5 . 1  Semiconductor  Crystal  Srowth  from  the  7 apor  Phase: 
Ion-Surface,  Plasma,  and  Laser  Stimulated  Reactions 


We  are  studying  energetic  particle-surface  interactions  and  stimulated 

gas  phase  reactions  which  control  the  nucleation  and  growth  kinetics, 
chemistry,  and  physical  properties  of  compound  and  alloy  semiconductors  grown 
from  the  vapor  phase  by  UHV  ion  beam  sputtering,  plasma- assisted  chemical 
vapor  deposition,  and  laser- induced  chemical  vapor  deposition.  The  common 
feature  of  these  techniques  is  that  crystal  growth  proceeds  under 
non- equilibrium  thermodynamic  conditions  through  the  production  of  highly 
reactive  gas  phase  species:  excited  atoms,  metastables,  radicals,  and  ions. 
Such  species  transfer  energy  to  the  growth  surface  upon  condensation  thereby 
altering  the  surface  reactivity  as  well  as  adsorption  and  adatom  diffusion 
kinetics  allowing  film  growth  at  lower  temperatures,  a  wider  range  in 
controlling  doping  concentrations  and  tailoring  film  properties,  and  the 
growth  of  unique  metastable  materials.  This  work  is  being  pursued  from  both 
an  analytical  and  an  experimental  point  of  view  to  establish  a  detailed 
understanding  of  fundamental  film  growth  mechanisms.  We  have  recently 
published  two  invited  review  papers  ^5.1,  5»6]  in  this  area.  Results  from 
this  research  have  a  wide  range  of  applications  in  addition  to  crystal  growth 
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including  the  active  research  fields  of  reactive  ion  etching,  ion  bean 
lithography,  nicrochemical  analysis,  plasma  chemistry,  and  laser-materials 
interactions . 

5. 1 .  1  1  rowth  of  high  Quality  Single  1  rystal  Sal 3  Films 

by  Sputter  Deposition 

epitaxial  SaAs  films  have  been  grown  on  ( 1 00 )  Tr- doped  SaAs  substrates  by 
3nutter  deposition  r5.2"!  from  an  unioned  SaA3  target.  Unintentional! y  iopei 

samples  were  high  resistivity  105  to  *  ft-cm  n-type  -with  room  temperature 

mobilities  as  high  as  5 000  cm-/7-sec.  Phis  represents  one  of  the  highest 
reported  mobilities  in  high  resistivity  SaA3.  Residual  oxygen  and  carbon 
contamination  was  minimised  through  the  use  of  a  liquid  nitrogen  cooled  shroud 
surrounding  the  discharge  during  deposition.  Yxcess  arsenic  was  provided  to 
the  growing  film  from  an  effusion  cell.  Film  doping  wa3  accomplished  using 
either  an  evaporated  source  (jin')  or  the  addition  of  a  gas  phase  impurity  ’Rog'i 
during  snutter  deposition.  ’in-doped  p-type  films  were  grown  with  room 

temperature  carrier  concentrations  between  1.7  x  ID17  and  2.5  x  ID15  cm-'1  ani 
corresponding  hole  mobilities  of  24-0  to  cm*  /’/-sec.  These  values  agree  very 
well  -.with  results  for  both  LPT  ani  733  films  doped  to  similar  concentrations . 
The  0  doped  films  were  n-type  (see  section  5.1.2).  In  all  cases,  the  dopant 
incorporation  probability  was  controllable  by  varying  the  negative  bias  on  the 
growing  film. 

5.1.2  I  on  Bombardment  Tf facts  on  3 lemental  Sticking  Probabilities 

We  have  shorn  previously  that  low  energy  (<  700  eV)  ion  bombardment  of  a 
growing  film  can  play  a  large  role  in  controlling  elemental  sticking 


probabilities,  allowing  an 

extension 

in 

the  growth 

temperature 

range  over 

which  stoichiometric  III  -7 

compounds 

and 

alloys  can 

be  formed 

as  veil  as 

allowing  the  growth  of  non- 

equilibrium 

phas 

es  ( 3ee  sec 

tion  5.2'). 

^his  occurs 

through  a  variety  of  mechanisms  including  self- ion  trapping,  enhanced  surface 
reactivity,  sputter  removal  of  adatoms,  and  impact-stimulated  dissociative 
chemisorption.  For  a  discussion  of  these  mechanisms  and  related  equations  see 
the  review  article  by  Ireene  _5.71, 
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Sulfur  concentrations  of  up  to  1.2  x  1 021  cm“^ ,  as  measured  by  secondary 
ion  mass  3pectroscopy,  were  obtained  in  sputter  deposited  single  crystal  CaAs 
layers  using  an  H93  glow  discharge  source.  The  incorporated  3  concentration 
increased  linearly  with  increasing  H?s  partial  pressure  but  was  not  affected 
by  varying  the  As /3a  flux  ratio  between  12  and  4-0.  The  S  incorporation 
probability  varied  from  0.02  for  essentially  thermal  species  to  nearly  ’unity 
for  sulfur  containing  ions  accelerated  by  a  negative  substrate  potential  of 
325  V.  The  enhancement  in  S  incorporation  probability  was  not  due  only  to 
implantation  processes,  but  was  also  related  to  the  strength  of  the  Ca-3  bond. 
Very  3harp  doping  profiles  with  no  indication  of  surface  segregation  were 
obtained  for  programmed  doping  steps  (for  comparison,  see  results  for  Sn 
segregation  in  MBE  OaAs,  from  a  study  carried  out  in  collaboration  'with 
Prof.  Morkop,  summarized  in  Unit  l). 

5 . 2  1 rowth  of  Single  C rystal  Metastable  Semiconductors* 

As  reported  last  year  [31,  52],  we  have  carried  out  the  first  detailed 
study  of  the  growth  of  new  single  crystal  metastable  semiconductors.  The  key 
feature  in  stabilizing  the  growth  of  these  materials  was  the  controlled  use  of 
low  energy  ion  bombardment  during  deposition  to  modify  elemental  sticking 
probabilities  and  to  enhance  adatom  diffusivities .  During  the  past  year,  we 
have  investigated  the  crystal  growth,  thermal  and  temporal  stability, 
electrical  properties,  and  the  thermodynamics  of  a  new  subclass  of  metastable 
semiconductors,  (tII-V)1 _x(lV)x. 

Epitaxial  metastable  (3aAs)1_xSex  alloys  with  compositions  between  x  = 
0.10  and  0.75  have  been  grown  on  (100)  laAs  substrates  by  ion  beam  sputtering 
in  an  ultrahigh  vacuum  system.  Electron  channeling,  double  crystal  X-ray 
diffrac tometry ,  and  X-ray  topography  analyses  indi^r^e  that  the  films  are  of 
very  high  crystalline  perfection.  Layers  with  x-values  up  to  at  least  x  = 
0.4?  were  in  the  zinc  blende  structure.  Either  n-type  or  p-type  conduction 
•with  n  and  p  varying  over  several  orders  of  magnitude  could  be  obtained  by 
varying  the  film  composition,  the  growth  temperature,  and  the  A3  overpressure 
during  deposition.  The  equilibrium  0aA3-0e  pseudobinary  phase  diagram  has 

*  Thi3  research  was  supported  by  the  Joint  Services  Electronics  Program 
(’J.3.  Army,  TJ.S.  Navy,  and  U.S.  Air  Force-)  under  Contract  N0001  4-79-C -0424, 
and  by  the  Department  of  Energy. 
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been  determined . 

Recent  thermodynamic  calculations  combined  with  differential  scanning 
calorimetry  measurements  have  provided  a  more  detailed  understanding  of  both 
the  growth  and  the  thermal  stability  of  these  unique  materials.  EXAFS 

measurements,  carried  out  in  collaboration  with  Prof.  Then  of  the  Metallurgy 
Department  at  HINT,  combined  -with  STEM  and  thermal  annealing  results  have 

shorn  conclusively  that  eollisional  mixing  due  to  low  energy  ion  bombardment 
can  be  used  to  retard  phase  transformations.  finally,  we  have  begun 
collaborations  -.nth  Prof.  Klein  of  the  Physics  Department  at  7173  and 

Prof.  Raccah  of  the  Physics  Department  at  7 ITT  to  investigate  optical 
properties  o:  metastaole  '  la  As''  ^  _x3ex  and  [  3a3b)  i  _x3e.-  alloys  using  Raman 
soectroscopv ,  electro- reflectance ,  and  photo  luminescence . 

5 . ’  Reactive  Ion  Etching  of  3 aA s* 

The  first  detailed  study  of  tne  reactive  ion  etching  of  TaA3  was 

carried  out  by  our  group  during  the  past  year.  Df  primary  interest  wa3  the 
effect  of  ion-surface  interactions  on  the  mechanisms  and  kinetics  of  surface 
reactions.  Experiments  were  carried  out  in  pure  TT1,?,  ry  discharges,  as 
w>ll  as  in  dilute  CCl^,  TCloFo,  and  T?^  discharges  containing  ID  model  Ar. 
Anisotropic  etching  -with  removal  rates,  R,  of  up  to  D.3  ym/ min  have  been 

obtained  in  reactive  discharges  operated  at  10  mTorr  and  -3  kV,  whereas  the 
physical  sputtering  rate  in  pure  Ar  discharges  operated  under  the  same 
conditions  was  only  D.  D1  ym/min.  A  combination  of  optical  emission  and 
absorption  spectroscopies  have  been  used  to  show  that  in  both  pure  and  dilute 
chlorine-containing  reactive  discharges,  physical  sputtering  of  atomic  Ta  and 
A3  is  not  the  primary  etching  mechanism  for  Da As  although  ion  bombardment  does 
play  an  important  role  in  the  overall  process.  Rather,  etching  occurs  by  the 
formation  and  desorption  of  volatile  3a  and  As  chlorides.  In  reactive 
fluorine-containing  discharges,  the  rate  limiting  3tep  is  the  removal  of 
non-volatile  3aFx. 


*  Thi3  research  wa3  supported  by  the  Joint  Services  Jl:'tronics  Program 
,J.3.  Army.  7.3.  Navy(  and  7.3.  Air  Force)  -under  Ton  tree  .i'.DTl  l-’o-T-DITi, 
and  by  the  I3M  Researcn  Laboratory,  East  Fishkill,  New  York. 
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5.4  Laser-Surface  Interactions 

In  collaboration  with  Prof.  Eden  of  the  Electrical  Engineering  Department 
at  UIUC  we  have  initiated  a  program  to  investigate  film  growth  by  laser 
induced  chemical  vapor  deposition  (LCVD).  In  initial  experiments  we  have 
studied  the  growth  of  amorphous  hydrogenated  Si  and  Ge  thin  films,  ~  1  cm2  in 
area,  by  photolysis.  Dilute  SiH^,  (CH-^)^Si,  or  GeH^  mixtures  in  an  inert 
carrier  gas  were  photodissociated  using  ArF  (193  nm)  or  KrF  (243  nm)  lasers. 
Film  thicknesses  of  up  to  0.5  pm  were  obtained  on  3iC>2  substrates  at 
deposition  rates  ranging  from  4  to  1 0  A/sec  for  active  gas  partial  pressures 
between  5  and  30  Torr  and  incident  laser  intensities  up  to  25  NW/cm-. 
Absorption  coefficients  of  as-deposited  Oe:H  and  Si:H  layers  were  2.3  x  1 05 
cm-1  and  5*9  x  10“^  cm"',  respectively,  and  sheet  resistivities  _>  10^  ft/o  have 
been  measured.  Incubation  times  and  initial  growth  rates  have  been  determined 
in-si tu  as  a  function  of  OT  photon  flux  from  He-Ne  laser  absorption 
measurements.  A1  concentrations  of  up  to  several  atomic  percent  have  been 
incorporated  into  Ge  films  by  the  simultaneous  photodissociation  of  GeH^  and 

A  parallel  effort,  carried  out  in  collaboration  with  the  Physics  Division 
of  Eastman  Kodak  Research  Laboratories,  Rochester,  New  York,  has  been 
established  to  investigate  liquid  phase  regrowth  of  semiconductors  by  OW  laser 
annealing.  We  have  recently  demonstrated  the  growth  of  epitaxial  Ge/GaAs 
heterostructures  by  scanned  CW  argon  laser  annealing  of  440  nm  thick  amorphous 
Ge  films  on  (100)  3emi- insulating  GaAs  substrates  [5.3].  Depending  on  the 
incident  laser  power  and  scan  rate,  two  modes  of  film  regrowth  ware  observed. 
At  low  powers  (between  a,  1.6  and  4.0  'tf  for  a  beam  diameter  of  a,  40  um)  and 
scan  rates  between  1  and  400  cm/s,  polycrystalline  Ge  with  a  (100)  preferred 
orientation  was  formed  by  an  "explosive”  crystallization  mechanism.  At  higher 
powers,  and  over  a  scan  rate  range  of  20  to  400  cm/s,  single-crystal  films 
containing  some  dissolved  GaAs  in  solution  were  obtained  by  liquid-phase 
regrowth.  Typical  film  resistivities  p  were  as  follows:  as-deposited,  130 
Q-cm;  polycrystalline  films,  3  x  1 0“-  ft-cm;  single- crystal  films,  9  x  1 0"'1 
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6. 1  Introduction 

The  objective  of  the  microwave  acoustic  research  is  five  fold: 

1 .  to  use  the  interactions  between  acoustic  waves  and  free  carriers  in 
CaAs  to  make  charge  transport  devices, 

2.  to  use  traveling  wave  electric  fields  similar  to  those  produced  by 
acoustic  waves  to  evaluate  semiconductor  electronic  properties, 

to  integrate  semiconductor  active  devices  and  microwave  acoustic 
devices , 

4.  to  analyze  surface  acoustic  wave  transducers  and  reflectors  from 
basic  material  parameters,  and 

5.  to  investigate  mode  conversion  techniques  of  line  acoustic  wave 
generation. 


S.  Mahon 

R.  Miller 

S.  Vilkus 


Progress  toward  this  objective  is  described  in  the  following  paragraphs. 
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5.2  Electron  Transport  in  Piezoelectric  Semiconductors* 

Acoustic  surface  and  line  waves  propagating  in  piezoelectric 

semiconductors  drag  the  carriers  and  produce  an  acoustoelectric  current.  'When 
the  electric  fields  oroduced  by  the  waves  are  greater  than  the  30und  wave 
velocity  divided  by  the  mobility,  the  carriers  are  dragged  at  the  sound 
velocity.  Further  increases  in  acoustic  power  do  not  increase  velocity  or  the 
acoustoelectric  current.  The  carriers  are  then  3aid  to  propagate 
synchronously  and  the  acoustoelectric  current  is  saturated. 

The  charge  transport  devices  envisioned  in  this  research  depend 
exclusively  on  synchronous  carrier  transport.  Extensive  study  has  been 
carried  out  to  determine  the  nature  of  charge  transport  in  the  presence  of 
very  large  traveling  wave  fields.  The  effect  of  traps  on  synchronous  carrier 
transport  is  of  particular  concern  and  Cr  doped  UaA3  has  been  studied 

extensively  because  of  its  high  density  of  traps. 

The  saturated  acoustoelectric  current,  caused  by  the  nonlinear 
interaction  of  photo-generated  electrons  -with  large  amplitude  SAW  induced 
electric  fields  in  Cr  doped  UaAs,  i3  found  to  be  too  small  to  be  explained  on 
the  premise  that  all  electrons  move  synchronously  at  the  SAW  velocity.  Thi3 
finding  is  based  on  the  results  of  large  signal  acoustoelectric  measurements, 
for  SAW  both  internal  and  external  to  the  CaAs  substrate,  via  a  separated 
media  device  geometry.  The  expected  linear  dependence  of  saturated  current 
with  SAW  velocity  is  not  observed  when  SAW  with  different  velocities  are  used 
to  induce  the  saturation.  It  is  concluded  that  some  mechanism  other  than 
complete  carrier  bunching  is  responsible  for  acoustoelectric  current 

saturation  in  CrrSaAs  [R1 ]. 

During  the  last  year,  we  have  demonstrated  that  electrons  can  be 

entrained  by  an  acoustic  wave  propagating  along  the  surface  of  a  Tr-doped  SaAs 
crystal.  This  result  was  surprising  in  view  of  the  weak  piezoelectricity  of 
SaAs;  it  was  made  possible  by  the  unusual  property  of  high-frequency  acoustic 
surface  waves  of  providing  a  concentrated  field  over  a  narrow  cross  section. 


Electronics  Program  under 


*  This  work  was  sunoorted  by  the  Joint  Services 
conract  N000U-79-S -D‘124. 
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One  can  visualize  a  new  class  of  high-performance  charge  transfer  ievices 
based  on  acoustically  entrained  carriers.  Such  devices  would  have  no  gates 
and  no  complex  multi-layer  patterns  to  interconnect  the  gates;  the  region 
along  which  the  electrons  travel  would  be  free  for  pickup  electrodes.  Along 
with  the  simpler  structure,  these  devices  would  have  negligible  transfer  loss, 
making  it  possible  to  obtain  a  time- band width  product  orders  of  magnitude 
higher  than  that  of  present-day  CCD's. 

6 . 3  A_  New  Method  to  Measure  the  Mobility  of  .Amorphous  Semiconductor  Pilms* 

The  measurement  of  mobility  in  high  resistivity  low  mobility  amorphous 
semiconductor  films  can  be  difficult.  We  have  used  the  acoustoelectric 
current  described  in  the  last  section  as  a  measure  of  material  mobility. 

The  acoustoelectric  method  of  determining  carrier  mobility  u  in 
semiconductor  films  has  been  modified  to  permit  its  application  to 
high-resistance,  low-mobility  films.  The  conventional  method  requires 
knowledge  of  the  acoustic  loss  caused  by  the  mobile  carriers;  in 
high- resistance  films  this  loss  becomes  too  small  to  be  measured.  We  show 
that  the  required  information  may  be  derived  from  knowledge  of  the  acoustic 
power  and  the  device  geometry.  Our  samples  were  amorphous  hydrogenated  Si  and 
3iQ.g(JeQ  ^  films  on  nonpiezoelectric  substrates,  separated  by  a  convenient  air 
gap  (12.5  urn)  from  a  LiNbO-^  slab  carrying  surface  acoustic  waves.  One  sample 
had  1 0®  ft  /□  and  u  =  0.08  cm  - /V  sec;  another  sample  10'^  Q./O  and  u  =*  0.5  cm~/V 
sec  [o.l  ]. 

In  amorphous  materials  the  drift  mobility  of  carriers  is  usually  measured 
by  observing  the  time-of- flight  of  a  charge  packet  through  a  sample  in  an 
external  d.c.  electric  field.  The  results,  however,  cannot  be  interpreted  in 
terms  of  a  carrier  mobility  in  the  conventional  sense  due  to  the  spreading  of 
the  packet  by  multiple  trapping  and  release,  and  the  concept  of  time-dependent 
mobility  has  been  used  to  describe  such  dispersive  transport.  Recently  the 
mobility  was  also  deduced  from  transient  photoconductivity  and  photo-induced 
absorption  experiments.  U3ing  the  acou3tic-elec trie  method  to  measure 
mobility  provides  new  information.  \  simple  coherent  view  showing  the 
connection  among  all  these  different  measurement  techniques  ha3  been  written 

*  This  work  was  supoortei  by  the  Joint  Services  Electronics  Program  under 
contract  N000I  4 -79-1-3124. 


■iiiiai— - . 


MICROWAVE  ACOUSTICS 


o. 

[32]. 


6 . 4  Wideband  LAW  Transducer* 

The  fundamental  advantages  of  LAW  devices  are  veil  documented  but  there 
has  been  no  efficient  widebandwiith  transducer  vith  which  to  implement  them. 
We  have  recently  devised  a  mode  converter  that  converts  surface  acoustic  wave 
SAW  power  to  LAW  power.  This  device  consists  of  a  number  of  unconnected 
electroles  which  extend  across  the  SAW  and  LAW  propagation  paths.  The 
electrodes  are  bent  and  spaced  so  that  the  electrode  spacing,  measured  in 
phase  at  center  frequency,  is  the  same  for  the  LAW  and  the  SAW  beams.  The 
first  generation  device  had  a  o"o  bandwidth  and  a  loss  of  5  13  [o.Ml. 

The  second  generation  has  a  bandwidth  increased  to  10^  and  a  loss  reduced 
to  1  d3.  With  this  band-width  and  efficiency,  LAW  devices  may  become  practical 
because  efficient  SAW  transducers  can  be  U3ed  to  generate  the  SAW  beam  which 
in  turn  is  converted  to  the  LAW  beam  by  the  Multi3trip  Mode  Converter  (MSMC). 
The  mode  converter  is  reciprocal  and  LAW  waves  can  be  detected  by  conversion 
to  SAW,  and  then  a  SAW  transducer  is  used  to  ietect  them.  Efficient  moderate 
bandwidth  LAW  delay  lines  and  convolvers  are  now  possible. 

5 . 5  A_  Theoretical  Determination  of  Surface  Acoustic  Wave  Velocity 

and  Impedance  Differences  between  Metal  Strips  and  Free  Surface 

Regions  of  Metallic  P ratings** 

The  surface  acoustic  wave  i3  a  3ingle-moded ,  non- dispersive , 
one- dimensional  wave  that  is  frequently  modeled  using  a  transmission  line 
equivalent  circuit.  Thin  electrodes  and  shallow  grooves  located  in  the 
surface  acoustic  wave  propagation  path  produce  reflections  and  shift  the 
propagation  velocity  of  the  wave.  The  effect  of  the  electrodes  and  grooves 
are  modeled  a3  an  offset  in  characteristic  impedance  'DELTERE)  and  a  shift  in 
the  propagation  velocity  (DSLTSRV).  Thi3  impedance  offset  and  the  velocity 
shift  have  been  theoretically  determined  from  the  material  constants  and  the 

*  Thi3  work  was  sunoorted  by  the  Rome  Air  Development  lenter,  Deouty  for 
Electronic  Technology,  under  Air  Force  contract  F-l  >> 23-31  -P-DP?t  . 

**  This  work  was  3uooortei  bv  the  Rome  Air  Develooment  lantar,  Deouty  for 
Electronic  Technology,'  under  Air  Force  contract  F-l  AoT 3-7E-  1-DD1D. 
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geometrical  dimensions.  Since  the  materials  are  piezoelectric,  both  the 
elastic  and  electrical  loading  of  the  surface  have  been  considered.  The 
detailed  theoretical  derivations  for  the  impedance  offsets  and  velocity  shifts 
have  been  published  in  a  series  of  papers  during  the  course  of  the  development 
[5. 2, 6. 4-6].  A  report  has  been  written  to  tie  these  studies  together  and  to 
provide  a  simplified  analysis  of  the  surface  loading  produced  by  an  electrode 
that  will  be  useful  to  the  device  designer.  The  report  has  been  made  concise, 
with  all  derivations  taken  from  references  or  appendices  in  order  to  make  it  a 
better  design  tool. 
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7 . 1  Dif fusion  and  Reactivity  on  Met a la4 

The  decomposition  of  polyatomic  gases  on  30UI3  i3  of  interest  in  both 
chemical  a.ni  al ectronic3  technology.  We  have  been  examining  tha  interne tions 
of  methane  with  metal3  as  a  simple  but  manageable  example  of  auch  processes. 
Thi3  effort  has  been  concentrate!  upon  understani ing  the  mechanisms  by  which 
l  scorn  position  occurs;  of  particular  interest  is  -/ay  there  is  an  activation 
enargy  to  iecompo si tio 1  at  the  surface  and  how  thi3  i3  overcome. 

Measurements  reports!  last  year  ^Rt]  indicated  that  the  activation  energy, 
letarminei  when  the  temperature  of  both  the  surface  and  the  gases  i3  change!, 
is  considerably  smaller  than  the  value  found  some  time  ago  R2  !  in  beam 
experiments,  in  which  only  the  gas  i3  excite!.  Tnere  are  two  distinct 
possibilities  to  account  for  this: 


*7hi3  work  was  supports!  by  the  Rational  Science  Foundation  ander,  grant  TTR 
"3-2311  5.  as  well  ’as  by  tne  Joint  3ervices  Electronics  ^-ogram  '7.3.  Army, 
7.3.  'lavy,  U.3.  Ur  Force)  ander  contract  '1310!  -1 12 1. 
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1.  The  beam  experiments  were  done  on  single  crystal  tips,  rather  than  on 
evaporated  films,  which  served  as  substrate  in  the  more  recent  tom. 
There  may  be  a  significant  dependence  of  decomposition  rate  upon 
crystallography,  which  could  account  for  the  differences. 

2.  If  the  reaction  proceeds  through  a  precursor,  that  is  through  a  gas 
molecularly  bound  to  the  surface  prior  to  dissociation,  different 
activation  energies  are  expected  in  the  two  measurements. 

3oth  types  of  effects  are  likely,  and  we  have  carried  out  several  studies 
to  define  more  closely  the  important  events  at  the  surface.  Measurements  it 
low  temperatures  indicate  a  heat  of  physical  adsorption  of  5  kcal/mole.  Under 
the  conditions  of  the  decomposition  experiments,  a  reservoir  of  molecularly 
held  methane  therefore  exists  on  the  surface  and  dissociation  should  take 
place  via  these  molecules.  This  suggests  that  in  a  molecular  beam  experiment 
in  which  the  ga3  alone  i3  heated,  the  activation  energy  3'nould  be  3  kcal'mole 
higher  than  in  stardard  measurements,  and  should  amount  to  roughly  H 
kcal/mole  on  rhodium  film.  To  test  this  hypothesis  a  new  apparatus,  shown  in 
Fig.  7.1,  has  been  built.  In  this  equipment,  the  gas  can  be  heated  while 
keeping  the  film,  which  is  in  contact  with  a  large  thermal  mass,  at  a  constant 
temperature.  The  results  obtained  in  this  way  are  plotted  in  Fig.  t.o. 
After  correcting  for  contributions  from  gas  molecules  that  thermally 
accommodate  to  the  surface,  an  activation  energy  of  ID. 7  kcal/mole  is  found 
for  methane  decomposing  on  rhodium  films. 

This  i3  in  excellent  agreement  -with  expectations;  more  than  that,  it  is 
consistent  .with  other  phenomena  observed  in  the  interactions  of  methane  with 
rhodium.  It  has  previously  been  found  r_R7,R4]  that  exciting  the  v stretching 
mode,  which  has  an  energy  of  3.6  kcal/mole,  is  ineffective  in  accelle-ating 
decomposition.  This  is  now  iuite  -understandable .  An  energy  of  11  kcal/mole, 
much  more  than  i3  available  in  the  vibrations,  must  be  supplied  to  methane 
in  order  for  it  to  dissociate.  Most  likely  it  is  the  5V^  bending  mode  that 
must  be  excited  in  order  to  enhance  the  reaction  on  rhodium.  Detailed 
predictions  of  the  effect  cf  isotopic  substitution  have  been  made,  based  cn 
the  assumption  that  vibrational  excitation  is  important;  these  predictions 
are  in  good  accord  with  experiments. 
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rhodium  films  maintained  at  0®C. 
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The  work  carried  out  her  a  has  clarified  the  role  of  the  gas  in  the 
de-corn position  reaction.  Still  to  be  explored  now  are  the  contributions  of  the 
aolii,  and  the  extent  to  which  these  are  specifically  affected  by  the  atonic 
arrangement  at  the  surface. 

7 . 2  0 if fusion  and  Reactivity  on  Silicon* 

Although  the  importance  of  surface  effects  in  solid  state  devices  is 
generally  recognised,  relatively  little  is  known  about  the  properties  of 
surface  layers  on  semiconductors.  Ve  have  been  concentrating  our  efforts  on 
obtaining  information  on  the  spatial  stability  of  such  layers.  As  part  of 
these  investigations  we  have  developed  a  technique  for  direct  electron  beam 
writing  on  silicon  _7.2].  In  thi3  method,  molecular  nitrogen  is  physically 
adsorbed  on  silicon,  maintained  at  low  temperatures  at  which  a  monolayer  or 
more  i3  formed.  Molecular  nitrogen  on  the  surface  i3  then  dissociated  by 
electron  beam  impact.  In  warming  the  sample,  undissociated  molecular  nitrogen 
evaporates,  leaving  behind  a  nitrogen  layer  wherever  the  electron  beam  ha 3 
3 truck,  as  shown  in  Fig.  "\3 • 

Current  efforts  have  been  focused  upon  characterising  the  conditions 
under  which  these  nitrogen  layers  are  forms!.  Cur  studies  indicate  that 
electron  impact  i3  quite  efficient  in  promoting  dissociation  of  nitrogen. 
Ies3  than  5  x  ! electrons/cm-  suffice  to  create  a  saturation  layer; 
further  bombardment  does  not  increase  the  amount  of  nitrogen  chemically  bound 
to  the  silicon. 

These  saturation  layers  are  actually  nitrides,  not  just  chemisorbed 
layers  of  nitrogen.  It  i3  clear  from  Fig.  7.4-  that  the  Si'LW)  Auger  spectra 
of  the  layers  produced  by  elecron  bombardment  are  in  good  accord  with  the 
spectra  recorded  for  silicon  nitrile  produced  by  chemical  vapor  deposition 
Me  have  been  able  to  show,  however,  that  during  the  early  stages  of  the 
proce33,  what  appears  to  be  a  chemisorbed  layer  of  nitrogen  does  form.  This 
is  illustrated  in  Fig.  7o,  where  the  intensities  of  various  Auger  lines  are 
shown  after  increasing  exposures  of  nitrogen  and  concomitant  bombardment  with 
electrons.  The  Id-eV  line  charao teristi c  of  the  nitrile  does  not  appear  at 

*  This  work  was  supported  bv  the  Toint  "Services  Electronics  °-ogram 
v'J.3.  Army.  '1.1.  'lavy,  u.1.  Air  Force)  under  contract  'TOO  1*  l-7 1- d -"H2 i. 
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Fig.  7.3 


Secondary  electron  images  of  nitride  structures  created  by  electron 
beam  writing  on  molecular  nitrogen  layers  held  on  Si(100). 
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Fig.  7.5  Auger  intensities  observed  after  increasing  nitrogen  exposures 
and  electron  irradiation.  The  83-eV  nitride  line  only  develops 
after  an  induction  period,  during  which  a  chemisorbed  layer  forms 
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the  very  start.  However,  the  330 -eV  line  attribute!  to  nitrogen  is  appareat 
even  at  /ery  low  exposures,  indicating  that  there  is  a  nitrogen  layer  present, 
though  it  i3  not  a  nitrile.  The  nitride  layers  formed  after  extensive 
exposure  are  more  than  one  layer  thick.  The  mechanism  of  formation  of  these 
nitriles  is  not  clear,  but  it  is  quite  complex  and  may  involve  a  competition 
between  electron  impact  desorption  and  dissociation. 

Electron  beam  induced  formation  of  silicon  nitrile  is  a  particularly 
appealing  process.  It  gives  U3  the  capability  of  forming  nitrile  structures 
under  controlled  and  very  clean  conditions.  A  question  immediately  arises  - 
is  this  technique  applicable  to  the  formation  of  other  films  as  well?  We  have 
recently  examined  the  feasibility  of  creating  oxide  layers  on  silicon  in  this 
way.  On  exposing  silicon  to  molecular  oxygen,  a  chemisorbed  layer  forms  on 
the  surface.  At  room  temperature  and  below,  oxidation  i3  slow;  however, 
electron  bombardment  enhances  the  oxidation  rate.  As  i3  apparent  from  Fig. 
7.6.  it  is  therefore  possible  to  form  oxide  structures  by  electron  beam 
irradiation  of  an  oxygen  covered  surface.  7ne  most  advantageous  conditions 
for  this  process  are  now  being  explored. 

7 . 3  1 tomic  Exploration  of  3 rystal  Layers* 

Over  the  years  considerable  information  has  been  accumulated  in  thi3 
laboratory  about  the  behavior  of  individual  metal  atoms  on  metal  surfaces 
^7.l].  Thi3  has  been  possible  through  the  use  of  the  field  ion  microscope  and 
it3  ability  to  resolve  single  atoms.  While  work  on  single  atoms  has 
continued,  we  have  also  been  examining  the  behavior  of  overlayers,  with  the 
intent  of  contributing  to  a  better  understanding  of  crystal  and  film  growth 
phenomena;  in  the  past  these  phenomena  have  not  been  accessible  to 
observation  on  the  atomic  level. 

Initially,  we  have  concentrated  on  a  survey  of  overlayar  formation  on 
W(llO').  Metal3  from  the  platinum  family  form  two-dimensional  layers,  provided 
the  adatom  concentration  on  the  substrate  is  high  enough.  As  is  apparent  in 
the  field  ion  micrographs  in  Fig.  7.7,  individual  atoms  in  the  layers  are 
clearly  resolved.  Thi3  is  quite  different  from  whit  is  found  when  the 

*  Thi3  work  was  supported  by  the  National  Science  Foundation  -under  grant  DHR 
30-11361  and  by  thV  Joint  cervices  Electronics  Program  ('J .  3.  Army,  U.  3.  Navy, 
■J.3.  Air  Force)  under  contract  N003'  l-’3— 3-012A. 
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Fig.  7.6  Secondary  electron  image  of  oxide  structures  formed  by  electron 


beam  writing  on  Si(100). 
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Fig.  7.7  Platinum  metal  layers  on  W(110).  Shown  in  order  are  palladium, 
iridium,  and  platinum  monolayers,  as  well  as  a  schematic  of  a 
p(2  x  2)  structure. 


Fig.  7.8  Field  ion  image  of  silicon  layer  on  (110)  plane  of  tungsten 
(left),  and  proposed  structure  at  right. 
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substrata  itself  is  imaged.  On  a  W('1  3)  plane  of  the  same  size,  individual 
atoms  cannot  be  discerned.  The  spacing  in  the  layers  of  platinum,  palladium, 
ani  iridium  must  therefore  be  a  multiple  of  the  substrate  spacing.  This  is 
not  necessarily  true  for  other  metals;  for  rhenium,  iniiviiual  adatoms  are 
not  resolved  in  overlayers  on  V(!10),  ani  we  therefore  infer  a  much  denser 
packing . 

A  particularly  interesting  behavior  has  been  observed  for  iridium 

clusters.  The ir  arrangement  ani  stability  is  strongly  dependent  upon  size. 

At  450  X  for  example,  a  cluster  of  15  or  more  iridium  atoms  on  W(l  1  0)  is 
stable  in  a  rhomboheiral  configuration.  A  cluster  of  only  1  o  adatoms, 
however,  occurs  in  a  square  arrangement.  This  is  not  an  artifact;  the 
cluster  maintains  it3  shape  'under  repeated  heating  intervals.  Vhen  a  single 
atom  i3  removed  from  this  layer,  the  arrangement  becomes  'unstable  -  the 
cluster  breaks  up  into  iniiviiual  iridium  chains. 

Thi3  sort  of  behavior  i3  not  limited  to  iridium.  As  an  important  part  if 

our  studies  we  have  been  looking  at  silicon  layers  on  W(l!0).  Moderate  loses 

of  silicon  equilibrated  at  300  K  yield  a  very  open  overlayer,  shown  i.n  ?i*. 
7.3,  made  up  of  silicon  chains  oriented  along  r.''0|.  The  spacing  between 
chains  i3  close  to  9.5  A.  Closer  packed  layers  can  be  formed,  however,  by 
equilibrating  heavier  doses  at  higher  temperatures.  3y  field  evaporation  the 
concentration  of  silicon  adatoms  can  be  lowered  at  will.  Repeated  field 
evaporation  followed  by  equilibration  converts  these  ien3e  layers  into  chain 
arrays,  analogous  to  those  in  Fig.  7.  5.  On  further  reducing  the  alatom 
population,  these  overlayers  disorder  similarly  to  iridium  layers. 

The  formation  of  these  loosely  structured  layers  on  W(110)  has  not  been 
previously  observed  by  other  techniques  [Ro].  A  survey  has  been  made  of 
silicon  layers  on  other  planes.  3o  far  nothing  unusual  has  been  observed  on 
W(3P)»  (310)  or  (lOO).  It  is  already  clear,  however,  that  field  ion 
microscopic  examination  i3  capable  of  yielding  interesting  information  about 
overlayer3,  and  these  studies  a3  well  as  other  measurements  on  individual 
atoms  end  small  clusters  are  continuing. 
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3.1  Millimeter-Nave  Integrated  Circuits 

In  the  past  twelve-month  period,  we  have  made  significant  progress  in  the 
development  of  devices  and  structures  for  use  in  millimeter-wave  integrated 
circuits.  A  number  of  analytical  techniques,  applicable  to  the  dielectric 
waveguides  used  in  these  circuits,  have  al30  been  developed. 

\  particular  area  of  interest  has  been  that  of  dielectric  antennas.  We 
have  investigated  a  type  of  leaky-wave  antenna  which  consists  of  a  dielectric 
waveguide  on  which  are  placed  periodic  metallic  perturbations.  This  antenna 
i3  frequency  scannable  and  is  readily  integrated  with  a  receiver- transmitter 
system.  Prior  to  cur  investigation,  it  was  believed  that  the  perturbations 
had  to  be  'uniform.  Our  experiments  have  shown  that  ncnunifcrm  perturbations 
result  in  significantly  better  antennas.  Antennas  with  26  dB  suopresscn  cf 
siielcbe  and  endfire  radiation  levels  have  been  built.  These  antennas  can 
have  beamwidths  in  the  plane  parallel  to  the  antenna  axis  as  low  as  ^.0C. 


*  Thi3  work  was  supported  by  the  Joint  Services  Electronics  Program 
v’l.S.  Army,  U.3.  Navy,  'J.3.  Air  Force')  'under  contract  NOOO1  d-",R-G-'5 12  d. 
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For  some  applications,  the  radiation  pattern  of  this  antenna  in  the  plane 
transverse  to  the  guile  axis  may  be  too  broad.  To  reduce  this  beamwidth,  a 
new  structure  was  developed  in  which  the  dielectric  guide  with  periodic 
perturbations  on  top  is  embedded  in  a  metal  trough  with  a  metal  flare  attached 
on  each  side.  With  this  arrangement,  the  antenna  behaves  like  a  linear  array 
in  the  longitudinal  plane,  while  the  radiation  pattern  in  the  transverse  plane 
resembles  that  of  a  horn.  With  this  structure  we  can  achieve  7°  beamwidths  in 
the  transverse  plane. 

Many  passive  components  in  a  dielectric-based  millimeter- wave  integrated 
circuit,  such  as  ring  resonators  and  couplers,  are  comprised  of  curved 
waveguide  sections.  Therefore,  it  is  necessary  to  have  a  good  understanding 
of  the  behavior  of  the  field  in  the  neighborhood  of  the  bend  in  order  to 
predict  the  coupling  and  radiation  characteristics.  General  analytical 
approaches  to  this  problem  have  been  very  involved  to  data.  We  have  found  a 
simple,  approximate  technique  for  calculating  the  field  x.n  a  curved  dielectric 
waveguide  of  rectangular  cross-section. 

Dur  analysis  is  based  on  a  conformal  mapping  which  is  used  to  transform  a 
homogeneous  curved  waveguide  into  a  straight  waveguide  with  an  inhomogeneous 
refractive  index.  Solution  of  the  resulting  equations  yields  Airy  functions 
which  accurately  describe  the  field  both  inside  and  outside  the  waveguide. 

In  the  area  of  active  device  design  we  have  worked  on  the  problem  of 

IMPATT  oscillators  which  are  directly  integrable  with  dielectric  waveguides. 
Usually  this  results  in  excessive  mismatch  and  radiation  loss.  Through  the 
use  of  partial  shielding  structures  and  better  matching  networks  we  have 
reduced  these  losses  significantly  and  realized  oscillators  which  work  at  34 

GHz. 

A  new  method  of  solving  dielectric  waveguide  discontinuity  problems  i3 
al30  being  developed.  Unlike  a  closed  waveguide,  dielectric  waveguides 

support  radiation  modes  with  a  continuous  eigenvalue  spectrum.  This  makes  the 
application  of  mode-matching  and  related  techniques  very  difficult.  lur 

method  is  based  on  the  spectral  domain  method  and  utilizes  the  TFT  algorithm 
rather  than  matrix  inversion.  At  present  small  steps  on  dielectric  waveguides 
have  been  handled  successfully.  We  are  currently  inccrpo rating  a  variational 
expression  to  enhance  c cnvergence  for  larger  discontinuities. 
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3.2  Electromagnetic  Radiation  and  Scattering 

Thi3  area  of  our  research  is  primarily  concerned  with  the  efficient  and 
accurate  calculation  of  electromagnetic  radiation  and  scattering  problems.  We 
have  concentrated  our  efforts  in  the  high  frequency  region  where  moment  method 
solutions  become  impractical  and  where,  either  by  body  size  or  complexity, 
asymptotic  solutions  are  not  applicable.  In  the  past  a  spectral  domain 
approach  has  been  formulated  and  applied  to  thin  plates.  Recently  we  have 
been  able  to  accurately  extend  this  technique  to  two-dimensional  bodies  of 
significant  size,  which  may  be  either  perfectly  conducting  or  dielectric  in 
nature.  As  this  technique  is  iterative  and  uses  the  PFT  algorithm,  it 
contains  a  built-in  boundary  condition  checic  and  is  very  efficient. 

We  have  also  done  experimental  and  theoretical  studies  of  twist 
polarizers.  These  polarizers  can  be  used  in  radar  systems  where  traditional 
gimbaled  systems  are  too  slow,  as  the  twist  polarizer  i3  very  light  weight  and 
can  be  moved  rapidly.  The  theoretical  analysis  uses  the  spectral  domain 
technique  and  gives  good  agreement  with  experimental  results. 
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9 • t  Statistical  ? roperties  of  ? lasma  Turbulence 

The  object! re  of  this  research  is  to  3tuiy  the  statistical  properties  of 
plasma  turbulence.  The  system  under  investigation  is  the  positive  column  of  e 
1c- discharge  in  helium  which  is  'unstable  with  respect  to  the  spontaneous 
excitation  of  ion  acoustic  'saves .  Electric  field  fluctuations  cf  the 
turbulent  wave  field  are  picked  up  by  a  Langmuir  probe  and  their  probability 
distribution  is  measured  by  a  fa3t  sampling  method.  The  measured  deviations 
from  laussian  behavior  can  be  fitted  to  a  Irani- Iharlisr  expansion  of  the 
probability  distribution  including  terms  up  to  fourth  order.  The  non- 1  a U3siao 
character  of  the  probability  distribution  allows  conclusions  to  be  drawn  13  to 
the  nature  of  the  ensemble  of  the  waves  that  compose  the  turbulent  field. 
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10.1  Rarefied  Cas  Dynamics* 

Rarefied  gas  dynamics  deals  with  non- equilibrium  gas  flow  problems  in 
which  microscopic  treatment  according  to  kinetic  theory  i3  necessary  to 
determine  the  effect  of  intermolecular  collisions  and  gas  surface  interactions 
on  both  microscopic  and  macroscopic  gas  flow  properties.  Such  rarefied  gas 
flow  problems  occur  not  only  in  aerodynamics,  but  also  in  electronics, 
aeronomy,  environmental  fluid  dynamics,  and  other  related  fields. 

The  aim  of  this  research  program  is  to  develop  numerical  methods  to  solve 
a  -wide  range  of  problems  under  conditions  far  from  and  near  thermal 
equilibrium.  A  Monte  Carlo  method  has  been  developed  at  the  Coordinated 
Science  Laboratory  |_31  ]  to  solve  directly  the  3oltzmann  equation  and  has  been 
used  by  the  Boltzmann  group  to  solve  the  3oltzmann  equation  for  several 
rarefied  gas  flow  problems  under  a  wide  range  of  nonequilibrium  and  boundary 
conditions  [R2-S7].  The  solutions  we  have  obtained  yielded  detailed 
microscopic  and  macroscopic  non- equilibrium  properties,  most  of  which  have 
never  been  treated  and  studied  before.  We  have  also  studied  numerical 
solutions  of  other  kinetic  equations  and  other  numerical  methods  to  3olv 
rarefied  gas  flow  problems,  including  the  direct  simulation  technique. 

*  This  work  wa3  sunnorted  by  NATO  Research  Grant  1075  and  by  the  University  of 
Illinois. 
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Nonlinear  evaporation-condensation  problems  are  encountered  in  such 
iiversified  areas  as  upper  atmosphere  meteorology,  the  cooling  of  nuclear 
reactors,  design  of  space  experiments,  petrochemical  engineering,  vacuum 
technology,  and  the  interaction  of  high  power  laser  radiation  with  metal 
surfaces.  The  treatment  of  these  problems  requires  first  the  consideration  of 
vapor  kinetics  problems  characterised  by  the  nonequilbrium  vapor  motion  in  a 
Ynuisen  layer  at  the  interphase  boundary.  Under  a  .joint  research  effort,  we 
have  successfully  studied  the  evaporation  problem  [rs].  Our  Boltzmann 
solutions  establish  the  validity  of  the  Krook  solution  and  Ytrehus'  kinetic 
theory  approach  f R9 1 •  This  approach  can  be  used  to  calculate  the  jump 
conditions  across  the  Xnudsen  layer  and  the  net  mass,  momentum  and  heat  fluxes 
rR10l  and,  together  with  the  conventional  continuum  method,  to  calculate  the 
flow  parameters  at  the  outer  edge  of  the  Xnud3en  layer  for  problems  in  which 
evaporation  occurs  at  the  interphase  boundary. 

Studies  have  also  been  made  for  the  condensation  problem  r_10.2J.  The 
non- equilibrium  vapor  near  a  condensing  surface  differs  from  that  near  an 
evaporating  surface  because  of  the  difference  in  edge  condition  at  the 
interpase  boundary  and  its  behavior  is  more  complex.  Our  3oltzmann  solutions 
yield  distinct  non-equilibrium  behavior  in  the  Knudsen  layer  of  condensing 
vapor.  For  example,  the  temperature  gradient  becomes  negative  near  the 
interphase  -under  certain  conditions.  The  condensation  rate  as  a  function  of 
pressure  ratio  was  found  to  depend  strongly  on  a  substance  parameter  which  is 
proportional  to  the  ratio  of  the  latent  heat  and  its  liquid  surface 
temperature.  Some  of  our  results  were  found  to  be  in  agreement  with  those 
obtained  by  using  the  simple  apnroaches  of  Oguchi  ^ R 1 1  1  and  Ytreus  R 1  2 . 

Ne  plan  to  3tudy  further  the  evaporation  and  the  condensation  problems  by 
considering  more  complex  flow  and  boundary  conditions. 

10.2  Naval  Hydrodynamics* 

Free  surface  wave  problems  encountered  in  Naval  hydrodynamics  are 

characterized  by  complexities  in  flow  geometry,  flow  features  and  boundary 

conditions.  The  flow  has  an  unknown  free  surface  and  it  is  propagative  an  i 

*  This  work  wa3  supported  by  the  Office  of  Naval  Research  under  Oontract 
N0001 4-30-C-0740  and  the  National  Science  Foundation  under  Irani  NSF  ENl 
77-20436. 
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transient.  The  boundary  condition  at  the  free  surface  is  of  a  nixed, 
parabolic  type  and  it  contains  highly  nonlinear  terns.  In  the  steady  state, 

there  also  exists  a  radiation  boundary  condition,  since  the  waves,  once 

generated,  propagate  downstream.  These  complexities  have  led  to  several 
computational  difficulties:  accurately  accommodating  the  free  surface 

geometry,  satisfying  the  boundary  condition  uniformly  over  the  free  surface, 
treating  the  radiation  boundary  condition  and  implementing  the  boundary 
condition  at  the  contact  line  of  the  free  surface  with  the  3olid  surface  of  a 
partially  submerged  body. 

In  the  initial  phase  of  our  study,  we  have  developed  two  time-dependent 
numerical  schemes  with  Eulerian  grid  systems  for  solving  steady  and  unsteady 

potential  flows  for  nonlinear  free  surface  problems.  In  one  method 

[RIS.RU, R16],  we  use  the  finite  element  method  to  deal  -with  the  geometrical 
complexity  and  the  free  surface  boundary  condition  of  the  nonlinear  free 
surface  problems.  In  this  scheme,  the  finite  element  method  is  used  to  make 
the  field  calculation  and  the  finite  difference  method  is  used  for  the  time 
evolution.  We  have  used  thi3  method  to  solve  two  problems:  a  pressure 
distribution  moving  with  a  constant  speed  and  a  moving  submerged  elliptical 
cyliner  or  a  symmetrical  hydrofoil. 

We  have  also  developed  an  explicitly  time-dependent  finite  difference 
scheme  R 1 T 1 .  Explicit  schemes  may  be  more  favorable  for  solving  large  scale 
problems  on  "supercomputers"  which  have  more  stringent  implementation 
requirements.  We  have  used  this  method  to  3olve  two  problems:  a  pressure 
distribution  moving  with  a  constant  3peed  and  an  elliptic  cylindrical, 
surface-piercing  strut  accelerating  from  rest.  Our  solutions  serve  to 
demonstrate  the  feasibility  of  using  the  method  to  3olve  two-dimensional  as 
well  as  three-dimensional  problems. 

The  results  of  the  initial  phase  of  our  study  were  presented  at  the 
Second  International  Conference  on  Numerical  Ship  Hydrodynamics  R 1  5  j - 

In  both  numerical  schemes  mentioned  above,  the  computational  domain  is 
expanded  downstream  periodically  during  the  computation  as  the  disturbance  on 
the  free  3-urface  i3  propagated  close  to  the  downstream  boundary-  Therefore, 
the  ’undisturbed  condition  is  applied  on  the  cut-off  downstream  boundary. 
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In  order  to  obtain  the  longer  time  and  steady  state  solutions,  we  have  to 
use  a  fixed  downstream  boundary  set  olose  to  the  disturbance  so  that  we  can 
increase  the  computation  efficiency  by  dealing  with  a  small  computational 
domain.  The  present  phase  of  our  study  has  been  focused  on  the  application  of 
our  numerical  schemes  to  the  nonlinear  free  surface  problems  with  a  fixed, 
small  computational  domain.  There  are  three  major  tasks  involved  in  this 
3tudy:  (l  )  investigate  the  implementation  of  the  open  boundary  condition  at 
the  outflow  boundary  to  allow  the  waves  to  pass  through  it,  (2)  improve  our 
numerical  schemes  so  that  accurate  numerical  solutions  can  be  obtained  using 
the  open  boundary  condition  at  the  outflow  boundary,  and  (3)  apply  the  methods 
developed  to  selected  nonlinear  free  surface  problems. 

We  have  studied  in  detail  the  method  to  implement  the  open  boundary 
condition  used  by  Chan  [313].  It  is  based  on  Orlanski's  scheme  R 1 9  ]  in  which 
the  Summerfeld  radiation  condition  is  applied  at  the  outflow  boundary  and  the 
phase  velocity  is  calculated  numerically.  We  have  conducted  numerical 
experiments  to  study  systematically  the  errors  of  Chan's  numerical  scheme  so 
implement  the  advection  equation  and  to  study  way3  to  minimize  these  errors. 

We  have  studied  two  approaches  for  solving  potential  flows  for  nonlinear 
free  surface  wave  problems  using  a  fixed,  3mall  computational  domain.  In  each 
approach,  we  have  used  Orlanski's  method  to  implement  the  open  boundary 
condition  at  the  outflow  boundary.  The  two  approaches  differ  in  the  treatment 
of  high  frequency  errors  that  appear  in  the  solution  and  that  inhibit  accurate 
treatment  of  the  open  boundary  condition.  In  the  first  approach,  we  used 
filtering  to  minimize  high  frequency  errors,  however,  it  was  found  to  be 
effective  only  at  an  early  time;  therefore,  it  is  suitable  only  for  obtaining 
transient  solutions.  In  the  second  approach,  we  modified  our  time  marching 
schemes  so  that  the  schemes  are  characterized  by  damping.  This  damping  is 
more  effective  in  controlling  high  frequency  errors.  Transient  and  3teady 
state  solutions  have  been  obtained  for  the  pressure  distribution  and  the 
accelerating  strut  problems. 

Pig.  10.1  3 hows  the  pressure  distribution  problem  ’with  the  outflow 
boundary  set  close  to  the  disturbance.  The  appearance  of  a  two-grid  interval 
wave  in  the  solution  of  this  problem  near  the  outflow  boundary  is  shown  in 
Fig.  10.2.  The  eliminations  of  high  frequency  waves  in  the  solution  are 
shown  respectively  in  Figs.  10.3  and  10.4  for  two  methods,  using  filtering 


Fig.  10.3  Elimination  of  2Ax  waves  in  the  solution  of  the  pressure 

distribution  problem  (Fr  =  1//2tT)  at  t  =  1  by  using  filtering 


Fig.  10.4  Solution  of  the  pressure  distribution  problem  by  using  damping 


in  the  time-integration  scheme. 
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and  damping.  Fig.  10.5  3hows  the  steady  state  solution  obtained  by  using 
damping  for  a  fixed  computational  domain.  Fig.  10.5  shows  a  transient 
solution  of  wave  height  for  the  accelerating  strut  problem. 

The  results  of  the  second  phase  of  our  work  were  presented  at  the  Third 
International  Conference  on  Numerical  Ship  Hydrodynamics  [l0.4"j. 
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Research  in  thi3  area  is  concerned  with  computer  systems  which  have  both 
high  performance  and  high  reliability,  and  -which  can  be  developed  using  Very 
High  Speed  Integrated  Slrcuit  (7H3I3')  technology. 

1 ' . 1  Perfc rmance  of  Instruction  ? ipelines 

An  instruction  unit  pipeline  differs  from  a  functional  unit  oipeline  in 
that  it  i3  3ub  d  ec  t  to  delays  from  system  phenomena  external  tc  itself 
■/e.  g.  memory  and  functional  units'*.  Thus  a  segment  in  an  instruction  unit 

*  This  work  wa3  3uoocrted  by  the  Joint  3ervices  Electronics  Program  '7.3. 
Army,  7.3.  Navy,  7.3.  Air  Pores'*  under  contract  T3001  i-"3-3 -did  i,  by  the 
7.3.  Navy  under  contract  700037-37-3-J55C,  and  by  the  Hewlett-Packard  3cmpany. 
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pipeline  may  take  mere  than  one  time  unit  in  the  event  of  cache  miss,  memory 
collision,  lata  dependency,  etc. 

An  exact  solution  for  pipeline  performance  was  obtained  by  considering 
the  execution  time  of  any  given  segment  to  be  a  discrete  random  variable  which 
is  governed  by  a  multinomial  distribution.  The  expected  wait  time  of  an 
instruction  in  a  given  segment  can  be  recursively  defined  in  terms  of  the 
expected  wait  time  of  the  instruction  immediately  preceding  it.  If  the 
difference  of  these  wait  times  is  negative,  the  actual  wait  time  will  be  zero. 
Unfortunately,  this  solution  is  not  in  closed  form. 

The  effective  bandwidth  of  the  pipe  can  be  given  as  a  function  of  the 
expected  time  gap  between  the  ith  and  (i+1  Hh  instructions.  The  stationary 
value  of  this  expectation  requires  that  i  be  infinite,  but  it  was  found  that 
for  values  of  i  as  small  as  3.  estimators  can  be  obtained  for  the  parameter 
ranges  of  interest  -with  a  high  degree  of  confidence. 

A  closed  form  solution  was  obtained  heuristically  for  a  pipe  consisting 
of  N  +  2  segments.  The  two  end  segments  are  variable- time  segments,  and  the 
middle  N  segments  are  all  constant-time  segments.  The  purpose  of  this  study 
was  to  determine  the  degree  to  which  a  finite-length  buffer  could  decouple  the 
bad  behavior  of  two  variable-time  segments.  This  solution  was  found  by 
matching  boundary  conditions  (on  M)  against  3n  exponential  function  of  N,  and 
comparing  against  a  wide  range  of  simulated  values  (729  points  in  all).  This 
solution  had  an  error  range  <  1*.  with  an  RMS  error  of  <  2*  ever  the  simulated 
data  used  in  the  comparison. 

11.2  Test  Peneraticn  for  Pcmplex  Digital  Systems 

A  new  and  general  fault  model  for  microprocessors  has  been  developed,  and 
test  techniques  have  been  derived  tc  detect  all  the  faults  in  the  fault  model 
,11.3].  A  methodology  was  found  for  generating  combinational  structures  from 
high-level  descriptions  (using  assignment  statements,  "if”  statements,  and 
single-nested  loops’)  of  register-transfer  level  operators  r_  1  1  .1  .  The 
generated  structures  are  cellular,  and  are  interconnected  in  a  tree  structure. 
A  general  algorithm  has  been  developed  to  test  cellular  tree  structures  with  a 
test  length  which  grows  only  linearly  with  the  size  of  the  tree.  It  can  be 
proved  that  thi3  test  length  is  optimal  to  within  a  constant  factor.  Ways  of 
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making  the  structures  self-checking  hare  al3o  been  found. 


11.3  Tcncurrent  Trrcr  Ietection  in  Arithmetic  ani  Logic  'Jnits 

A  new  method  of  concurrent  error  ietection  in  Computer  Arithmetic  ani 
Logic  unit3  ha3  been  devised.  Thi3  methoi  exploits  the  inherent  physical 
redundancy  of  components  in  the  time  icmain.  Thi3  technique  i3  applicable  to 
all  operations  of  a  typical  ALL.  in  this  sense,  our  method  is  superior  to 
Arithmetic  Cedes,  which  are  applicable  to  only  the  arithmetic  cneraticns  of  an 
ALU.  furthermore ,  cur  methoi  also  works  with  many  types  of  carry  lock  aheai 
circuits.  In  terms  of  the  error  coverage  capability,  cur  mechod  compares 
favorably  with  a  iuplicatei  ALU  methoi  cf  error  ietection.  In  terms  of  the 
hariware  cost,  our  methoi  is  superior  to  both  the  Arithmetic  Ccies  ani  the 
Iuplicatei  ALU. 

■  1  .  1  High  ? erfc trance  and  Reliability  in  .Systolic  3 tructures 

3v3tolic  arrays  have  been  proposed  for  convenient  ani  efficient 
implementation  on  a  VLSI  chip.  However,  systolic  system  algorithms  orcocsei 
ic  net  use  all  the  inherent  parallelism  in  the  hardware .  He  have  found 
techniques  which  increase  by  a  factor  of  two  or  three  the  performance  cf 
3ystclic  systems  by  rearranging  the  sequence  cf  the  data  flow  ani  making 
3light  modifications  to  the  original  systolic  system.  In  addition,  a  high 
level  error  detection  technique  has  been  found  for  use  with  svstolic  systems 
performing  matrix  multiplication.  This  method  ices  net  use  much  extra 
hariware,  but  will  provide  protection  against  both  transient  ani  permanent 
erro  r3 . 

1 t . 5  Heme ry  Addressing  Architecture 

Tor  the  reduction  of  processor-memory  address  bandwidth,  we  have 
developed  an  address  oredicticn  mechanism.  An  address  predictionstack  holds 
the  addresses  cf  instructions  ani  data  which  have  a  high  orebability  cf  being 
referenced  soon.  The  reduction  cf  address  bandwidth  is  achieved  by  encoding 
positions  in  the  stack  structure  ani  sending  an  encoded  oc3iticn  tc  the  memory 
subsystem  as  an  address.  The  memory  system  decodes  thi3  3tack  ocsition  and 
aoplies  it  tc  an  identical  3tack  structure  maintained  -within  the  memory,  '-'sing 


the  information  in  it3  stack,  the  memory  system 


can 


K 


a  complete  address 
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A  full  address  is  transmitted  cn  a  3tack  miss. 

An  LRU  replacement  policy  and  several  line  cpdate  schemes  were  evaluated 
by  simulation.  The  simulation  results  provide  an  information  theoretic  lower 
bound  cn  the  average  number  of  address  bits  needed  per  memory  reference.  This 
bound  requires  frequency  encoding  of  the  stack  positions  and  mi3S  indicator. 
Conventional  fixed  length  binary  encoding,  which  is  mere  easily  implementable , 
requires  only  one  or  two  extra  bits  per  address.  For  the  frequency  encoding, 
as  the  stack  depth  increased  and  the  stack  width  remained  relatively  small, 
the  average  number  of  bits  needed  per  memory  reference  continued  to  decrease. 
This  ohenomenon  did  not  occur  in  the  case  of  the  conventional  binary  encoding 
because  after  a  certain  ooi.nt  a  price,  in  terms  of  extra  bits  needed  for 
encoding,  was  being  paid  for  the  greater  number  of  stack  locations.  Stacks  of 
modest  size  using  conventional  encoding  can  reduce  the  average  number  of 
address  bits  needed  oer  memory  reference  to  below  8  bit3,  even  when  the  total 
address  space  is  very  large. 

11.6  Shared  C ache  Organizations  for  Multiple-Stream  Ocmputer  Systems 

Organisations  cf  parallel-pipelined  multiple  instruction  stream 
processors  with  shared  two-level  memory  hierarchies  have  been  studied.  The 
cache  memory  interference  and  the  shared  cache  hit  ratios  in  such  systems  have 
been  investigated.  Performance  analysis  is  carried  out  by  using  discrete 
Markov  Ohain  and  probability  based  theorems.  Performance  is  evaluated  as  a 
function  cf  hit  ratio,  the  number  cf  processors,  the  number  of  pipeline 
segments  and  the  cache  organization  characterized  by  the  number  of  lines,  the 
number  cf  cache  modules  per  line  and  the  cache  cycle  time.  Come  design 
tradeoffs  are  discussed  and  examples  are  given  to  illustrate  a  variety  cf 
design  options  r 1 1 . 1 2  1 . 

11.7  Multiprocessors  with  ? rivate  Cache  Memo ries 

We  have  developed  an  approximate  analytical  model  for  the  performance  cf 
multiprocessors  -with  private  cache  memories.  All  private  cache  memories  are 
connected  to  a  single  shared  main  memory  via  an  interconnection  network.  We 
have  extended  an  earlier  analytical  model  to  include  complex  cache 
organizations  like  write- through,  load-through  and  buffered  .write-back. 
Considering  the  complexity  of  such  organizations,  cur  approximate  analysis  is 
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quite  good  ever  a  -file  range  cf  system  perimeters.  These 
cache  mi33  retie,  block  transfer  time,  eni  3ice 
interconnection  networks  l’.}1. 

'  1  .  3  Memo  ry  System  Archit ecture  fc  r  3  erellel  3  recessing 

Research  is  in  prepress  which  explores  the  design  cf  e  memcrv 
architecture  far  the  support  cf  parallel  processing.  Ve  have  defined  a 
preliminary  memory  system  architecture  which  ccn3i3ts  cf  a  set  cf  memory 
access  requests,  responses,  ani  the  semantic  action  associate!  with  each 
transaction.  This  architecture  will  serve  oo  define  a  communication  protocol 
which  can  he  applied  to  a  variety  cf  parallel  memory  structures  and  parallel 
interconnection  technologies.  The  parallel  architecture  allows  the  sharing  cf 
data  in  a  very  intimate  fashion  as  in  the  traditional  shared  memory 
multiprocessor,  while  also  supporting  the  importing  ''caching's  cf  shared  data 
for  temporary  use  to  exploit  data  reference  locality.  1  paper  describing  this 
shared  memory  architecture  has  been  submitted. 

11.9  Dcmputer  Architecture  for  high  1  eve 1  Languages 

Research  has  been  in  progress  in  the  area  cf  designing  directly 
interpretable  languages  'Til's')  for  high  level  language  support.  Techniques 
for  the  formal  derivation  cf  Dll's  designed  to  efficiently  support  high  level 
languages  are  described  in  a  paper  which  has  been  submitted  for  publication. 
Results  indicate  that  it  may  in  the  future  be  possible  to  systematically 
design  efficient  DIL's  from  a  formal  specification  cf  the  high  level  language 
to  be  supported.  Ve  are  currently  exploring  the  nossibility  cf  systematically 
designing  the  Dll  Interpreter  from  a  formal  specification  cf  the  syntax  and 
semantics  of  the  Dll.  The  objective  is  to  explore  techniques  for  a  complete 
systematic  design  of  a  high  level  language  3uopcrt  processor.  Tn us ,  the 
design  process  would  net  only  yield  efficient  candidate  Dll's,  but  also 
hard-ware  structures  for  the  efficient  interpretation  of  a  selected  Dll. 


oarameter3  include 
and  cyme  cf  the 


* i . 1 D  logic  Validation 


An  automatic  design  verification  facility  is  being  developed  for  cur  7131 
DAD  system.  Ve  have  made  progress  in  defining  the  orcblem  cf  design 


verification . 
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(l  )  determining  appropriate  design  language  restrictions  for  the  high-level 
functional  design  language,  (2)  translating  the  high-level  procedural 
description  of  the  system  function  to  an  appropriate  lev-level  representation 
which  is  suitable  for  comparison,  (5)  translating  the  hierarchical  description 
of  the  system  structure  to  the  low-level  representation,  and  (!)  performing 
the  comparison  between  the  low-level  descriptions  of  the  system  function  and 
the  system  structure.  Work  on  the  fourth  sub-topic  may  be  aided  by  recent 
worse  in  the  literature  on  an  analytic  verification  algorithm  'which  uses  binary 
decision  diagrams  to  represent  digital  functions.  We  hope  to  be  able  tc 

incorporate  such  a  verification  method  as  part  of  our  hierarchical  design 
verification  system,  since  this  approach  could  potentially  eliminate  the  need 
for  simulation  in  the  design  verification  process. 

11.11  Simulation  0 ackages  for  Complex  System  Design 

In  any  design  work,  simulation  becomes  an  important  part  of  the  design 
cycle.  Since  today's  logic  simulators  are  based  on  TTL  gate- level  structures, 
many  MOS  circuits  cannot  be  faithfully  modeled  using  them.  Circuit  level 

simulators  require  a  high  computational  and  memory  space  penalty  and  cannot  be 
used  for  large  circuits.  A  logic  level  simulator  called  7ETSIM  for  MOS 

integrated  circuits  has  been  developed.  The  simulator  takes  into  account 
slowdown  due  to  capacitive  loading  and  therefore  gives  accurate  tuning 

estimates  in  addition  to  logic  values. 

A  preprocessor  based  simulation  system  has  also  been  developed  that 
allows  the  specification  and  system  level  simulation  of  digital  systems.  The 
simulator  allows  hierarchical  specification  and  logical  connection  of  modules 
and  comprises  about  5000  lines  of  code  (including  direct  support  codeK  It 
has  been  used  to  examine  multiprocessor  bus  protocols,  memory  eganizatien  in 
multiprocessor  systems,  and  failure  modes  for  a  multiprocessor  system.  Since 
the  system  can  deal  with  a  multilevel  specification,  it  is  fairly  efficient 
and  a  real  time  to  simulation  time  ratio  of  50,000  for  systems  with  a  miner 
clock  cycle  of  10  ns  has  been  observed.  Theoretical  models  for  such  systems 
are  -under  development,  and  will  be  calibrated  using  the  simulator. 
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*2.1  Introduction 


Tie  primary  objective  of  the  research  reported  in  this  section  is  the 
development  of  efficient  computational  techniques  and  the  analysis  of  the 
capabilities  of  various  models  of  computation.  We  are  concerned  with  the 
resources  -  such  as  time,  equipment,  memory,  interconnection  -  either  used  or 
needed  in  the  algorithmic  solution  of  given  problems.  This  dynamic  discipline 
-  concrete  computational  complexity  -  is  not  only  contributing  to  our  basic 
understanding  of  computing,  but  it  is  extremely  relevant  to  actual  practice, 
both  in  hardware  and  in  software.  The  relevance  is  greater,  the  closer  the 
adopted  computation  models  are  tc  current  or  projected  computing  systems. 
Therefore,  it  i3  only  natural  that  the  great  technological  innovation 
represented  by  Tery-large-Scaie-Integratei  '7L3T)  Circuitry  ha 3  already  had  a 
substantial  impact  on  the  discipline,  opening  new  horizons  and  peeing 
challenging  problems.  Indeed  a  substantial  part  of  our  research,  organized 


below  in  four  subsections  '’’Parallel  Tomputation  in 
Geometr.",  "Approximation  Algorithms",  and  "storage 
Information  Retrieval"')  draws  its  motivation  from  this 


'/I.' I"  ,  "  Computational 

and  Access  Gosts  in 
Important  technological 


*  This  worm  was  suooortei  by  the  Joint  Services  Electronics  Program 
'7.3.  Army,  7.3.  Wavy.' 'I.  3.  Air  Force'  under  contract  *1001  i-~’ ?-C -0  42  4,  and  by 
the  National  Science  Foundation  under  grants  MG  3  “"3-1  5642,  MGS  3!  -05552,  MGS 
30-33354,  1ST  30-12240,  and  I3T  30-12242. 
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revolution. 

12.2  °arallel  Oomputation  in  VL3 1 

Parallel  computation  provides  a  fundamental  answer  to  the  need  of 
increasing  computer  processing  power.  Technological  advances  have  made  it 
possible  to  realize  parallel  computing  systems.  Our  studies  in  this  area  are 
concerned  with  the  development  of  realistic  computation  models,  the  conception 
of  powerful  interconnections  of  modules  (architectures)  realizable  in  the 
chosen  models,  and  the  analysis  and  design  of  algorithms  to  be  executed  on 
those  architectures.  Hi thin  the  commonly  accepted  "synchronous  model"  of  7L3I 
computation,  -ye  have  studied  the  realization  of  algorithms  and  have  developed 
area-time  optimal  schemes  for  matrix  operations  ( I  2. 1 3, 1 2. 1 9],  for  the  TFT, 
and  for  the  multiplication  of  integers  (l2.20].  He  have  also  developed  a 
critique  and  an  appraisal  of  current  VLSI  models  -  based  on  a  physical 
analysis  of  propagation  delays  -  and  have  begun  the  exploration  of 
three-dimensional  VLSI  techniques  (papers  in  preparation]. 

! 2. 3  Computational  Isometry 

Computational  geometry  is  relevant  to  a  number  of  applicational  areas, 
3uch  a3  operations  research,  statistical  analysis,  and  design  automation.  He 
have  continued  our  lines  of  investigation,  and  have  developed  time-optimal 
algorithms  for  a  variety  of  problems  concerning  two  families  of  orthogonal 
segments  (l  2.6,12.7];  this  problem  i3  very  relevant  to  VL3I  design.  He  have 
also  begun  the  exploration  of  a  class  of  problems  of  planar  realizability  in 
the  presence  of  obstacles,  and  have  developed  a  time-space  optimal  algorithm 
for  the  case  where  the  obstacles  are  modelled  as  parallel  segments  r_!2.16]. 


croximation  .Algorithms 


Motivated  by  the  intractability  of  many  important  problems,  we  are 
concerned  with  the  development  of  efficient  approximation  algorithms.  Cur 
wor'x  with  two-dimensional  bin  packing  (l  2. 1 , 12. 4, 12. 21  ]  models  certain 
scheduling  applications  and  also  relates  to  VL3I  placement  problems.  Ye  have 
recently  studied  7L3T  3treet  routing  problems  and  have  obtained  both  a  lower 
bound  on  the  number  of  tracks  required  and  have  developed  an  algorithm  which 
i3  at  most  431  worse  than  optimal  (papers  in  preparation!. 
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13.3  3 toraga  ani  A ccess  los £3  for  I nfo rmation  Retrieval 

la  any  info  rmation- retrieval  system,  efficient  storage  ani  accessing  of 
lata  are  essential.  Ye  are  studying  particular  lata  structures  ani  tine 
storage-access  traieoffs  inherent  in  performing  associate!  operations.  Tn i3 
is  leading  tee  the  development  of  analysis  techniques  which  can  he  applied  to 
various  retrieval  problems,  using  coat  measures  appropriate  to  the  particular 
applications . 

Ye  have  also  3tuiiei  r_1 2. 2, 1 2. 15  ]  a  general  class  of  resource  traieoffs 
that  arise  in  such  problems  as  parallel  sorting  algorithms,  linear  recursion 
schemata,  7L3I  layouts,  ani  searching  algorithms.  "These  traieoffs  span  a 
variety  of  resources,  allowing  us  to  traie  time  against  space,  area,  program 
3ic3,  ani  number  of  processors. 

A  lata  structure  with  local  access  permits  an  access  pointer  into  the 
structure  to  change  it3  value  from  one  location  only  to  an  adjacent  location 
in  the  structure.  Multiheai  Turing  machines  moiel  storage  ani  retrieval 
operations  in  these  lata  structures:  tree  machines  for  trees, 

multiiimensional  machines  for  arrays,  a h  have  u3«i  these  machines  to  ievelop 
new  methois  for  minimising  access  pointers  iato  trees  ani  arrays.  Tver y 
multiheai  tree  machine  of  time  complexity  n)  can  be  simulate!  on-line  by  a 

tree  machine  with  only  two  access  heals  in  time  3(t'n'  log  t(n)'log  log  n'l '  . 

Tver y  multiheai  e- iimensional  machine  can  be  simulate!  on-line  by  a 

i- dimensional  machine  with  only  two  access  heals  in  time 
3 ' r/  n)**(  1  *-  1  t  i  -  1  /ie'!  log  t'o>)  r_1  2. 1  1  2.  22 1 .  The  simulation  for  trees  i3 

optimal . 
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13.1  Natural  Language  'J riders tending* 

Research  has  continued  in  both  theoretical  and  practical  areas  relatei  to 
computer  understanding  of  natural  language  and  the  obtaining  of  information 
from  data  bases.  Attention  '.fas  focused  on  exploring  and  formalizing  the 
relationship  between  language  and  perception.  "Svent  simulation",  the  process 
of  constructing  a  sequence  of  sets  of  assertions  to  represent  the  meaning  of 
sentences,  was  proposed  and  investigated  as  a  method  for  .judging  the  relative 
plausibility  of  various  real  world  interpretations  of  sentences  describing 

*  This  work  was  3unported  by  the  Office  of  Naval  Research  under  jontract  No. 
N0QOU-75-C-061 1. 
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scenes  -and  actions;  event  simulation  was  also  3hown  to  be  useful  for  filling 
in  gaps  in  text  in  the  absence  of  "scripts".  Primitive  event  simulation  terms 
are  being  designed  to  provide  a  30und  basi3  for  models  of  language 
acquisition,  concept  learning,  language  generation,  and  "common  sense" 
_ 13.P9].  The  primitives  are  based  in  part  on  cur  study  of  Southeast  Asian  and 
Native  American  languages,  ani  the  methods  they  use  to  express  spatial  and 
event-related  information. 

The  problem  of  evaluating  natural  language  systems  wa3  investigated,  and 
a  number  of  novel  evaluation  techniques  and  dimensions  for  evaluation  were 
developed.  Evaluation  methods  were  designed  to  help  in  the  grinding 
improvement  of  the  tested  systems;  methods  were  tested  on  PLANES  ani  other 
natural  language  systems,  and  an  extensive  analysis  of  the  test  results  was 
carried  out  ani  documented  [13.3?]. 

•J9L,  the  relational  data  base  query  language  system  U3ed  in  PLANES  vas 
completely  rewritten  to  allow  "cooperative  responses",  probability,  and  the 
planning  of  optimal  data  base  search  strategies  :_!3.30l. 

Work  vas  completed  on  BROWSER,  a  system  which  automatically  attempts  to 
find  interesting  patterns  and  facts  in  a  data  base  '>3.311.  3R0WSER  is  driven 
by  a  model  of  what  is  interesting,  in  terms  of  statistical  patterns,  measures 
of  significance ,  and  models  of  the  real  world  actions  and  objects  which  are 
reflected  in  the  data  base. 

We  have  al30  been  working  on  a  computer  system  that  can  learn  from  it3 
experiences  with  the  world.  While  designed  primarily  for  a  natural  language 
processing  3y3tem,  the  research  i3  applicable  to  other  AT  domains  a3  well. 
The  technique  is  called  Explanatory  Schema  Acquisition  and  is  a  Wind  of  one 
trial  or  insight  learning;  with  this  technique,  a  natural  language  processing 
system  ‘will  be  able  to  learn  both  the  vocabularies  and  concepts  of  new 
knowledge  domains  directly  from  it3  input  text.  Tor  example,  a  natural 
language  3y3tem  that  knows  nothing  about,  say,  blackmail  will  be  able  to  learn 
that  concept  from  a  input  story  describing  a  blackmail  event  in  some  detail. 

Theoretical  work  has  been  carried  out  on  the  use  of  novel  "message 
passing"  algorithms  on  a  hexagonal  array  for  computing  high  level  visual 
features  'e.g.  .symmetry  md  shape")  of  objects,  and  for  defining  appropriate 
neighborhoods  for  objects  in  scenes  1  13.33"'.  Results  here  have  broad 
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significance ,  e.g.  for  VLSI  design  of  high  level  vision  processing  hardware. 

1  3 . 2  A_  Knowledge- based  System  for  Intelligent  Monitoring  and  Diagnosis* 

The  goal  of  this  research  is  to  develop  a  knowledge-based  expert  system 
capable  of  performing  intelligent  diagnosis  on  engineer-designed  systems.  As 
traditional  expert  systems  have  been  criticised  for  lack  of  understanding 
about  their  working  domains,  our  effort  is  to  .join  the  power  of  expert 
heuristics  with  the  completeness  of  general  domain  knowledge. 

In  our  approach,  an  expert  system  acts  as  the  front-end  for  a  diagnostic 
process.  It  efficiently  postulates  plausible  failures  to  explain  an  observed 
abnormality  based  on  the  powerful,  yet  not  perfect,  diagnostic  heuristics 
acquired  from  human  experts.  In  the  background,  postulated  fault  hypotheses 
are  scrutinized  by  a  rationalization  process.  Using  a  constraint  model,  the 
global  effects  of  each  hypothesis  are  thoroughly  deduced  and  then  compared  to 
observed  sensory  data  for  consistency. 

By  emphasizing  coordination  between  proposal  and  verification  processes, 
a  basic  expert  system  can  be  enhanced  by  the  deep- thinking  intelligence.  The 
philosophy  of  our  intelligent  expert  approach  is  useful  not  only  for 
diagnosis,  but  for  other  expert  problem-solving  systems  as  well. 

13.3  Automated  Air  Traffic  Contro 1** 

Air  traffic  control  can  be  subdivided  into  a  set  of  tasks  each  having 
only  one  primary  concern.  3uch  tasks  include  the  avoiiance  of  collisions, 
dangerous  weather  and  restricted  airspaces,  and  the  control  of  flow  rate  to 
airports.  However,  a  3et  of  techniques  to  perform  each  of  these  tasks 
independently  does  not  produce  a  system  which  performs  like  a  human 
controller.  What  i3  required  is  an  architecture  that  supports  a  constructive 
dialogue  among  the  experts  comprising  the  system.  Such  a  system  has  behavior 
that  depends  on  the  power  of  the  expert  interaction  language. 

*  Thi3  work  wa3  supported  by  the  National  Aeronautics  and  Soace  Administration 
under  contract  NASA  NAD  1 -20. 

**  Thi3  work  wa3  3uonorted  by  the  'J.3.  Department  of  Transportation  under 
Contract  XT  FA72-VA-T3oO. 
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Oar  research  has  produced  techniques  like  those  us  eh  by  the  experts,  a3 
'.fell  a3  a  planning  process  involving  the  use  of  criticism  and  suggestions  from 
a  boiy  of  experts.  Experts  constantly  monitor  traffic  for  problems  in 
parallel.  Should  an  expert  detect  a  problem,  it  plans  a  solution  and  submits 
it  to  the  other  experts  for  analysis.  This  process  may  cycle  before  a  final 
solution  is  adopted.  The  experts  share  a  common  conceptual  language 
describing  air  traffic  and  each  knows  its  own  relative  importance  in  the 
control  task.  Thus  the  experts  can  "argue”  and  decide  priorities.  Thi3 
architecture  emphasises  modularity  while  maintaining  the  cohesive  nature  of 
human  controller  performance. 

Still  the  system  ’will  not  improve  with  time  -unless  it  has  a  learning 
capability.  >ur  extensive  review  of  literature  on  learning  yielded  no  general 
theory,  although  some  simple  learning  has  been  demonstrated.  Vork  i3 
continuing  to  develop  learning  techniques  to  improve  expert  performance 
through  experience. 

1  S .  1  Theoretical  lomputer  7 ision* 

A  book  on  the  general  subject  of  modeling  spatial  patterns  was  completed 
^'"5.1  1.  Most  of  the  book  is  devoted  to  a  discussion  of  the  properties, 
33mthesi3,  extensions  and  perspective  display  of  mosaics,  coverage  patterns 
and  long  crested  wave  patterns,  which  are  explored  as  models  to  planar 
textures.  Model  based  texture  measures  are  also  discussed. 

Approaches  to  planar  decomposition  for  hierarchical  image  representation 
were  studied  Square  ani  triangular  quad  trees  were  found  to  be  the 

only  feasible  methods,  having  the  same  computational  complexity.  Tor  grii 
images,  the  type  of  grii  (triangular  or  square)  determines  the  appropriate 
quad  tree  type  (triangular  or  square). 

A  new  classification  of  low  level  image  models  was  proposed  under  che 
categories  of  pixel  based  and  region  ba3ed  models  I5.1!'.  Traditional  models 
were  surf  eyed  within  this  framework. 


*  Thi3  work  was  supported  by  NS?  under  grant  S0S-S1  lolOS. 
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1  3 . 5  ? Hot  Interaction  with  Automated  Airborne  Decision  Making  Systems* 

The  purpose  of  this  program  of  research  is  to  explore  the  issues 
surrounding  problems  associated  with  human- computer  interaction  in  flight 
management .  During  the  last  year,  three  projects  have  been  pursued: 
multi-task  decision  making,  measures  and  models  of  planning  behavior,  and 
onboard  computer-based  information  systems.  The  multi-task  decision  making 
project  has  been  involved  with  applying  the  modeling  and  function  allocation 
notions  developed  earlier  to  realistic  full-mission  simulations  13.131.  On 
the  basis  of  a  previous  review  paper  r32l,  the  planning  project  has  evolved  to 
a  series  of  experiments  to  study  the  planning  process  of  aircraft  pilots  in 
emergency  and  abnormal  situations  1)3.25,13.26].  The  computer-based 
information  system  project  started  with  an  initial  fairly  simple  experiment 
^13.1'7]  and  resulted  in  the  design  of  a  system  for  full-mission  evaluations 
J 3.27,13.23]. 

15.o  Human  Problem  Solving  Performance  in  Fault  Diagnosis  Tasks** 

The  purpose  of  this  project  is  to  develop  an  'understanding  of  human 
performance  in  fault  diagnosis  tasks  as  well  as  aid3  for  U3e  in  training 
humans  in  such  tasks.  In  the  last  year,  two  large  experimental  studies  of 
mantenance  trainees  were  performed.  The  first  experiment  was  the  third  in  a 
series  of  studies  of  human  performance  using  computer  simulations  of  aircraft 
power  plants.  The  goal  of  these  three  experiments  was  to  3tudy  how  suboptimal 
performance  is  influenced  by  a  variety  of  factors  r_  31 , 1 3. 35 1 .  The  second 
experiment  performed  this  year  involved  a  second  study  of  transfer  of  training 
from  computer  simulations  to  real  equipment. 

13.7  Modeling  of  Human  Behavior  in  Seeking  and  lenerating  Information*** 

This  project  is  studying  online  use  of  computers  by  scientific  personnel 
in  two  areas.  The  first  involves  the  use  of  online  editors  for  preparation  of 

*  This  work  wa3  sunnorted  by  the  National  Aeronautics  and  Brace  Administration 
under  Brant  No.  NS5-2119. 

**  Tnis  work  was  supported  by  the  U.  3.  Army  under  contract  MDA  9D3-''9-T-342!  . 

***  This  work  wa3  suonorted  by  the  National  Science  Foundation  under  grant 
I3T-791 6171  . 
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programs  and  documents  r_1  5.34-1.  The  second  involves  the  online  searching  of 
full- text  bibliographic  data  bases  rM.3ol.  The  goals  of  both  of  these 
subprojects  include  developing  mathematical  models  of  human  behavior  that  are 
suitable  for  use  as  design  aid3. 
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11.1  Introduction* 

During  the  1 930-1  931  time  period  the  Information  Retrieval  Research 
Laboratory  ( IRRL)  conducted  a  number  of  research  and  development  projects  and 
directed  the  operation  of  the  University  of  Illinois'  online  search  service. 
Major  activities  included  analysis  of  database  data;  a  test  of  the 
feasibility  of  automatically  determining  the  overlap  between  bibliographic 
databases;  development  of  a  computer-readable  database  directory;  and 
management  and  direction  of  the  University  of  Illinois'  online  search  service. 

11.2  Analysis  o f  Database  Data* 

The  IRRL  maintains  a  body  of  information  about  commercially  available 
databases  called  the  Database  of  Databases.  He  have  continued  to  study 
various  statistical  characteristics  of  this  population  based  on  the  material 
in  our  database,  ouch  statistics  ani  analyses  are  useful  both  to  researchers 
in  the  area  ani  to  the  U3er3  and  producers  of  databases. 

The  particular  information  developed  in  the  studies  performed  in  the  last 
year  considered  the  age,  size,  type  (scientific,  medical,  etc.'1,  and  source 
( government ,  private,  etc.)  of  the  databases  covered.  it  present,  ~T3 

*  This  worlt  was  supported  by  the  University  of  Illinois. 
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databases  are  covered  by  our  database  of  databases,  but  any  given  3tudy  tight 
include  only  a  specified  subset  of  those  depending  on  the  needs  of  the  3tudy. 
The  most  recent  analyses  are  for  the  year  !  930. 

The  software  used  in  this  research  consisted  of  a  program  framework 
containing  3lots  into  which  logic  could  be  inserted  to  perform  specific  data 
selection,  cleanup,  and  analysis.  The  iata  was  first  extracted  from  the 
database  by  our  data  access  program  USSR  and  written  into  a  file  of  (tag, 
value)  pairs.  The  analysis  programs  read  that  file  and  selected  the  tags  of 
interest.  In  effort  was  made  to  overlap  multiple  analyses  by  re-using  iata 
selected  and  cleaned  for  prior  use. 

Cross  tabulations  were  made  based  on  defined  groupings  or  actual  values, 
depending  on  the  iata  involved.  Thus,  cross  tabulations  were  done  involving 
both  raw  year  of  origin  versus  grouped  size  and  grouped  year  of  origin  versus 
grouped  3ize.  statistical  summaries  included  mean3.  ranges,  and  standard 
deviations  of  raw  iata  and  correlation  coefficients  between  raw  iata  and 
grouped  iata.  Detection  of  bad  iata  included  identification  of  records 
containing  illegal  or  missing  values. 

The  work  involved  in  thi3  project  also  developed  a  framework  for  further 
studies  of  the  databases  available  to  researchers  today.  'Ye  are  currently 
updating  the  database  and  improving  the  coverage  and  consistency  of  data 
fields  that  are  likely  to  be  important  in  future  studies. 

11.3  Integrated  Man/Machine  Interface  to  Facilitate  Network  Resource 

Utilization* 

For  the  first  time  in  history  computerized  information  retrieval  i3 
.widespread  and  economically  viable.  Evidence  of  this  is  shown  in  the  volume 
of  machine- readable  records,  databases,  online  software  packages,  vendors, 
search  services,  and  searches.  The  major  portion  of  the  currently  published 
scientific  and  technical  literature  can  be  identified  through  computer 
searches,  because  the  references  are  in  computer- read  able  form.  The  majority 
of  the  world's  currently  published  abstracting  and  indexing  literature  is  in 

*  This  work  .was  suooortei  by  the  National  Science  Foundation  under  Irani 
No.  T3F  031 
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computer- read  able  fora.  There  are  more  than  SCO  publicly  available  databases 
containing  approximately  more  than  ISO  million  recoris  (citations).  More  than 
"’5  million  of  these  recoris  are  now  available  through  the  various  online 
system  vendors  in  the  0.3. ,  Canada,  and  Europe.  Some  of  the  major  vendors 
include:  the  National  Library  of  Medicine  (MLM) ,  Lockheed  Information  Service 

(LI3),  System  Development  Corporation  (3DC),  Bibliographic  Retrieval  Service 
(3R3),  the  Canada  Institute  for  Scientific  Information  (CT3TI\  and  the 
European  Space  Agency  (ESA).  These  online  service  organizations  use  a  variety 
of  sophisticated  online  software  packages  such  as  ELHILL,  DIALOG!,  ORBIT, 
STAIRS,  CAM-OLE,  SRI?,  and  RECON.  Services  from  the  iatabases,  through  the 
online  vendors,  are  provided  by  dozens  of  centers,  libraries  and  brokers  and 
they  conduct  nearly  6  million  retrospective  searches  per  year. 

Unfortunately,  there  are  no  standards  governing  the  format  of  iatabases, 
the  elements  and  vocabularies  included  in  them,  or  the  online  systems  for 
searching  them.  The  lack  of  standards  forces  searchers  to  become  familiar 
with  each  existing  database's  coverage  and  vocabulary  and  with  each  system's 
access  protocol,  system  features,  and  command  language  and  to  keep  up  to  iate 
with  changes  made  in  all  of  them. 

Many  of  the  differences  and  variations  can  be  maie  less  apparent  by 
developing  translators  or  transducers  for  converting  the  procedures, 
conventions,  and  terminology  of  one  system  into  equivalent  procedures, 
conventions,  and  terminology  of  other  systems.  These  converters  would  make 
systems  appear  alike  to  the  searcher  and  make  databases  function  as  if  they 
used  the  same  vocabulary.  Such  aids  would  make  the  differences  transparent  to 
the  users;  thus  they  can  be  called  Transparency  Aids.  \n  integrated 
collection  of  transparency  aids  can  be  said  to  comprise  a  Transparent  System. 
The  aid3  then  are  the  Transparency  Subsystems. 

Mith  the  reduction  of  differences  the  iatabases  and  systems  become  easier 
to  use;  the  searching  environment  becomes  more  user-oriented.  In  a  more 
u3er-oriente d  environment  the  requirement  for  intermediary  searchers  decreases 
and  we  begin  to  approach  the  day  when  mo3t  searching  will  be  done  by  end  U3er3 
(those  that  po3e  the  search  questions) . 
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This  research  orogram  studied  net  ho i 3  to  increase  the  transparency  of 
retrieval  steps  by  making  the  discrecsness  of  3aparata  task3  lass  obvious  and 
by  x arcing  the  variety  of  system  and  database  differences  less  obvious.  Within 
this  project  we:  'i  )  designed  a  generalized  xodel  for  a  Transparent  System, 
(Rj  examined  alternative  factors  affecting  the  Transparency  Interface,  .3) 
built  coxponsncs  of  a  user-oriented  terminal,  and  (l)  integrated  example 
Transparency  Subsystems. 


The  Transparent  System  is  comprised  of:  user  interface  programs; 

languages  used  by  the  Transparent  System;  transparency  aids 
-ous  functions,  such  as  conversion  of  command  languages  or 
databases;  data  files  containing  descriptive  and  statistical 
ections  to  databases  and  network  resources, 
t  System  iay3  out  the  relationships  among  these 
iy  explores  advantages  and  disadvantages  of 
s.  Development  of  a  prototype  subsystem 
a  Transparent  System  by  showing  how  3ubsy3tem3 
can  interact  in  performing  their  functions. 
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The  generalized  model  includes  a  comprehensive  list 


existing , 


proposed,  and  potential  transparency  aii3  along  -with  estimates  of  their 
intercommunication  requirements.  The  alternatives  review  then  considers  the 
band  width  of  communication  paths  required  for  various  configurations  and  the 
effects  on  response  time  and  performance  of  changes  in  configuration .  The 
Prototype  for  the  Searcher's  Workbench  ' PTSW d  implements  a  structured  dialogue 
between  the  user  and  a  model  search  system.  The  dialogue  is  then  translated 
into  the  form  required  by  real-world  3y3tems  to  be  used.  The  structured 
dialogue  integrates  access  to  two  existing  transparency  aids,  a  database 
selection  adviser  and  a  vocabulary  transformer. 


1  !■ .  i  Directory  0  f  Tom  put  er-Reai  able  Bibliogra  phi  c  Databases* 

The  Tomputer-Readable  Bibliographic  Databases  -  Directory  and  Data 
Sourcebook,  compiled  and  edited  by  Professor  Rartha  E.  Williams  and  Sandra 
Rouse,  wa3  published  in  IT's  by  Knowledge  Industry  Publications ,  Inc.  'Ipdates 
to  the  Directory  were  issued  in  April  '  D"77  and  April  ’  ^7B.  A  new  ii-ect 

*  This  work  is  partially  supported  by  Knowledge  Industry  Publications,  Inc 
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including  data  that  was  relevant  as  of  December  ! 973,  was  published  in  autumn, 
! 37b.  The  computer  generation  of  the  printed  product  i3  supported  by  the 
American  Society  for  Information  Science.  'Jpdating  of  the  underlying  database 
in  preparation  for  a  new  edition  to  appear  late  in  1931  ia  now  underway. 

Data  Collection 

1 .  Database  processors  in  Europe  and  North  America  were  contacted  by 
mail  for  information  on  the  databases  they  provided  access  to. 

2.  The  17  page  questionnaire  designed  for  the  first  directory  va3 
revised  to  achieve  a  more  consistent  format  and  expanded  to  include 
non-'J.  3.  data.  A  turnaround  questionnaire  i3  being 
computer-generated ,  containing  the  information  in  our  files  for  the 
databases  covered  by  the  first  directory.  Existing  data  entries  need 
only  to  be  verified  or  changed  and  missing  information  supplied  by 
the  database  producer. 

3.  The  questionnaires  were  mailed  to  the  database  producers  and  :J.3.  and 
Canadian  producers  were  subsequently  contacted  by  telephone  for 
verification  of  the  data  and  assistance  with  questions  they  might 
have . 

d.  Information  on  the  availability  of  new  databases  was  obtained  from 
trie  published  literature  and  contacts  with  database 
processo rs/ producers .  The  producers  of  new  databases  were 

interviewed  first  by  phone.  The  completed  questionnaire  *as  then 
mailed  to  them  for  verification  and  if  necessary,  followed  up  by 
telephone.  Dnly  online  databases  that  were  publicly  available  were 
treated  fully.  Databases  available  for  internal  use  only  or  which 
are  only  accessible  offline  will  be  listed  separately. 

9.  Returned  questionnaires  were  checked  and  data  entries  standardized. 
Data  Entry 

1  .  Data  will  be  entered  online  into  the  database  of  databases  file  using 
the  03IN  program.  Entries  can  then  be  printed  out,  verified,  and 
corrected  and  edited  as  necessary  u3ing  the  DBEDIT  program. 


ITPORHATIH  337RI3VA.L  R35TA3UH 


1  i. 

Organization 

The  directory  is  organize!  alphabetically  by  database  name  and 
an  intro  lac tion  ani  four  indexes:  Sub  ject  Index, 

Name/ Acronym/ Synonym  Index,  Producer  Index  and  Processor  Index. 

Information  on  each  database  follows  the  general  format: 

(1  )  Basic  Information 

;2)  Producer /'Distributor''lenerator  Information 
(5)  Availability  and  Charges  for  Database  Tapes 
( i )  Subject  "latter  and  Scope  of  Data  on  Tape 
'5)  Subject  Analysis ''Indexing  Data 
(5j  Data  Clements  Present  on  Tape 
(7j  Database  Services  Offered 
(S')  "Jser  Aids  Available 

Pro  due  tion 

The  bibliography  will  be  computer  typeset  using  IRRL  software  ani  a 
soecial  printer  on  IRRL's  microcomputer  system.  Ail  indexes  will  then  be 
generated  automatically  from  the  iatabase  ani  typeset  similarly. 

1 i . B  University  of  I lino  is  Online  Search  Be rv ice* 

During  the  pa3t  year,  the  University  of  Illinois  Online  Search  Service, 
under  the  general  supervision  ani  coordination  of  IRRL,  expanded  its  services 
to  the  students,  faculty  and  staff  of  the  university.  In  cooperation  with  the 
University  Library,  searching  is  now  available  at  11  different  locations 
around  the  campus,  with  further  growth  pending.  Locations  where  searching  is 
offered  to  university  students,  faculty  and  staff  include  IRRL,  Reference  Room 
of  the  "lain  Library,  Chemistry  Library,  Biology  Library,  Labor  ani  Industrial 
Relations  Library,  Education  and  Social  Sciences  Library,  Agriculture  Library, 
Natural  History  Survey  Library,  Physics  Library,  Veterinary  Hedicine  Library, 
Health  Sciences  Library,  Agricultural  Economics  Reading  Room,  Sec logy  Library, 
and  the  Deologicai  Survey  Library. 

*  Thi3  wor!<  va3  supported  by  the  State  of  Illinois. 
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Online  systems  available  through  the  Search  Service  include  Lockheed' 3 
DIALD3,  System  Development  Corporation's  DR3IT,  the  Bibliographic  Retrieval 
Service  system,  and  the  National  Library  of  Medicine  MEDLARS  system. 
Individual  databases  available  from  these  systems  total  176. 

:Jse  of  all  systems  during  ?Y  31  is  given  in  the  following  table.  The 
costs  given  represent  the  total  amount  paid  by  the  University  to  the  online 
vendors  concerned.  Approximately  90*  of  this  was  recovered  from  the  users  of 
the  Search  Service  through  the  charging  of  fees. 

USE  DP  ALL  SYSTEMS  ?Y  31 


Connect 

Hours 

Costs 

3R3 

364.  0? 

23,566.  hi 

1R3 

1  26 . 36 

1!  ,423.91 

3  DC 

1 1  .70 

1 ,510.00 

NLM 

6.  SO 

156.00 

TOTAL 

1 007.16 

11  ,450.35 

The  total  connect  hours  U3ed  represent  approximately  2100  individual 
searches,  a  search  being  defined  as  a  single  session  at  the  terminal, 
regardless  of  the  number  of  databases  or  systems  accessed.  Df  these,  veil 
over  SDo  vere  conducted  and  paid  for,  at  least  in  part,  by  University  faculty, 
stuient3  or  staff.  Che  remainder  vere  conducted  by  IRRL  or  library  3taff 
either  as  practice  as  they  learned  to  a3e  the  system  or  for  general  reference 
or  'verification  vork. 

3eyond  providing  general  direction  and  coordination,  IRRL  performs  the 
following  specific  tasks  related  to  the  Dnline  Search  Service:  negotiation  of 
contracts  with  online  vendors;  training  of  new  searchers;  consultation  on 
system  problems  or  difficult  searches;  keeping  statistics  on  all  aspects  of 
the  Search  Service;  accounting  vork  associated  with  collecting  fees  and 
paying  bills. 
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1  1 . 6  Database  Analys.  e* 

IRRL  is  conducting  ongoing  analyses  of  the  MARC  database.  This  work 
involves  processing  the  annua'  files  distributed  by  the  Library  of  Congress 
and  generating  a  detailed  statistical  analysis  of  the  rate  of  occurrence  and 
sise  distribution  of  each  field  of  the  records.  The  resulting  tables,  which 
have  appeared  in  the  Journal  of  Library  Automation,  3 ho w  for  each  field  (i  ) 
the  percentage  of  records  in  which  it  appears,  [l)  the  rate  of  occurrence 
within  records  where  it  appears,  and  (O  the  length  of  the  field  in 
characters.  For  (’)  and  (5)  we  report  average,  minimum,  and  maximum  values 
and  standard  deviations. 

The  distribution  of  records  in  Dewey  and  LC  classes  is  al30  reported 
along  with  average  record  lengths  for  each  class.  Summary  statistics  are  also 
listed  separately.  Finally,  we  report  the  changes  between  reporting  periods 
(showing,  for  instance,  which  fields  are  used  more  now  than  at  the  last 
report)  . 

Me  have  also  begun  a  second  project  using  the  3ame  computer  programs  to 
analyse  the  D0J3ER  database  for  the  Douncil  on  Library  Resources. 

These  data  allow  users  of  the  databases  to  optimise  their  applications 
for  the  actual  characteristics  of  the  databases. 

11.7  Automatic  Detection  of  Database  Overlap** 

This  research  program  addresses  the  problem  of  overlap  coverage  of 
journal  articles  among  multiple  databases.  Overlap  leads  to  duplicative 
processing  of  material  by  database  producers;  duplicative  processing  and 
storage  of  material  by  online  and  batch  vendors  of  database  services;  and 
retrieval  of  duplicative  material  by  users  of  databases.  The  net  result  is 
reflected  in  cost3  associated  with  the  databases  and  an  unnecessary 
expenditure  of  time  on  the  part  of  users. 


*  This  work  wa3  supported  by  the  council  on  Library  Resources. 

**  Thi3  work  was  supported  by  the  National  Science  Foundation  under  Jrant 
Ho?  IDT  71 -V  31  3. 
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The  true  extent  of  the  problem  is  unknown.  Revere!  3tuiie3  heve 
eiires3ei  the  question  of  journal  coverage  overlaps,  some  heve  addressed  the 
question  of  article  overlaps  in  narrow  subject  areas,  and  one  has  attacked  the 
problem  of  identifying  maximum  possible  overlaps  among  specific  databases. 
None  has  studied  the  problem  of  determining  the  actual  overlaps  among  database 
services . 

'rfe  are  developing  a  methodology  for  automatically  determining  overlap 
among  databases  and  testing  and  implementing  this  methodology  on  seven 
databases:  AORI70LA,  3I0SI3  Previews,  OACon,  Txcerpta  Medica,  Medline, 

Science  Citation  Index  (301),  and  Social  Science  Citation  Index  (3371).  The 
methodology  involves  analysis  of  the  databases  to  determine  the  lata  elements 
present  in  and  common  to  multiple  databases  in  order  to  determine  effective 
keys  for  matching  records  in  multiple  files.  The  matching  process  will  be 
done  in  at  least  two  passes.  The  first  pass  or  broad  screen  will  use  a  match 
key  with  high  iiscrimination  power  to  identify  candidate  duplicate  records.  1 
set  of  finer  screens  will  be  used  for  verifying  pairs  of  records  as 
dunlicates.  The  second  pass  will  employ  Harrison  keys  (which  are  bit  string 
representations  of  titles)  and  Hamming  distance  techniques  (for  overcoming 
slight  variations  in  representation  of  titles)  together  with  matches  on  other 
elements  as  needed. 

Results  will  be  tested:  (l  )  by  using  manual  verification  of  a 
statistically  representative  sample  in  order  to  determine  mismatches,  and  (?'' 
by  using  a  known  set  of  duplicate  records  to  see  if  they  are  identified  by  the 
screens.  The  evaluation  in  terms  of  mismatches  (wrong  matches)  and  missed 
matches  will  provide  a  measure  of  the  effectiveness  of  the  techniques. 
Further  evaluation  in  terms  of  gross  costs  will  be  carried  out  in  order  to 
determine  the  efficiency  of  the  techniques. 

The  3tudy  will  produce  overlap  statistics  among  seven  databases,  data 
element  frequency  statistics,  and  a  methodology  for  duplicate  detection  that 
might  be  U3ed  by  database  producers  or  database  vendors.  The  resultant 
methodology  could  be  used  as  i  Dreprocessing  step  to  reduce  storage  and 
processing  cost3  or  as  a  post  retrieval  step  to  eliminate  the  duplicative 
outout  delivered  to  users. 
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1  t. 

Turin'?  thi3  first  year  of  the  project  tepee  have  been  acquire!  from  the 
iatabaae  producers  coverin?  the  i  R73  -  '  R7!  perioi.  cample  subsets  have  been 
real  ani  prelminary  planning  for  the  file  structures  ml  processing  flow  has 
been  lone.  The  contents  of  the  samples  have  been  stuiiei  statistically  to 
iientify  lata  elements  useful  in  the  matching  process  ani  to  allow  estimation 
of  storage  requirements. 
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15-1  Multiple-Terminal  Digital  Communications* 

This  ia  a  major  research  area  which  includes  various  problems  in 
spread-spectrum  communications  and  random-access  communications  networks.  The 
topics  that  were  investigated  during  the  past  year  include  the  performance  of 
spread-spectrum  multiple-access  communications  systems,  the  effects  of  fading 
on  spread-spectrum  communications,  communications  in  the  presence  of  jamming, 
conflict  resolution  in  random-access  communications  systems,  packet 
communications  systems,  and  signal  detection  in  multiple-access  channels. 
Progress  in  each  of  these  areas  is  described  in  the  subsections  that  follow. 


*  This  research  was  supported  by  the  Army  Research  Office  under  grant 
DAAG29-73-3-01  U  and  contracts  DAAG29-31  -'<-0062  and  DAAG29-31 -<-0064,  by  the 
•Joint  Services  Slectronics  Program  under  contract  M0001  d-7?-T-0421,  by  the 
Naval  Research  Laboratory  under  contract  NOO0H-30-C-0802,  and  by  the  National 
Science  Foundation  under  grant  3NG  73-06630. 
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'5. 

1  5 •  1  . 1  3 p real- 3  pec t  rum  Multiple-Access  lommunioation3  Systems 

luring  the  past  year  ve  have  obtained  many  new  results  on  the  problem  of 
evaluating  the  average  probability  of  error  for  asynchronous  direct -sequence 
spread  -  spectrum  multiple-access  (33*14)  communications  systems.  Because  of  the 
nature  of  the  multiple-access  interference,  standard  numerical  integration 
techniques  require  an  excessive  amount  of  computation  in  order  to  provide 
approximations  that  meet  a  prespecified  degree  of  accuracy.  We  have  developed 
upper  and  lover  bounds  vhich  can  be  evaluated  with  a  moderate  amount  of 
computation  for  direct-sequence  33*14  systems  vith  up  to  four  or  five 
simultaneous  transmitters.  The  bounds  are  presented  in  '.15]  along  vith 
numerical  results  vhich  illustrate  their  application. 

Ve  have  also  developed  accurate  approximations  to  the  average  probability 
of  error  for  direct-sequence  33*1.4  systems.  These  approximations  are  based  on 
the  integration  of  the  characteristic  function  of  the  multiple-access 
interference .  The  primary  advantage  of  these  approximations  is  that  the 
amount  of  computation  required  for  their  evaluation  grows  only  linearly  a3  the 
.number  of  simultaneous  transmitters  increases,  vhereas  the  amount  of 
computation  required  by  standard  numerical  integration  'and  most  other  known 
approximation  techniques)  increases  exponentially.  In  [13.3b”1  ve  present 
results  on  the  approximation  of  the  average  error  probability  for  binary  ?3X 
iirect-sequence  33*14  systems.  Extensions  to  quaternary  systems  vith  possibly 
offset  in-phase  and  quadrature  components  are  given  in  _3l].  In  ’.31  i  ve  also 
aliov  chip  waveforms  of  arbitrary  shapes  30  that  the  model  includes  as  special 
cases  the  forms  of  direct-sequence  modulation  known  a3  quadriphase-shift 
keying  (}?3K),  offset  QPSX,  and  minimum-shift  keying  (MSX).  The  approximation 
method  that  ve  present  gives  sufficient  accuracy  for  nearly  all  applications . 
However,  if  necessary,  even  greater  accuracy  can  be  achieved  by  using  this 
approximation  to  obtain  an  expansion  point  for  a  Taylor  series  representation 
of  the  actual  probability  of  error.  3y  employing  this  combination  of  the 
characteristic  function  method  followed  by  a  series  -  expansion  method,  ve  can 
obtain  any  prespecified  degree  of  accuracy. 


Eurther  results  have  also  been  obtained  on  the  performance 
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double- frequency  terms  ere  consiiere'i.  In  ^15.2]  the  signal- to-noise  ratio  is 
analyzed  for  general  quaternary  D3/SSMA  systems.  Numerical  results  are  given 
for  QPSK,  offset  Q?5K,  and  NSX  systems.  In  !j5.10]  optimal  phases  (with 
respect  to  signal- to-noise  ratio)  are  determined  for  m-sequence3  of  periods  51 
and  55.  "Jsing  the  phases  determined  in  [)5.10]  and  the  system  analysis  of 
M  5-2]  we  have  shown  that  the  signal- to-noise  ratio  can  be  significantly 
increased  by  careful  selection  of  the  phases  of  a  set  of  signature  sequences. 
Our  subsequent  worlc  on  the  probability  of  error  (e.g.,  ^1  5.591,  [Rl],  and 
[R5])  has  3hown  that  similar  improvements  in  the  bit  error  rate  are  also 
obtained.  The  effects  of  the  3hape  of  the  chip  waveform  on  the  performance  of 
direct-sequence  33NA  systems  has  been  investigated.  Results  on  the  sine  pulse 
vs.  the  rectangular  pulse  are  given  in  ^15.2]  and  [ri],  and  a  more  extensive 
3tud.y  covering  various  pulse  shapes  is  presented  in  lJ5«64]. 

Nearly  all  of  the  analytical  results  on  the  performance  of 
lirsc t-sequence  33MA  communications  systems  are  restricted  to  systems  for 
which  the  number  of  chips  per  data  bit  is  an  integer-multiple  of  the  period  p 
of  each  signature  sequence.  That  is,  there  are  an  integral  number  of  periods 
of  the  signature  sequence  in  each  data  bit.  7or  certain  applications, 
however,  the  necessity  to  discriminate  against  multipath  and  to  reduce  the 
probability  of  interception  dictates  the  use  of  very  long  signature  sequences. 
Typically  thi3  will  in  turn  necessitate  the  use  of  direct-sequence  modulation 
with  multiple  data  bit3  per  period  of  the  signature  sequence,  in  which  case  N 
may  be  much  smaller  than  p.  In  a  DS/SSNA  system  -with  multiple  data  bits  per 
period,  the  performance  depends  on  the  partial  correlation  properties  of  the 
signature  sequences  'rather  than  the  aperiodic  correlation  properties  as  in 
the  special  case  p  *  N.  In  l]l5.52]  we  analyze  the  performance  of  a  TT'TTdA 
system  for  which  p  >  N  and  N  and  p  are  relatively  prime.  This  analysis  taxes 
into  account  the  effects  of  partial  correlation.  In  particular  we  compare  the 
multiple-access  capability  of  two  classes  of  D5/33NA  systems  with  N  chips  per 
bit  when  the  performance  measure  i3  the  mean-3qu3re  value  of  the 
cro33Correlatio.a.  We  find  that  the  multiple-access  capability  i3  less  for  the 
cla33  in  which  p  and  N  are  relatively  prime  than  for  the  class  with  one  hit 
per  period.  Thi3  i3  primarily  because  of  improvements  that  result  from  the 
optimization  of  the  phases  of  the  signature  sequences  in  the  case  p  *  N. 
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In  addition  to  the  worse  described  above  on  direct-sequence  37d  A 
communications  systems,  several  investigations  of  freqency-hopped  33"?  4 
communications  systems  have  also  been  completed.  Both  fast -frequency-hopping 
',15.33]  and  slow- frequency-hopping  have  been  considered.  Bor  slow- frequency¬ 
hopping  we  have  investigated  DP3X  [l5.35]  and  BSX  J5.13]  data  modulation,  and 
results  have  been  obtained  for  hybrid  direc t- sequence  frequency-hopped 
spread- spectrum  systems  ',15.36]. 

15.1.2  Effects  of  Bading  on  5 p read-3 pact  rum  Communications 

Several  of  our  current  research  projects  are  concerned  with  the  effects 
of  fading  on  the  performance  of  spread-spectrum  communicat lon3  systems.  The 
fading  may  be  either  selective  or  nonsalective  and  either  Rician  or  Rayleigh. 
Both  frequency-hopping  and  direct- sequence  spread- spectrum  systems  have  been 
considered . 

In  U  5.4.0]  bounds  and  approximations  are  derived  for  the  average 
probability  of  error  in  an  asynchronous  BSK  slow-frequency-hopped  33’U 
communications  system.  Three  types  of  channels  are  considered:  the  nonfading 
additive  white  0iu33ian  noise  channel,  the  nonselective  Rician  fading  channel, 
and  the  doubly-selective  wiie-sense-3tationary  uncorrelated-scattering  ]‘V55TJ3' 
fading  channel.  The  two  types  of  hopping  patterns  consiierei  are 
deterministic  frequency-hopping  patterns  based  on  Reed-3olomon  codes  and 
rnniom  frequency-hopping  patterns.  Results  for  nonselective  Rayleigh  fading 
are  obtained  as  a  special  case  of  the  results  for  nonselective  Rician  fading. 

In  rJ5.35]  we  investigate  the  performance  of  3low-frequency-hopped  33 "1 4 
U3ing  3?3X  data  modulation.  The  channels  consiierei  are  Rayleigh  /T33U3  fading 
channels.  Brequency-selective  as  wall  as  time-selective  fading  channels  are 
considered  for  3y3tam3  with  various  data  pulse  waveforms  and  channel 
correlation  functions.  Analytical  results  for  the  average  probability  of 
error  are  presented,  and  numerical  examples  of  the  most  important  cases  are 
given. 

The  performance  of  direct-sequence  spread- spectrum  commun Lcations  in  a 
specular-multipath  fading  environment  is  investigated  in  ,!5.t1.  fhe  bounding 
methods  of  R5j,  which  were  developed  for  mul tipi  e-access  interference,  are 
applied  to  specular-multipath  interference  in  15.V'.  It  i3  shown  that  our 
approximation  to  the  probability  of  error  that  is  cased  on  the  signal- to  -  no ise 
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ratio  [|S.2]  also  gives  very  good  results  when  applied  to  the 

specular-multipath  problem. 

The  improvements  that  can  be  obtained  by  employing  error-correcting  codes 
on  a  slow-frequency-hopped  3SMA  system  have  been  investigated.  One  conclusion 
of  our  results  (e.g.,  [15.40])  is  that  in  a  Rayleigh  fading  channel,  bit  error 
rates  on  the  order  of  lO-’  to  10"">  cannot  be  attained  in  most 
3low- frequency-hopped  3SMA  communications  systems  without  the  use  of 
error-correcting  codes.  However,  the  application  of  certain  error-correcting 
codes  -/rill  bring  the  bit  error  rate  down  to  this  range. 

1 1 . 1 . 3  Communications  in  the  0 resence  of  J amming 

Last  year  we  began  an  investigation  of  the  problem  of  communicating  in 
the  presence  of  hostile  .jamming.  Ye  have  formulated  the  problem  in 
game-theoretic  terms,  using  Shannon's  mutual  information  as  the  payoff 
function.  Ve  showed  that  under  certain  restrictions  the  players  in  this  game 
(the  coder  and  the  jammer)  have  simultaneously  optimal  ( saddlepoint) 
strategies.  These  strategies  are  generally  memoryless,  and  so  in  an 
information- theoretic  sense  the  oommonly-u3ed  interleaving  strategies  are 
superfluous.  Ye  have  al30  determined  that  in  the  presence  of  binary  coding, 
uniform  Caussian  jamming  is  in  general  inferior  to  pulsed  jamming.  However, 
we  found  that  provided  a  molest  amount  of  band-width  expansion  is  tolerable, 
suitable  coding  can  neutralize  a  pulsei  jammer.  These  results  are  presented 
in  J5.  19;  for  systems  which  utilize  phase- shif t-’.ceying  modulation. 

For  the  case  of  frequency  shift  Ixeying  (FSX)  modulation  we  investigated 
the  capacity  in  the  presence  of  partial-band  jamming  and  nonselective  fading. 
Ye  found  that  there  i3  an  optimal  positive  code  rate  that  minimizes  the  energy 
per  bit  required  to  achieve  reliable  communication .  AI30,  we  showed  that  for 

low  code  rates  (C  .50  suffices)  ’uniform  jamming  is  optimal  for  the  jammer. 
Finally  we  analyzed  the  performance  of  a  specific  coding  scheme  for 
nartial-band  jamming.  These  results  are  presented  in  15.551. 

1 5. 1 • 1  Oonflict  Resolution  in  Random-iccess  Oommunication  Networks 

Conflict  resolution  consists  of  strategies  which  allow  remote  stations  to 
share  a  communication  channel  of  limited  capacity  in  a  decentralized  fashion. 
Which  strategy  i3  appropriate  depends  on  the  iemanis  of  the  stations  and  on 
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the  feedback  information  available  to  the  stations.  Ve  have  derived  upper 
bo'inis  on  random- access  throughput,  and  we  have  ieveloped  random- access 
strategies  for  situations  when  feedback  is  severely  limited  and  'under  several 
types  of  station  demands. 

1  bounding  technique  inspired  by  the  state-space  approach  of  stochastic 
control  is  introduced  in  [l5.3j.  It  is  used  to  establish  an  upper  bound  of 
.5125  packets/ slot  for  the  throughput  of  a  time-slotted  multi-access  broadcast 
channel  subject  to  an  infinite  population  of  user  stations  (modeled  by  a 
Poisson  process)  using  feedback  information  consisting  of  a  3,  1  ,  or  e  to 
denote  a  slot  with  aero,  one,  or  at  least  two  packets,  respectively. 

In  (15.15]  and  (15.4.2]  a  class  of  simple  recursive  retransmission  control 
strategies  is  given  for  the  random-access  broadcast  channel ,  again  subject  to 
an  infinite  population  of  users.  They  provide  stable  throughput  up  to  e*' 
packets/3lot  under  a  variety  of  feedback  models,  and  they  are  insensitive  to 
fluctuations  in  traffic  intensity.  Analytical  and  numerical  performance 
evaluations  are  given. 

A  clas3  of  random- access  strategies,  termed  acknowledgement  based 
retransmission  control  (A3RC)  strategies,  are  analysed  in  ]R3]  for  both  finite 
and  infinite  populations  of  users.  The  only  feedback  needed  to  implement  an 
A3R3  strategy  is  that  acknowledgement  of  a  successfully  transmitted  packet 
mu3t  be  given  to  the  transmitter.  3n  the  basis  of  an  equilibrium  analysis  and 
the  local  Poisson  approximation  introiued  in  ^’5.15]  it  is  shown  that  the  A3RT 
strategies  provide  satisfactory  stable  throughput  if  there  is  an  appropriately 
chosen  upper  bound  on  the  number  of  retransmission  attempts  allowed  to  any  one 
user.  Impressions  for  delay,  throughput,  and  probability  of  rejection  are 
found  and  suggestions  for  the  choice  of  parameters  are  given. 

Another  type  of  user  iemand  is  addressed  in  ]l5.44j.  Here  the  users  each 
have  a  steady  stream  of  packets  to  transmit  and  thus  a  TTH4  access  mode  is 
appropriate.  However,  if  there  is  limited  communication  among  the  users,  and 
moreover,  if  the  population  of  active  users  is  time-varying  (as,  for  example, 
in  a  mobile  communication  network)  then  3ome  strategy  for  locking  into  TT4A 
with  limited  information  must  be  used.  Preliminary  versions  of  such 
strategies  and  performance  evaluations  are  given  in  15. Id]. 
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15.1.5  Peneral  '1  umerical  Methods  and  Bounding  Techniques 
for  ? acket  Communication  Bys terns 

One  numerical  technique  we  have  investigate!  is  the  method  of  phase- type 
Markov  processes.  Phase- type  processes  are  Markov  processes  t  >  0 

where  the  component  (■?*.)  is  itself  a  Markov  process,  called  the  phase 
process.  For  example,  may  represent  the  number  of  packets  queued  for 

transmission  at  a  communication  station,  and  3,.  may  represent  the  phase  of  a 
TDM  A  protocol,  or  3.,.  may  model  fluctuations  in  incoming  traffic  intensity,  or 
both.  Numerical  result'3  of  this  application  are  reported  in  r56]. 

New  basic  results  in  the  computational  and  conceptual  aspects  of 
phase- type  Markov  processes  are  presented  in  ;15.1t].  In  particular,  a 
systematic  method  for  treating  general  boundary  states  is  found. 

The  method  of  phase-type  Markov  processes  is  also  applied  in  [m6]  to  find 
the  average  queueing  delay  in  a  link  of  a  communication  network  when  the 
incoming  traffic  is  obtained  from  a  Poisson  stream  by  a  generalized  round 
robin  (BRR)  routing  rule.  'Jnder  1R3.  routing,  successive  packets  coming  into  a 
node  are  routed  out  of  the  node  in  a  pre- determined  sequence  according  to  the 
order  of  arrival.  The  delay  for  3RH  strategies  is  often  signif icantly  smaller 
than  for  the  standard  independent  Bernoulli  splitting  routing  strategy. 

A  general  bounding  technique  based  on  drift  analysis  is  given  in 
Dimple  conditions  are  given  to  ensure  that  drift  implies  convergence.  One 
application  in  _R2j  is  a  proof  of  channel  stability  under  the  retransmission 
policies  introduced  in  1 1  5 - 1 5 ] •  A  stability  result  and  bounds  on  waiting 
times  in  DT/l/l  queues  are  also  presented  in  [ M2]. 

In  ■.  -14 ]  it  i3  shown  that  among  all  arrival  processes  (not  necessarily 
stationary  or  renewal  typed  for  an  exponential  server  queue  with  specified 
arrival  and  service  rates,  that  the  arrival  process  which  minimizes  the 
average  delay  and  related  quantities  is  the  process  with  constant  intararrlval 
times.  The  proof  13  based  on  a  newly  discovered  convexitv  property  of 
exponential  server  queues  which  i3  of  independent  interest.  Bince  the  traffic 
need  not  be  stationary  or  renewal  type,  the  theorem  provides  lower  bounds 
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(which  are  readily  computable  by  existing  methods')  to  the  average  delay  in  a 
network  of  queues  under  rather  general  routing  disciplines.  Some  related 
bounds  and  comparison  with  some  numerical  results  are  presented  in  3.6 1 . 

15.1.6  Signal  Detection  in  Multiple- Access  Noise 

A  project  in  multi-user  communications  that  has  been  initiated  recently 
is  the  study  of  optimum  signal  detection  procedures  for  application  in  the 
multi-user  environment.  Such  procedures  are  inherently  nonlinear  because  of 
the  non-0aus3ian  nature  of  the  multiple-access  interference  in  the  channel. 
Thus  the  initial  goal  of  this  study  is  to  determine  the  degree  of  performance 
improvement  exhibited  by  optimum  receivers  over  the  mors  conventional  linear 
(matched  filter)  receivers,  and  to  identify  those  situations  in  which 
performance  is  sufficiently  enhanced  to  justify  the  use  of  the  more  complex 
nonlinear  systems.  Early  progress  in  thi3  study  includes  the  derivation  of  an 
optimum  (minimum-error- probability)  receiver  structure  for  a  two-user  channel 
and  the  development  of  simplifying  approximations  to  thi3  optimum  receiver 
that  will  facilitate  the  desired  performance  analysis. 

15.2  Signal  Detection,  Estimation,  and  Filtering* 

Progress  in  thi3  general  area  ha3  involved  two  primary  continuing  areas 
of  study  as  well  as  the  initiation  of  a  new  area.  The  two  continuing  topics 
of  study  are:  (l  )  robustness  in  detection,  estimation,  and  filtering,  and 
(2)  the  design  and  analysis  of  digital  detection  systems.  Progress  in  these 
two  areas  is  described  in  the  following  subsections.  The  new  area  of  study  is 
signal  detection  in  multiple-access  channels,  and  early  progress  in  thi3  study 
is  reported  briefly  in  Section  15-1.6. 

15.2.1  Robustness  in  Detection,  Estimation,  and  Pilfering 

In  general,  the  area  of  robustness  deals  with  the  design  of  systems  and 
procedures  which  are  relatively  insensitive  (in  terms  of  performance)  to  small 
deviations  from  an  assumed  model.  Robust  techniques  are  of  interest  in  a  -tide 

*  T'ui3  research  was  3upoorted  by  the  Army  Research  Office  under  contract 
0AA329-31  -'<-0062,  the  Joint  Services  Electronics  Program  under  contract 
'10001 1-73-0—0124,  the  Office  of  Naval  Research  under  contract 
M0001 1-31  -0-001  d,  and  the  National  Science  Foundation  -under  grant 
ECS-79-l 6455. 
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variety  of  applications  including  radar,  3onar,  and  3ei3mology,  inasmuch  as 
inaccurate  models  are  the  norm  for  such  situations;  and  robustness  ha3  been 
studied  widely  in  the  contexts  of  communications ,  control,  and  inferential 
statistics.  The  basic  assumption  for  the  analytical  study  of  robustness  in 
detection  and  filtering  is  that  the  statistics  of  signals  and/or  noise  are  not 
'clown  exactly,  but  rather  are  known  to  be  within  some  (usually  nonparametric) 
classes  representing  ’uncertainty  in  the  'underlying  model.  The  primary 
technique  for  designing  a  robust  system  is  to  seek  a  system  achieving  the  best 
worst-case  performance  over  the  relevant  uncertainty  classes;  i.e.,  the 
primary  design  philosophy  is  minimax.  'rfithin  this  context  several  detection 
and  estimation  problems  have  been  considered  'under  this  project,  and  these  are 
described  in  the  following  paragraphs. 

The  problem  of  robust  linear  smoothing  of  a  stationary  random  signal  with 
uncertain  spectrum  observed  in  additive  noise  -with  uncertain  spectrum  is 
considered  in  f  15.21  1.  Here,  a  general  solution  to  this  problem  is  given  for 
spectral  ’uncertainty  classes  of  a  general  type  based  on  Tnoquet  capacities. 
Thi3  type  of  model  includes  standard  -uncertainty  models  such  as  contaminated 
mixtures  as  well  as  several  topological  models  of  ’uncertainty.  The  usefulness 
of  these  results  is  extended  in  (l5.5H]  by  generalizing  the  class  of  models 
that  can  be  treated  in  the  framework  of  fl5.2l].  Moreover,  an  extensive 
numerical  analysis  [tcoo]  indicates  that  the  worst-case  performance  of  the 
proposed  techniques  i3  generally  much  better  than  that  resulting  from  designs 
which  3imply  ignore  the  presence  of  uncertainty.  More  general  problems  of 
robust  estimation  of  stationary  signals  (including  smoothing,  filtering,  and 
prediction)  also  have  been  consiierei  for  the  case  of  discrete  time  J5-5-7]. 
Here,  a  general  miaimax  result  is  given  from  which  robust  solutions  to  a 
variety  of  problems  follow  straightforwardly. 

AI30  consiierei  in  this  general  context  is  the  problem  of  3tate 
estimation  in  linear  stochastic  systems  with  uncertainties.  Two  approaches 
have  been  considered  for  this  problem.  Dne  of  them  is  to  consiier  minimax 
designs  for  systems  with  -uncertain  state  and  observation  noise  statistics 
(|  5.29,15.51  1-  Here,  under  the  assumption  that  the  noise  is  white  with 
uncertain  componentwise  correlation,  it  i3  3hown  that  the  minimax - 
mean-square-error  state  estimator  is  the  Kalman  filter  for  a  least-favorable 
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model,  and  several  examples  are  considered.  Related  results  for  minimax 
linear-quadratic  control  are  found  in  'J5-47]  and  are  discussed  in  the 
Decision  and  Control  'unit  of  this  report.  Another  approach  to  state 
estimation  with  -uncertainties  concerns  systems  wi th  random  .jumps  which  are 
modeled  by  considering  the  standard  linear  system  model  wi th  adiitional  marked 
?oi3son  input  noise.  4  suboptimal  sequential  smoothing  algorithm  for  this 
model  is  developed  and  analysed  in  [15.51.  This  scheme  involves  detecting  the 
incidence  times  and  estimating  the  mark3  of  the  Poisson  input  process  and  then 
reconstructing  the  state  wi  th  a  (conditional)  Kalman  smoother.  The 
performance  of  this  scheme  is  seen  to  be  superior  to  that  of  optimum  causal 
estimators  and  linear  noncausal  estimators,  and  the  implementation  of  the 
scheme  is  superior  to  the  optimum  (nonlinear)  smoother,  which  i3  described  by 
an  infinite  sequence  of  filtering  equations.  .4  related  study  i3  found  in 
Ji.f]  which  considers  the  problems  of  minimax  (or  robust)  filtering  and 
state  estimation  for  observed  nonhomogeneous  ?oi3son  processes  with  uncertain 
rate  functions.  T.n  this  situation,  analogies  -.with  the  case  of  continuous 
observations  can  be  exploited  to  derive  minimax  designs. 

.several  topics  in  robust  signal  detection  have  also  been  considered. 

Three  basic  problems  have  been  treated  in  this  context.  Pirst,  the  role  of  a 
generalised  measure  of  signal- to-noise  ratio  (known  as  the  deflection)  in 
designing  robust  detection  systems  is  consiiered  in  [l5.23j.  It  i3 

demonstrated  that  this  criterion  is  often  much  simpler  than  earlier  criteria 
for  u3S  in  robust  design  because  a  constructive  method  for  finding  robust 
solutions  is  available  for  the  new  criterion  whereas  no  such  general  method  is 
available  for  earlier  formulations.  .4  second  problem  involves  the  design  of 
robust  systems  for  detecting  signals  in  the  presence  of  dependent  noise 

[15.26].  The  proposed  technique  uses  a  moving- average  model  for  noise 
dependence  and  considers  designs  that  are  minimax  robust  to  first  order  in  the 
degree  of  dependence.  The  third  effort  in  thi3  area  involves  finding 

solutions  to  the  problem  of  robust  discrete- time  matched  filtering  _'5.55]. 
The  problem  of  robust  matched  filtering  has  been  considered  previously  within 
a  general  Hilbert  space  formulation  J5.217],  and  J5.53;  derives  specific 
solutions  for  the  particular  problem  of  discrete- time  matchel  filtering  within 
a  variety  of  uncertainty  models  for  signals  and  noise. 


'  ... — - 
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15.2.2  Digital  Detection  Systems 

The  second  continuing  area  of  study  in  detection,  estimation,  and 
filtering  involves  the  analysis  and  optimum  design  of  digital  signal-detection 
systems.  This  problem  is  of  interest  in  a  vide  variety  of  applications  in 
which  receiver  digitization  is  increasingly  important  for  efficient 
implementation.  Recent  progress  in  this  general  area  can  be  divided  roughly 
into  two  main  topics:  optimum  data  quantization  for  detection  systems  and 
optimum  nonlinear  designs  for  signal  detection  in  dependent  noise. 

Progress  in  the  first  of  these  areas  is  reported  in  ^ 1 5.4],  _15.29j,  and 
1 1  5 . "50 1 .  In  1  5 . 4 1 ,  results  of  an  earlier  study  on  optimum  data  quantization 
for  coherent  detection  systems  are  extended  to  systems  for  the  detection  of 
purely  stochastic  3ignal3.  Dptimum  and  3uboptimum  schemes  are  compared,  and 
it  is  seen  that  significant  detection  efficiency  can  be  gained  by  implementing 
an  optimum  quantization  scheme  rather  than,  for  example,  by  using  standard 
minimum- ii3to rtion  quantization  in  a  quadratic  receiver.  These  results  are  of 
interest  in  sonar,  radio  astronomy,  and  seismological  applications  in  which 
signals  are  often  modeled  as  aondeterministic  stochastic  processes.  A  further 
study  of  quantization,  [l5.50],  considers  the  design  and  analysis  of 
quantizers  for  use  in  memoryless  (recursive)  signal  detection  systems  for 
detection  in  m-dependent  noise.  It  is  demonstrated  here  that,  by  utilizing 
the  dependence  structure  among  noise  samples,  one  can  improve  memoryless 
detection  performance  considerably  without  increasing  complexity.  Moreover  it 
i3  demonstrated  that,  in  the  presence  of  strong  dependence,  quantizer- 
detectors  utilizing  even  small  numbers  of  quantization  levels  can  be  nearly  as 
efficient  as  the  optimum  (Neyman-’earson)  detector  for  the  case  of  Daussian 
noise.  A  third  problem  considered  is  that  of  designing  optimum  input 
amplitude  compressors  for  coherent  detection  systems  utilizing  uniform 
quantization.  It  i3  shown  in  ^15- 29 ]  that  the  optimum  such  input  compressor 
i3  a  scaled  version  of  the  locally  optimum  detection  nonlinearity,  and  that 
the  performance  degradation  (over  locally  optimum  detection)  resulting  from 
the  proposed  quantization  scheme  is  negligible  for  only  a  moierate  number  of 
quantization  bits. 
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Procedures  for  designing  efficient  systems  for  detecting  signals  in 
dependent  noise  backgrounds  have  been  discussed  in  two  of  the  above  contexts, 
notably  M  5.26]  which  deals  with  robust  detection  in  weakly  dependent  noise 
and  [15. so]  "which  deals  with  optimum  quantization  for  detection  in  m- dependent 
noise.  Optimum  detection  for  the  weak-dependence  model  of  J5.26]  is 
considered  in  [lp.25].  It  is  shown  here  that,  for  a  moving- average  dependence 
model,  efficient  detection  can  be  achieved  by  reshaping  the  independent-noise 
optimum  detection  nonlinearity  with  an  additive  linear  correction  term.  in 
investigation  of  several  commonly  used  noise  models  indicates  that  the 
resulting  performance  improvement  i3  most  significant  for  impulsive  types  of 
noise.  In  [15.46]  similar  modified  structures  are  considered  for  detecting 
signals  in  4>-mixing  noise.  It  is  demonstrated  that  the  4>-mixing  model  with 
only  the  mixing  coefficients  given  is  not  sufficiently  descriptive  of 
dependence  structure  to  admit  a  design  which  is  uniformly  better  than  the 
independent-noi3e  design  (This  is  in  contrast  to  the  moving-average  model  of 
'  1  5.25]).  Ho  waver,  optimum  memoryless  designs  for  ^-mixing  and  other 
dependence  models  can  be  derived  using  general  results  of  1 5 » 22  ]  provided 
that  the  second-order  distributions  of  the  noise-process  are  known.  The 
optimum  system  in  this  case  is  the  solution  to  a  Fredholm  operator  equation  of 
the  second  kind  and,  as  such,  can  be  characterized  in  terms  of  the 

3arrett-uampard  expansion  of  the  noise  process.  A  further  result  in  this 
general  area  involves  the  comparison  of  memoryless  detection  systems  to 

systems  with  linear  memory  (such  as  the  Neyman-Pe arson  detector  for  the  case 
of  laussia.n  noise).  In  [15.34]  it  is  shown  that,  for  a  general  situation 
involving  stationary  noise,  memoryless  systems  are  asymptotically  as  efficient 
as  systems  "with  linear  memory.  This  indicates  that  implementation  complexity 
can  be  reduced  ‘.without  sacrificing  performance  for  3uch  situations  involving 
large  sample  3izes. 

15.5  Tata  Compression  Theory  and  Techniques 

lur  research  on  data  compression  has  been  extended  into  several  new 
cla33es  of  problems  including  the  construction  of  universal  source  codes, 
bounds  on  the  redundancy  of  universal  data  compression  techniques,  and  an 

investigation  of  the  redundancy  of  Huffman  codes.  Existence  results, 

redundancy  bounds,  and  construction  techniques  for  universal  source  codes  for 
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cla3383  of  memoryless  sources  are  given  in  [15.9].  The  codes  considered  are 
fixed-length- to-variable-length  (PL-VL)  source  codes.  New  results  have  been 
obtained  on  other  types  of  source  codes  including  variable- length- to- 
fixed-length  (VL-PL)  codes.  For  short  block  lengths  we  have  found  that 
properly  designed  VL-PL  universal  source  codes  provide  lower  redundancy  than 
our  best  PL-VL  source  codes  previously  reported  in  [l5*9]. 

Our  investigation  of  PL-VL  universal  source  codes  for  unifilar  Markov 
sources  has  provided  several  new  results.  A  computational  formula  was  found 
for  evaluating  upper  and  lower  bounds  on  the  redundancy  of  weighted  'universal 
coding  for  a  binary  two-state  Markov  source,  where  the  weighting  distribution 
of  the  transition  probabilities  is  ’uniform.  This  formula  was  also  analyzed 
for  large  block  lengths  and  found  to  agree  asymptotically  with  the  upper  bound 
on  the  redundancy  for  minimax  coding.  Parlier  results  on  upper  bounds  on  the 
redundancy  for  minimax  coding  for  'unifilar  sources  were  extended  to  cases 
where  the  structure  of  the  unifilar  source  was  unknown. 

Another  problem  -.which  has  been  investigated  is  that  of  the  redundancy  of 
Huffman  coding.  Huffman's  algorithm  is  a  well-known  procedure  for 
constructing  a  variable-length  code  of  minimum  possible  average  length  for  a 
source  -with  precisely  known  statistics.  The  redundancy  of  a  Huffman  code  is 
defined  to  be  the  difference  between  its  average  length  and  the  source 
entropy.  We  have  developed  a  powerful  general  technique  for  estimating  this 
redundancy.  It  implies  for  example  that  most  sources  have  redundancy  near  the 
value  logj log-j  *  .02976657 3,  and  th3t  apart  from  a  few  exceptional 

cases,  no  code  has  redundancy  exceeding  I  -  log,  e  *  log,  e  log,  e  * 
.736071302.  This  work  is  reported  in  ^15.43]. 
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16.1  Computer  Algorithms  for  Integrated  Circuit  Design* 

In  this  research  a  study  was  carried  out  on  the  use  of  partitioning  and 

tearing  algorithms  for  the  analysis  of  large-scale  circuits  3nd  systems.  In 

this  approach  the  circuit  is  partitioned  into  3ubcircuits,  each  subcircuit  is 

solved  separately,  and  then  the  subcircuit  solutions  3rs'  combined  to  obtain 

the  solution  of  the  entire  circuit.  Tearing  provides  many  computational 

advantages,  especially  in  the  simulation  of  large-scale  integrated  circuits. 

First,  computer  storage  can  be  reduced  if  identical  subcircuits  are  created. 

Secondly,  computer  time  can  be  saved  by  exploiting  subcircuit  latency;  that 

is,  subcircuits  that  are  inactive  during  a  given  time  interval  can  easily  be 

by-passed  in  the  solution  process.  Thirdly,  the  use  of  tearing  methods  allows 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
('J.3.  Army,  U.3.  Navy,’ ’  and  U.3.  Air  Force)  under  Contract  NOOOl  d-79-C-012i, 
and  by  the  National  Science  Foundation  under  Grant  SCS  30-0’70'75. 
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for  the  parallel  processing  of  subcircuit  equations.  Finally,  tearing 
facilitates  the  use  of  reduced-order  or  macromodel3  in  the  analysis  and 
permits  the  application  of  mixed-mode  simulation. 

In  our  research,  equation  solution  techniques  for  the  analysis  of 
circuits  and  systems  were  investigated.  In  particular,  when  tearing  methods 

are  employed,  there  are  a  number  of  solution  techniques  that  can  be  used  to 

exploit  latency.  These  techniques  have  been  studied  and  for  electronic 

circuits  the  most  efficient  solution  technique  has  been  identified  [l6.3S]. 
In  conjunction  with  this  work  a  fast  and  reliable  ordering  technique  for  the 
modified  nodal  equation  formulation  was  found  flS.S].  Our  ordering  procedure 
eliminates  the  time-consuming  pivot  selection  process  that  is  usually  required 
for  the  solution  of  the  modified  nodal  equations.  In  addition  to  the  above 
work  on  the  solution  methods,  various  latency  criteria  for  the  transient 
analysis  of  digital  circuits  'were  studied.  In  order  to  test  our  algorithms 

the  3PICE2  circuit  analysis  program  was  modified  to  include  tearing  methods. 
This  new  circuit  analysis  program  is  called  SLAT'S  [l6.3l]:  in  general,  thi3 
program  runs  faster  than  the  3PICS2  program  for  the  transient  analysis  of 
digital  logic  circuits  and  yields  the  same  accuracy. 

In  related  work,  algorithms  for  efficient  solution  updating  when  large 
changes  in  the  values  of  some  system  parameters  occur,  have  been  derived 
[16.4.].  This  problem  has  also  been  investigated  when  changes  occur  in  large, 
partitioned  systems  [l6.22]. 

In  addition,  the  problem  of  finding  the  shortest  path  in  a  network,  which 
occurs  in  network  routing  among  other  applications,  has  been  studied  and  an 
efficient  updating  algorithm  when  changes  occur  in  the  network  parameters  has 
been  derived  fl6.2ll. 

In  many  cases  electronic  device  models  are  represented  in  tabular  form. 
This  tabular  form  can  be  used  to  approximate  the  nonlinearities  by 
piecewise-linear  functions.  In  our  research  the  properties  of  the  solutions 
of  piecewise-linear  circuits  have  been  investigated  in  detail  ^1 6.5, 15.20 j. 

Research  has  also  been  done  on  the  use  of  piecewise  nonlinear  models. 
The  use  of  piecewise  nonlinear  models  for  the  exponential  nonlinearity  of  a  pn 
junction  offers  great  flexibility  in  modeling  1  6. 29, 1 5.  ^3].  However  these 
models  do  not  satisfy  the  requirements  of  the  standard  iterative  solution 
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techniques.  In  this  work  we  have  explored  new  solution  algorithms  tailored  to 
the  piecewise  nonlinear  models  and  compared,  computationally,  the  new 
algorithms  and  models  to  many  of  the  traditional  methods  used  in  programs  for 
the  calculation  of  the  DC  operating  point  of  bipolar  circuits.  Our  results 
indicate  that  the  piecewise  nonlinear  models  and  algorithms  compare  very  well 
to  the  best  methods  that  are  available  today  [l 6.30, 1 6. 36]. 

Future  work  in  this  area  will  include  theoretical  studies  of  the 
algorithms  for  piecewise  nonlinear  models.  Also  the  applicability  of  this 
type  of  modeling  for  use  in  MOS  circuits  will  be  studied. 

1  5.2  A_  New  Structured  Logic  A rray* 

The  use  of  highly  structured  architectures  is  becoming  necessary  to  take 
advantage  of  the  advent  of  VLSI  circuits.  One  area  that  lends  itself  to 
structuring  is  that  of  logic  arrays.  These  arrays,  such  as  ROM's,  PLA's,  gate 
arrays,  and  custom  logic  cells,  are  -widely  used  throughout  the  10  industry. 
The  success  of  these  structured  logic  arrays  is  in  large  part  due  to  the 
development  of  successful  GAD  tools  for  their  design,  layout,  and  routing. 

There  are  two  broad  classes  of  structured  arrays:  customer 

reprogrammable  and  customer  nonreprogrammable .  In  this  work  [l 5.341  we  are 
investigating  a  new  architecture  for  a  customer  reprogrammable  logic  array. 
This  array  structure,  known  as  REPLICA  (Random,  Electronically  Programmable, 
Logic  Integrated  Custom  Array ^  is  a  regular  structure  of  NOR  gates  with 
restricted  fan-in  and  fan-out.  We  have  developed  CAD  tools  for  the 
minimization  of  logic  realizations  in  REPLICA  similar  to  the  folding 
algorithms  used  for  PLA's.  We  have  also  developed  guidelines  for  the  optimal 
form  for  the  logic  equations  to  be  Implemented  in  REPLICA.  Several  examples 
of  logic  functions  implemented  in  REPLICA  indicate  that  it  competes  well  with 
PLA's  and  has  some  features  that  may  make  it  more  iesireable  than  PLA'3  in 
situations  such  as  realizations  of  sequential  logic.  Work  is  progressing  on 
comparisons  of  REPLICA  with  existing  structured  arrays  and  on  improving  the 
CAD  tools  necessary  for  realizing  larger  test  examples  in  the  REPLICA 
architecture . 

*  Thi3  work  was  supported  by  the  Joint  Services  Electronics  Program 
/-J.3.  Army,  ’J.5.  Navy,  ’  and  U.S.  Air  Fores'!  under  Contract  N0001  A-h-C-OiCl. 
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1  o  •  3  Hierarchical  Design  Techniques  for  VLSI  Systems* 

With  the  complexity  of  the  systems  that  can  be  realized  on  a  single  VLSI 
chip  the  design  time  is,  using  conventional  techniques,  on  the  order  of  100 
man-years.  In  order  to  reduce  this  unacceptable  design  time,  structured 
design  techniques  are  being  investigated.  In  this  work  [16.37]  the  impact  of 
a  structured  design  methodology  on  the  design  of  a  VLSI  system  was  studied. 
The  design  and  testing  of  the  chip  were  both  considered  from  the  beginning  of 
the  design  process.  As  a  test  example  a  CORDIC  function  (used  for  realizing 
trigonometric  functions)  was  designed.  The  function  to  be  implemented  was 
studied  by  simulation  in  order  to  finalize  the  architecture  of  the  design. 
Then  followed  a  structural  design,  a  physical  design  including  layout,  and  a 
verification  using  different  simulation  tools. 

The  investigation  clearly  showed  that,  with  some  sacrifice  in  chip  area, 
a  structured  design  approach  was  feasible  for  VL3I.  Various  other  systems  are 
now  being  studied  for  implementation  in  VLSI  using  structured  design 
techniques.  The  goal  of  these  studies  is  the  refinement  of  the  tools  and 
techniques  necessary  in  the  design  hierarcy  for  VLSI  systems. 

1 6. 4  The  Use  of  Multiple  Criterion  Optimization  for  the  Design  of 

Digital  Filters** 

The  design  of  digital  filters,  both  FIR  [l  6. 12, 16.32]  and  HR,  using 
competing  criteria  f  1 6.3]  is  being  studied.  In  particular,  various 
formulations  of  the  design  problem  which  trade-off  between  magnitude  and  phase 
performance  are  being  considered.  These  formulations  involve  the 
specification  of  a  desired  phase  behavior  as  well  as  the  standard  desired 
magnitude  behavior. 

Performance  constraints,  such  as  monotonic  magnitude  response  in  the 
transition  regions,  are  also  Imposed  on  the  behavior  of  the  filter.  In  order 
to  30lve  these  rather  difficult  optimization  problems  we  are  modifying  the 
sequential  quadratic  method  for  constrained  optimization  proposed  by  Han  and 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
(7.3.  Army,  7.3.  Navy,  and  U.3.  Air  Force)  under  Contract  M0001 4-'r9-C-OU24-. 

**  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
(7.3.  Army,  7.3.  Navy,  and  U.3.  Air  Force)  under  Contract  N0001 1-79— C-D421, 
and  by  the  National  ocience  Foundation  under  Grant  ECS  30-070"'5. 
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Powell.  We  are  also  exploring  the  use  of  a  dynamic  frequency  grid  which  will 
allow,  theoretically,  convergence  to  the  true  solution  of  the  design  problem 
and  not  the  approximation  that  is  found  in  all  existing  digital  filter  design 
programs  today. 

16.5  Analytical  and  Computer-Aided  Design  of  New  Circuit  Structures  for 

Switched  Capacitor  Filters* 

Research  is  being  conducted  in  optimum  switched  capacitor  (3C)  filter 
structures  and  in  computer-aided  analysis  and  design  (CAD)  of  SC  circuits. 
Recent  studies  of  optimum  SC  filter  structures  resulted  in  two  different 
designs  using  the  bilinear  z-transformation.  One  design  uses  a  state  space 
technique  [l  6. 1 1 , to. 14, 16. 17]  while  the  other  uses  a  follow- the-leader 
feedback  (FL?)  topology  [l6.l0].  Both  techniques  allow  SC  filters  to  be 
implemented  with  stray  insensitive  first  order  modules  (such  as  forward  Euler, 
backward  Euler,  or  bilinear  descrete  integrators)  or  second  order  biquadratic 
modules.  Optimization  procedures  were  applied  to  maximize  dynamic  range  while 
minimizing  sensitivity,  capacitance  ratios,  and  total  capacitance.  Impressive 
improvements  in  sensitivity,  as  analyzed  by  Monte  Carlo  analyses,  have  been 
obtained  for  an  example  of  an  optimized  fifth-order  Chebychev  LP  filter 
designed  with  the  state  space  and  FLF  approaches. 

In  the  area  of  CAD,  a  user-oriented  program,  SCAP  II,  based  on  modified 
nodal  analysis,  was  completed.  Frequency  and  sensitivity  analyses  were 
demonstrated  on  a  wide  variety  of  useful  SC  circuits  taken  from  the  literature 
[i6.25,16.26,16.55].  SCAP  II  is  capable  of  performing  frequency  analysis  and 
frequency  domain  sensitivity  analysis  on  SC  circuits  that  contain  M03 
switches,  capacitors,  and  voltage  controlled  voltage  sources,  and  which 
operate  -with  a  2-phase  50^  duty  cycle  switching  sequence.  The  analysis  is 
efficiently  implemented  with  sparse  matrix  numerical  methods  and  a  specially 
designed  equation  ordering  strategy  which  minimizes  fill-in  during  computation 
of  the  LU  factors.  Sensitivity  analysis,  which  is  an  implementation  of  first 
order  relations  obtained  by  differentiating  the  original  modified  nodal 
equations,  is  efficiently  implemented  for  capacitors,  capacitor  ratios,  and  op 

*  This  work  is  supported  by  the  National  Science  Foundation  unier  Grant 
3NG-7B-1  V31  5. 
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amp  gains  by  using  the  same  LU  factors  as  already  computed  for  the  frequency 
analysis.  3CAP  II  requires  approximately  twice  the  machine  time  as  the 
previously  reported  block-partitioning  technique  of  5CAP  I,  and  hence  appears 
to  be  a  useful  aid  in  SC  circuit  design  [ 1 6. 7, 1 5. 23] . 

The  most  recent  activity  in  SC  CAD  has  been  the  development  of  a  very 
general  N-phase  SC  analysis  routine,  referred  to  as  SCAPN  [ 1 6 . 13]-  In 
addition  to  using  sparse  matrix  methods  and  proper  reordering  techniques  SCAPN 
uses  a  specially  designed  "pre-L'J  factorization  algorithm"  which  minimizes  the 
amount  of  computation  that  must  be  done  to  update  the  LU  factors  at  each  new 
frequency.  SCAPN  also  features  a  very  efficient  sensitivity  analysis  based  on 
the  adjoint  network  concept  [lo.l-t],  and  a  nonlinear  distortion  analysis  based 
on  the  Volterra  series  technique  for  analyzing  the  effects  of  driving  MOS 
capacitors  and  amplifiers  into  nonlinear  operation  [ 15.1 5 1  - 

15.5  Automatic  Tuning  Algorithms  for  Analog  Filters* 

To  meet  filter  response  specifications  analog  filters  usually  must  be 
tuned  or  adjusted,  preferably  by  computer  automation  if  the  production  level 
is  high.  In  this  research  three  generalized  tuning  algorithms  have  been 
studied  on  the  basis  of  their  architecture,  computational  complexity,  and 
effectiveness.  Further ,  a  method  has  been  found  for  the  tuning  element  and 
frequency  selection  problem,  a  problem  relevant  to  all  three  methods.  Monte 
Carlo  simulations  have  been  done  on  several  active  filter  circuits  in  order  to 
enhance  the  comparison  and  provide  a  demonstration  [l5.23|« 

15.7  Fault  Isolation  in  Analog  Circuits** 

The  objective  of  this  research  is  to  be  able  to  isolate  faulty  components 
in  analog  circuits  from  test  point  measurements.  Presently  the  only  approach 
that  has  had  limited  success  in  practice  is  the  fault  dictionary  methoi. 

*  This  work  was  supported  by  the  National  Science  Foundation  under  Grant  3NG 
79-117*58. 

**This  work  was  supported  by  the  U.3.  Naval  Air  Fngineering  Center  under 
Contract  M63T55-30-C-0505  and  by  the  Joint  Services  Slectronics  Program 
(U.3.  Army,  U.S.  Navy,  and  U.3.  Air  Force''  under  Contract  N0001  i-'?— C-0421. 
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Unfortunately  this  method  requires  that  one  anticipate  all  possible  fault 
conditions  in  advance.  This  is  usually  an  impossible  task,  and  reliable  field 
lata  is  required  in  order  to  update  the  dictionary.  Typically  it  is  also 
necessary  to  confine  the  dictionary  to  single  catastrophic  faults  in  order  to 
avoid  excessive  storage  requirements  and  long  search  times  for  the  fault 
conditions.  Finally,  in  the  fault  dictionary  it  is  not  easy  to  determine  the 
testability  of  the  circuit  from  a  given  set  of  test  points. 

In  our  work  we  have  shown  that  a  sensitivity  matrix  can  be  used  to 
isolate  both  hard  and  soft  faults.  The  idea  behind  the  method  is  that  if  the 
columns  of  the  sensitivity  matrix  are  linearly  independent,  then  any  fault  in 
a  component  will  cause  a  unique  shift  in  the  test  point  voltages.  The  method 
can  locate  single  as  well  as  multiple  faults  in  analog  circuits.  Furthermore, 
the  independence  of  the  columns  of  the  sensitivity  matrix  are  a  good  measure 
of  the  testability  of  the  circuit  -with  respect  to  the  given  test  points.  The 
sensitivity  matrix  is  computed  for  the  linearized  circuit  model  at  the  nominal 
parameter  values  in  the  pre-test  phase.  Thus,  only  a  few  simple  mathematical 
operations  are  required  in  the  test  equipment.  A  sensitivity  matrix  must  be 
computed  for  each  test  signal,  and  assuming  that  the  rank  of  the  sensitivity 
matrix  is  equal  to  the  number  of  test  points  for  any  column  ordering,  then  the 
number  of  faults  in  the  circuit  must  be  less  than  the  number  of  test  points  in 
order  for  the  method  to  give  reliable  results. 
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17.1  Control  and  Decision  Strategies  for  Systems  Under  Imperfect  Information* 

Uncertainties  arise  because  of  unknown  system  parameters,  unknown  signal 
environments ,  and  hardware  tolerances.  Related  complexities  ari3e  in 
situations  involving  multiple  system  performance  criteria  and  multiple 
decision  makers.  The  objective  of  this  project  is  to  gain  a  basic 
understanding  of  the  behavior  and  control  of  complex  3ystem3.  During  the  past 
year,  attention  has  been  given  to  several  topics.  Highlights  are  summarised 
in  the  following. 

*  Thi3  work  wa3  supported  by  the  Joint  Services  Electronics  Program 
('J.3.  Army,  U.3.  Navy,  and  U.5.  Air  Force)  under  Contract  N0001  l-7?-C -042d . 
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Several  results  concerning  the  sensitivity  and  robustness  of  control 
systems  to  'uncertainties  in  the  plant  model  have  been  obtained.  We  have 
developed  bounds  on  the  magnitude  of  allowed  perturbations  which  guarantee 
that  the  perturbed  system  performs  better  than  an  equivalent  open  loop  system. 
These  bounds  are  expressed  in  terms  of  the  singular  values  of  the  return 
difference  matrix.  We  have  also  developed  the  singular  value  sensitivity 
function  as  a  robustness  analysis  tool.  Singular  value  sensitivities  can  be 
used  to  analyze  perturbations  'such  as  multiplying  incident  perturbations  or 
structured  perturbations)  for  -which  the  guaranteed  bounds  do  not  apply. 

In  the  area  of  fundamental  system  theory  basic  to  control  strategy 
development,  three  general  topics  were  the  subject  of  investigation  thi3  past 
year:  (l  )  the  input-output  behavior  of  3ome  nonlinear  feedback  loops,  (2)  the 

invertibility  of  input-output  systems,  and  (3)  the  existence  and  properties  of 
"strange  attractors"  in  system  theory.  We  give  below  a  brief  account  of  these 
studies,  concentrating  on  the  last  one,  since  it  is  new. 

1  .  We  have  treated  the  "locking"  of  a  class  of  nonlinear  feedback  loops 
to  a  certain  set  of  input  3ignal3;  sufficient  conditions  for  the 
locking  have  been  obtained.  Further,  sufficient  conditions  for  the 
"tracking"  of  the  input  signal  by  the  output  have  been  found.  These 
results  are  being  extended;  a  manuscript  is  in  preparation. 

2.  liven  an  output,  what  wa3  the  input  to  the  system?  It  is  thi3  which 
constitutes  the  question  of  invertibility.  Some  partial  results  have 
been  achieved,  in  the  sense  that  sufficient  conditions  for  the 
invertibility  of  a  class  of  input-output  systems  have  been  found. 
These  conditions  are  fairly  straight-forward  extensions  of  previous 
work.  Currently,  efforts  are  being  made  to  generalize  the  conditions 
30  that  they  are  applicable  to  a  wide  class  of  systems. 

3.  A  variety  of  3y3tems  are  suspected  of  exhibiting  "strange"  behavior. 
3y  this  i3  meant  something  of  the  following  sort  (the  definitions 
varying  according  to  the  investigator):  There  exist  trajectories  of 
the  system  which,  starting  close  to  each  other,  diverge  exponentially 
initially;  these  trajectories  remain  bounded,  however,  for  all  time. 
Finally,  the  trajectories  approach  a  3et  of  points  ''the  "strange 
attractor"1  which  are  not  periodic  orbits  or  a  finite  set  of  singular 
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?o  i.ats . 

The  forced  van  ier  Pol  aquation  with  vary  large  damping  and 
proportionately  large  sinusoidal  forcing  and  similar  such  equations 
are  the  only  differential  equations  for  which  it  has  been  proven  that 
a  strange  attractor  exists.  Numerical  solutions  of  the  Lorenz  system 
and  of  a  version  of  the  forced  Duffing  equation  exhibit  behavior  much 
like  that  of  strange  attractors,  but  proofs  of  this  are  absent. 

It  has  been  suggested  that  some  of  the  properties  of  turbulent 
mixing  and  of  a  number  of  other  phenomena  could  be  explained  if 
strange  attractors  existed  for  the  governing  equations. 

We  have  started  to  study  strange  attractors,  by  repeating  some 
of  the  numerical  work  reported  in  the  literature,  and  by  a  review  of 
the  theory.  We  were  unable  to  duplicate  the  results  on  the  forced 
Duffing- like  equation  for  a  long  time  -until  we  used  an 
Adams-Sashforth  predictor-corrector  integration  routine  for  the 
numerical  solutions  of  the  differential  equation;  now  our  results 
agree  with  those  published.  We  are  in  the  initial  stage  of  what  may 
be  a  long  term  project. 


Well-posedness  of  leader- follower  strategies  obtained  from  singularly 
perturbed  reduced  order  models  has  been  investigated  in  a  manner  similar  to 
our  approach  previously  developed  for  Nash  strategies  for  reduced  order 
models.  We  have  derived  order  reduction  procedures  which  lead  to  well-posed 
formulations  of  the  leader-follower  control  of  singularly  perturbed  dynamic 
systems.  We  have  investigated  the  strong  connection  between  information 
structure  and  well-posedness  of  singularly  perturbed  Nash  games  and  we  have 
identified  why  the  natural  order  reduction  is  ill-posed  as  we  reported 
earlier.  Ither  structures  which  preserved  information  were  demonstrated  to 
lead  to  well-posed  solutions.  Details  are  described  in  the  journal  and 
meeting  papers. 

The  problem  of  controlling  a  linear  stochastic  system  so  a3  to  minimize  a 
quadratic  steady  3tite  cost  functional  with  uncertain  process  and  observation 
noi3e  statistics  was  investigated.  The  problem  was  formulated  as  a  minimax 
optimization.  We  showed  that  the  minimax  solution  is  obtained  by  the  usual 
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linear  quadratic  deterministic  gain  multiplying  the  state  estimate  provided  by 
a  linear  filter.  The  linear  filter  is  a  Kalman-Busy  filter  for  a  particular 
worst  case  pair  of  spectral  density  matrices.  These  results  parallel  those  of 
the  minimax  state  estimation  problem  discussed  in  Area  15. 

17.2  Implementation  Constrained  Decomposition  and  Hierarchical  Control* 

In  large-scale  systems,  control  task3  may  be  decomposed  and  hierarchical 
levels  may  be  imposed  on  the  basis  of  analytically  derived  strategies.  Thi3 
project  is  devoted  to  the  fundamental  problem  of  investigating  an  analytical 
framework  for  incorporating  realistic  constraints  of  a  computer  network. 
Questions  concerning  information  flow,  time  scales,  and  loss  of  feedback  loops 
are  examined  in  coordinating  the  multiple  controllers  in  a  network. 

We  have  investigated  and  developed  several  classes  of  hierarchical 
optimization  algorithms.  Classes  of  algorithms  which  exploit  a  weak 
structural  coupling  between  variables  of  an  optimization  problem  were  analyzed 
using  local  convergence  theory.  Natural  decompositions  were  defined  in  terms 
of  splittings  which  directly  exploit  the  structure  of  the  problem  in  the 
decomposition.  We  demonstrated  that  iterations  defined  by  natural 
decompositions  would  converge  whenever  the  coupling  is  sufficiently  weak.  We 
also  have  developed  a  hierarchical  optimization  algorithm  which  can  be  used  to 
solve  the  linear  quadratic  robust  control  problem.  The  algorithm  exploits  the 
structure  of  the  dual  problem  to  define  a  simple  hierarchical  algorithm  which 
does  not  require  a  stabilizing  initial  point  and  which  has  guaranteed 
convergence  properties. 

A  common  feature  in  large  scale  system  practice  i3  that  different 

decision  makers  assume  different  simplified  low-order  models  of  the  same  large 

scale  system.  This  may  be  due  to  the  necessity  of  easing  the  computational 

burden  in  evaluating  strategies  or  the  lack  of  adequately  modeled  dynamics  of 

3ome  part3  of  the  large  scale  3y3tem.  Based  on  these  simplified  models,  the 

decision  makers  would  then  attempt  to  arrive  at  low-order  decentralize i 

strategies  which  would  be  economically  feasible  to  implement,  at  the  same  time 

preserving  near-optimality.  The  problem  of  designing  control  strategies  for 

decision  makers  under  a  situation  where  they  have  different  models,  different 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program 
('J.3.  Army,  U.3.  Navy,  and  U.3.  Air  Force)  'under  Jontraet  N0001 4--79-T-0424. 
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information  sets,  and  different  objective  functional3,  has  been  addressed 
previously  within  the  framework  of  multiparameter  singular  perturbations.  In 
this  framework,  a  large  scale  system  i3  viewed  as  consisting  of  a  "slow"  core 
coupled  to  a  number  of  "fast”  subsystems.  A  multimodel  situation  results  when 
each  decision  maker  retains  the  dynamics  of  one  fast  subsystem  and  the  3low 
core,  while  neglecting  the  dynamics  of  the  remaining  fast  subsystems.  This 
year  the  general  problem,  wherein  the  fast  subsystems  need  not  be  weakly 
coupled  and  each  fast  subsystem  might  be  affected  by  more  than  one  decision 
maker,  has  been  addressed.  A  procedure  has  been  formulated  to  obtain  near 
optimal  decentralised  Nash  strategies  in  the  two  player  case  under  a 
multimodel  situation  'within  the  framework  of  multiparameter  singular 
perturbations. 

1 7 . 3  Control  Strategies  for  Complex  Systems  for  ’J3e  in  Aerospace  Avionics* 

Whenever  model  -uncertainty  i3  present  or  a  range  of  operating  conditions 
i3  anticipated,  engineering  analysis  and  design  must  deal  with  questions  of 
sensitivity.  Performances  of  manufactured  components  are  necessarily 
specified  with  nonzero  tolerances.  The  parameter  values  which  characterize 
these  components  are  therefore  uncertain.  These  parameter  values  may  change 
during  the  operating  lifetimes  of  the  components  due  to  aging  and  due  to 
changes  in  environmental  conditions,  3uc'n  a3  pressure  and  temperature.  In 
addition,  mathematical  models  used  for  analysis  3nd  design  of  actual  systems 
cannot  possibly  lead  to  predicted  performance  which  exactly  matches  the 
performance  of  the  actual  systems.  It  is  useful  to  regard  some  of  the 
parameters  of  the  models  as  uncertain  in  order  to  make  these  models  more 
realistic.  The  combined  effects  of  parameter  uncertainties  on  overall  system 
behavior  are  of  principal  consideration  in  any  system  design. 

We  investigate  uncertainty  using  three  techniques.  Dne  point  of  view 
regards  the  parameters  a3  unknown  but  deterministic,  and  studies  the 
uncertainty  situation  using  parameter  sensitivity  methods.  New  results  here 
include  establisment  of  some  relationships  between  singular  value  robustness 
approaches  and  the  now-standard  comparison  sensitivity  methods.  Design- 
oriented  techniques  are  being  developed  for  large  parameter  variations  in 
linear  and  nonlinear  systems. 

*  Thi3  work  was  supported  by  the  'J.3.  Air  Force  under  leant  AFD3R-73-d673. 
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A  second  approach  models  the  parameters  as  random  variables.  The  tools 
for  this  approach  are  stochastic  control  and  dual  control  theory,  including 
stochastic  multi-controller  situations.  The  interplay  between  information 
structure  and  control  has  been  emphasized  in  our  recent  work. 

Ve  have  expanded  our  basic  research  to  encompass  an  area  of  increasing 
importance  in  control  research  and  aerospace  applications  -  adaptive  control. 
Our  approach  focuses  on  the  role  of  singular  perturbation  theory  in  the 
development  of  model  reference  adaptive  systems  when  the  plant  and  reference 
model  have  different  orders.  An  important  possible  interaction  has  been 
discovered  between  frequencies  present  in  a  "suff iciently  rich"  adaptive  input 
signal  and  neglected  high  frequency  terms  in  a  reduced -order  model  in  a  model 
reference  adaptive  system.  Details  of  these  results  may  be  found  in  the 
journal  articles  and  conference  papers  listed  in  the  publications  section. 

17-4  Large  Beale  Bystems* 

Several  new  results  have  been  obtained  in  multimodeling,  that  is,  in 
modeling  and  control  situations  in  which  different  decision  makers  employ 
different  models  of  the  same  system.  Multiparameter  singular  perturbation 
techniques  are  employed  to  capture  the  multimodel  nature  of  fast  dynamic 
subsystems  interconnected  through  slow  variables.  In  one  problem  studied,  a 
decentralized  filtering  and  control  scheme  was  developed  for  obtaining  low 
order  Nash  equilibrium  strategies  in  a  multimodel  setting.  The  decision 
makers  have  decentralized  information  structures  and  are  constrained  to  use 
finite  dimensional  compensators.  In  another  study,  physical  examples  (s.g., 
networks')  are  used  to  reexamine  the  role  of  modeling  in  large  scale  system 
analysis  and  design.  A  general  property  of  the  systems  considered  is  that 
they  are  strongly  coupled  (physically)  in  the  slow  time  scale,  and  weakly 
coupled  in  the  fast  time  scale.  This  leads  to  a  multimodel  situation  in  which 
every  subsystem  controller  can  neglect  all  other  fast  subsystems  except  his 
own. 

In  the  area  of  stochastic  multiplayer  strategies,  an  indirect  method  to 

obtain  the  Btackelberg  solution  of  two-person  nonzero-sum  differential  games 

with  general  dynamics  and  cost  functionals  has  been  developed  with  the  leader 

*  This  work  wa3  supported  by  the  National  Science  Foundation  under  Brant  ECS 
79-1  9-596- 
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having  access  to  sampled  values  of  the  state.  This  new  method  involves  a 
three-step  optimization,  and  converts  the  original  dynamic  problem  into  two 
open-loop  optimal  control  problems  and  one  static  finite-dimensional 
optimization  problem.  Solutions  of  these  optimization  problems  are  obtained 
as  solutions  of  a  set  of  differential  equations  with  mixed  boundary 
conditions,  which  in  turn  leads  to  the  Stackelberg  cost  value  of  the  leader, 
as  wall  as  to  a  class  of  strategies  that  yield  that  cost  level. 

Chained  aggregation,  an  iterative  approach  to  large  3cale  system  model 
reduction,  has  been  developed  further.  The  resulting  system  representation, 
the  Generalized  Hessenberg  Representation  (GHR),  has  been  shown  to  exhibit 
explicitly  the  system  observability  structure.  Geometric  (in  the  sense  of 
Vonham)  properties  of  the  GHR  representation  have  been  exposed  and  exploited 
in  the  development  of  a  three-component  decentralized  hierarchical  feedback 
design  technique. 

Finally,  coherency,  a  concept  widely  used  in  power  system  studies,  has 
been  examined  from  an  aggregation  point  of  view.  Specifically,  the 
relationship  among  coherency,  aggregation,  and  observ ability  has  been 
clarified,  and  a  fundamental  relationship  between  "slow  coherency"  and  weak 
coupling  has  been  established.  Details  of  these  and  other  results  are  given 
in  the  list  of  .journal  articles. 

17.5  Power  System  Normal  and  Security  State  Assessment* 

Normal  and  security  state  assessment  of  large  scale  power  systems 
requires  the  analysis  of  power  system  performance  under  all  levels  of  likely 
load  and  contingencies.  Performance  is  measured  in  terms  of  bu3  voltage 
levels,  line  power  flows,  stability  margins,  and  their  sensitivities  to  large 
changes  in  network  conditions.  Explicit  load  flow  solutions  have  been 
ieveloped  to  provide  near  closed  form  analysis  of  network  response  to  total 
load.  These  explicit  forms  may  be  used  to  produce  time  varying  solution 
trajectories  and  exact  stochastic  load  flow  solutions. 


*  This  work  was  supported  by  the  National  Ocience  Foundation  under  Grant  ENG 
75-05594,  and  the  University  of  Illinois  Power  Affiliates  Program. 
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Many  optimal  load  flow  algorithms  require  stability  margin  constraints. 
These  constraints  have  typically  hinged  largely  on  bus  voltage  angle  limits. 
When  finite  VAR  control  is  included,  these  angle  limits  have  been  shown  to  be 
erroneous  and  in  some  oases  unfeasible.  Extension  of  these  realistic 

stability  margins  to  an  arbitrary  network  i3  currently  being  considered. 

17.6  Control  of  Stochastic  Systems  Containing  Parameter  Uncertainty* 

The  main  result  of  this  research  is  the  determination  of  algorithms  for 
Sensitivity  Adaptive  Feedback  with  Estimation  Redistribution  (SAFER)  control. 
This  type  of  control,  which  we  proposed  and  developed,  has  dual  properties  in 
the  3ense  that  the  control  takes  into  account  future  parameter  estimation 

accuracy.  The  parameter  sensitivities  are  altered  by  changing  the  sensitivity 
weighting  coefficients  in  a  composite  performance  index.  Our  early  work  on 
thi3  problem  dealt  with  a  general  multivariable  state  space  formulation  which 
led  to  quite  complicated  algorithms.  We  continually  simplified  the  control 
algorithms  by  first  considering  input-output  representations,  and  finally 
single- input  single-output  representations.  Much  of  the  application  of 
adaptive  control  is  for  this  latter  class  of  problems.  Very  substantial 
reductions  in  computation  are  achieved  not  only  because  of  the  minimal  number 
of  ’unknown  parameters  in  the  single- input  single-output  case,  but  also  because 
of  our  use  of  only  two  sensitivity  functions  to  capture  the  impact  of  all  the 

parameters  in  a  SAFER  control  setting.  Our  simpler  algorithm  for  SAFER 

control  has  been  demonstrated  to  yield  performances  which  are  superior  to  the 
certainty  equivalent  control. 

This  project  is  a  cooperative  effort  between  the  Coordinated  Science 
Laboratory  of  the  University  of  Illinois  and  the  Systems  Engineering 
Laboratory  of  the  Venezuelan  Institute  for  Scientific  Research. 


*  This  work  was  supported  by  the  National  Science  Foundation  under  Orant  INT 
77-20969. 
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17.7  Aggregated  Hierarchical  Control* 

The  objective  of  this  research  is  to  investigate  hierarchical 
decomposition  algorithms  for  which  the  problem  solved  at  each  level  of  the 
hierarchy  is  approximated.  Standard  hierarchical  decomposition  theory 
requires  that  an  optimization  problem  be  decomposed  exactly  and  that  the 
iteration  defined  by  the  decomposition  converge  to  the  exact  solution.  Ve 
have  relaxed  the  exactness  requirement  by  allowing  the  exact  problems  at  each 
hierarchical  level  to  be  replaced  by  an  operator  theoretic  aggregate.  The 
hierarchical  structure  is  determined  by  the  composition  of  an  exact 
decomposition  and  an  aggregation.  Conditions  'under  which  the  approximate 
decomposition  iterates  converge  to  an  approximation  of  the  exact  solution  have 
been  determined  for  linear  problems.  The  extension  of  these  results  to 
nonlinear  problems  through  the  use  of  nonlinear  aggregation  and  nonlinear 
splitting  functions  is  under  current  investigation. 


*  This  work  was  supported  by  the  National  Science 
ENG-79-03773. 
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15.1  Ari  Investigation  of  Residue  Number  Architectures  for  Digital  Filters 

with  Fault  Tolerance* 

This  project  involves  a  study  of  how  the  theory  of  redundant  residue 
number  systems  (RNS)  can  be  developed  and  applied  to  the  isolation  and 

detection  of  hardware  failures  in  digital  processors,  -with  subsequent  "30ft 
failure”,  that  allows  a  faulty  processor  to  continue  operating  with  reduced 
capabilities.  The  class  of  systematic  redundant  RNS  codes  was  analyzed  and 
simulated,  thereby  establishing  preferred  design  procedures  and  demonstrating 
failure  resistant  behavior  in  recursive  digital  filters  [l3. 121.  At  this 

time,  the  theory  of  failure  resistant  design  for  the  systematic  RNS  codes  i3 
relatively  complete,  including  bounds  on  the  necessary  redundancy  and 
procedures  for  identifying  faulty  modules  and  for  reallocating  resources 
following  soft  failure.  The  major  issues  now  concern  reliability  of  the  error 
checking  hardware  itself,  as  well  as  new  techniques  for  minimizing  the 

*  This  work  was  supported  by  the  National  Science  Foundation  under  Grant 

ENG-79-01 086. 
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complexity  and  cost  of  the  overhead  hardware. 

Progress  has  been  reported  in  the  design  of  complex  residue  arithmetic 
for  use  in  spectral  analysis,  processing  of  bandpass  signals  represented  in 
analytic  form,  and  in  complex  cepstral  analysis  for  speech 
[13.13,13.30,13.32].  Also  a  design  has  been  presented  which  permits  the  block 
elements  of  a  general  mixed  radix  converter  to  be  controlled  in  such  a  way 
that  the  general  converter  can  be  electronically  adapted  to  sequentially 
generate  all  projections  that  are  required  for  error  location  [ 13.26].  Since 
this  is  a  significant  step  toward  being  able  to  implement  an  entire  error 
checker  on  a  single  VLSI  integrated  circuit,  a  great  deal  of  effort  is  being 
expended  to  extend  this  concept  in  order  to  bring  failure  resistant  concepts 
into  practical  usage. 

A  new  effort  has  been  started  to  evaluate  the  advantages  and  to  extend 
the  theory  of  redundant  RNS  product  codes  for  digital  signal  processing. 
Rather  than  simply  adding  redundant  residue  digits,  data  is  encoded  in  a 
product  code  by  multiplication  with  a  code  generator.  .An  RNS  product  code  is 
an  arithmetic  code  which  is  preserved  'under  error  free  arithmetic  operations, 
’.fnile  not  too  much  is  known  about  the  theory  of  RNS  product  codes  in  signal 

processing,  initial  indications  are  that  they  have  useful  properties  worthy  of 

further  study. 

13.2  Roundoff  Noise  and  Limit  Cycles  in  Digital  Filtering* 

It  is  well  known  that  canonical  implementations  of  narrowband  recursive 
digital  filters  can  exhibit  large  error  due  to  arithmetic  roundoff.  However, 
this  error  can  be  reduced  at  the  expense  of  additional  hardware  complexity. 
For  example,  Mullis  and  Roberts  have  developed  a  procedure  for  finding  the 
minimum  roundoff  noi3e  linear  3tate-space  (LSS)  filter  structure  _R1  1. 
Another  realization  for  reducing  errors  in  digital  filters  i3  the  error 

spectrum  shaping  (ESS)  structure  [R2,R3].  The  ESS  structure  utilizes 
quantizer  feedback  and  generally  has  a  simpler  Implementation  than  the  L33 
structure. 

*  Thi3  work  was  supported  by  the  Joint  Services  Electronics  Program 

C * J . 3 .  Army,  H.S.  Navy,  j.S.  Air  Force)  under  Contract  N0001 4-79-C-04-24-. 
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We  have  compared  the  noise  reduction  performance  of  the  ESS  structure 
with  that  of  the  optimal  LSS  structure  for  second  order  digital  filter 
sections.  It  was  found  that  optimal  direct  form  1  and  direct  form  2  ESS 
realizations  always  have  a  higher  signal- to-noise  ratio  than  the  optimal  LSS 
structure.  Several  suboptimal  ESS  structures  with  simple  hardware 
implementations  were  also  considered  and  optimal  values  for  the  ESS 
coefficients  in  these  structures  were  derived.  For  filters  with  zeroes  at  z  * 
-1  ,  it  was  found  that  several  of  the  simple  ESS  structures  can  outperform  the 
optimal  LSS  structure.  For  elliptic  filters  it  was  found  that  the  suboptimal 
direct  form  2  ESS  structures  perform  poorly,  but  that  the  direct  form  1  ESS 
structures  perform  well.  We  are  currently  investigating  high-order  ESS 
structures  and  developing  ESS  implementations  with  low  coefficient 
sensitivity. 

Our  wori:  concerning  limit  cycles  covers  two  topics.  First,  we  have 
studied  the  concept  of  accessibility,  introduced  by  Claasen  _et  al.  ^34-1.  A 
zero- input  limit  cycle  is  accessible  if  it  can  be  reached  from  a  state  which 
is  not  a  state  of  the  limit  cycle.  It  is  important  to  study  accessible  limit 
cycles,  since  they  are  most  frequently  observed.  We  have  shown  that 
accessibility  places  restrictions  on  the  error  sequences  introduced  at  the 
filter  quantizers.  In  the  case  of  two  rounding  quantizers,  these  restrictions 
can  be  used  to  derive  a  lower  bound  on  the  minimum  amplitude  accessible  limit 
cycle.  A  related,  but  different,  result  holds  for  two  truncating  quantizers. 
We  have  also  shown  that  only  certain  states  can  directly  access  a  limit  cycle 
and  that  large  limit  cycles  tend  to  be  accessible. 

The  second  limit  cycle  topic  is  that  of  determining  exact  maximum 
amplitude  limit  cycles  in  digital  filters.  We  have  recently  devised  an 
algorithm  for  finding  the  maximum  amplitude  limit  cycle  of  a  particular 
period.  This  is  accomplished  by  considering  all  possible  initial  states, 
which  can  be  maximum  amplitude  states  of  a  limit  cycle.  Results  have  been 
obtained  which  significantly  reduce  the  number  of  states  requiring 
consideration,  thereby  resulting  in  a  computationally  feasible  algorithm. 
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1 3. 3  -lew  Techniques  in  Digits!  Signs!  Processing  for  Synthetic  Aperture  Radar* 

Research  in  synthetic  sperture  radar  (SAR'l  has  been  conducted  in  four 
distinct,  but  related  projects:  (l )  analysis  of  inherent  phase  distortion  in 
rectangular  format  FFT  algorithms,  (2)  a  SAR  computer  simulation  with  polar 
format  FFT  processing,  (3)  a  study  of  2D  interpolators  for  polar-to- 
rectangular  transformation,  and  (l)  an  investigation  of  number  theoretic 
concepts  for  high  speed  failure  resistant  digital  processors  required  in 
real-time  3AR.  Recently,  special  emphasis  was  placed  on  analysis  of  spotlight 
mode  operation,  including  requirements  for  2D  spatially-varying  interpolators 
for  polar- to-rectangular  coordinate  transformation  required  in  FFT  processing. 
Five  2D  interpolators  have  been  analyzed  and  simulated:  ( 1 )  1st  order  inverse 
distance  ( nonseparable)  ,  (21  1st  order  separable  bilinear,  ' 3 'I  3rd  order 
separable  Lagrange,  ( A )  separable  optimal  FIR  filtering,  and  (o'*  zero  order 
interpolation  with  over-sampling.  Attention  was  concentrated  on  low  order 
interpolators  because  they  are  simple  enough  to  implement  in  practice. 
Although  many  questions  remain,  it  appears  that  the  use  of  optimal 
space- invariant  FIR  interpolation  coupled  'with  nearest-neighbor  selection 
(zero  order  interpolation)  may  be  the  most  efficient  computational  means  to 
achieve  high  resolution  performance. 

It  was  discovered  that  in  polar  format  FFT  3AR  processing,  each  line  of 
recorded  data  is  the  Fourier  transform  of  a  density  integral  of  the  target 
area,  collected  at  a  specified  look  angle.  By  interpreting  the  data  recording 
in  thi3  'way,  it  was  possible  to  establish  that  the  SAR  principle  does  not 
depend  on  instantaneous  motion  between  the  target  and  the  radar,  but  rather  on 
the  total  view  angle  through  which  projections  are  taken.  This  principle 
seems  to  be  obscured  in  the  present  literature,  where  SAR  is  discussed  in 
terms  of  Doppler  frequency  analysis,  thereby  implying  a  dependence  on 
instantaneous  Doppler  frequency.  Recent  work  has  concentrated  on  learning  how 
reconstruction  can  be3t  be  accomplished  from  incomplete  frequency  domain  data, 
which  results  in  SAR  due  to  limited  view  angles. 


*  This  work  was  suDoorted  by  the  Air  Force  Office  of  Scientific  Research  under 

Drant  A?03R-79-D02§A. 
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1 3.4  New  Directions  in  Multidimensional  Signal  Processing* 

We  have  recently  begun  a  program  of  basic  research  in  multidimensional 
signal  processing  spanning  the  areas  of  digital  array  processing,  linear  space 
varying  processing,  and  nonlinear  processing  using  order  statistics.  In  array 
processing,  efforts  are  being  concentrated  on  the  development  of  a  common 
mathematical  framework  for  computer  aided  tomography  (CAT)  and  spotlight-mode 
synthetic  aperture  radar  (SAR).  The  underlying  principle  in  CAT  is  the 
projection-slice  theorem.  Ve  have  shown  that  this  same  theorem  also  forms  the 
basis  for  spotlight-mode  SAR;  hence  similar  processing  is  required  in  both 
systems.  Future  work  will  explore  the  possibility  of  applying  efficient  CAT 
algoritms  to  SAR  and  vice  versa.  In  the  area  of  linear  space  varying 
processing  we  have  extended  recent  results  concerning  required  sampling  rates 
for  time  varying  digital  systems.  Unlike  the  time  invariant  case,  the 
required  sampling  rate  depends  on  both  the  input  bandwidth  and  the  desired 
analog  impulse  response.  We  are  presently  investigating  the  use  of  simple 
time  varying  digital  filters  to  perform  time  invariant  filtering.  In  the  area 
of  nonlinear  filtering  we  have  studied  a  generalized  median  filter  whose 
output  is  a  weighted  sum  of  order  statistics.  Assuming  a  constant  signal  in 
white  background  noise,  a  procedure  has  been  developed  for  finding  optimal 
weight  values.  Deterministic  properties  of  the  generalized  median  filter  and 
various  possible  recursive  order  statistical  filters  are  also  being  studied. 

1 3. 5  Image  Sequence  Processing  and  Dynamic  Scene  Analysis** 

The  processing  of  images  involving  motion  has  become  increasingly 
important.  Applications  include  target  tracking,  dynamic  robot  vision,  image 
bandwidth  compression,  dynamic  medical  imaging,  and  highway  traffic 
monitoring.  Image  sequence  processing  involves  a  large  amount  of  data. 
However,  because  of  the  rapid  progress  in  computer,  LSI,  and  7  LSI 
technologies,  many  useful  processing  task3  for  image  sequences  can  now  be 
performed  in  a  reasonable  amount  of  time. 

*  This  -work  was  supported  by  the  Joint  Services  Electronics  Program 
(7.3.  Army,  U.S.  Navy,  IT.3.  Air  Force)  under  Gontract  N0001 4-79-C-042A. 

**  This  work  was  supported  by  the  National  Science  Foundation  'under  Grant 
MCS-73- 13906. 
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A  central  issue  in  image  sequence  processing  i3  motion  estimation. 
During  the  past  year  we  have  developed  a  number  of  efficient  algorithms  for 
estimating  the  three-dimensional  motion  parameters  of  rigid  bodies  (both 
translation  and  rotation)  from  two  time-sequential  image  frames.  The 
following  uniqueness  theorems  have  also  been  proven:  (i)  ?or  a  rigid  planar 
patch,  two  solutions  for  the  motion  parameters  will  exist  in  general, 
(ii)  Por  a  rigid  curved  surface,  the  solution  for  the  motion  parameters  i3 
unique . 


1 3. 5  Underwater  Image  Processing* 

We  have  been  collaborating  'With  the  Naval  Research  Laboratory  on  a  basic 
research  project  in  underwater  image  processing  using  both  optical  and  3onar 
sensors.  Applications  include  searching  for  man-made  objects  on  the  ocean 
floor  and  ocean  mining.  During  the  past  year,  our  research  has  concentrated 
on  studying  the  effects  of  ocean  water  on  optical  and  sonar  images  and  on 
finding  3hape  and  texture  features  of  images  of  the  ocean  floor  that  can 
distinguish  between  man-made  and  natural  objects.  Successful  algorithms  for 
extracting  3hape  features  nave  been  developed.  We  have  also  studied  the 
characteristics  of  various  acoustical  imaging  methods  and  concluded  that 
synthetic  aperture  sonar  and  holography  are  the  two  most  promising  approaches 
for  'underwater  imaging. 
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19-1  Summary  of  Progress 

This  group  has  two  primary  areas  of  interest.  The  first  is 
computer-assisted  education  in  population  and  food  and  their  impact  on  other 
social  and  economic  areas  of  society.  The  second  is  concerned  with  climate 
and  food  and  the  impact  of  large  interannual  variations  of  climate  on  the 
availability  of  food  for  the  world's  growing  population. 

In  the  first  area,  population  and  food,  we  have  designed  a 
computer-assisted  education  program  concerned  with  world  population  as  well  as 
specific  countries.  The  program  is  available  in  BASIC  and  has  been  adapted 
for  use  on  the  Apple  computer.  The  program  is  used  in  over  200  educational 
institutions  each  year.  Many  of  the  programs  are  used  3t  the  high  school 
level  where  distribution  of  the  programs  has  just  begun  during  the  past  year. 

In  the  second  area,  food  and  climate,  a  retrospective  study  has  been  male 
of  the  relationship  of  crop  production  to  climate  anomalies.  Various 
interesting  climate  sequences  have  been  observed  where  a  certain  type  of 
climate  anomaly  in  one  region  of  the  world  is  very  frequently  followed  by 
another  anomaly  in  another  part  of  the  world. 

*  Thi3  work  was  supoorted  by  the  National  Science  Foundation  under  Crant  SHD 
vq_i 3501  and  United  States  Department  of  Agriculture  Contract  SEA 
53-5193-0-334. 
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These  sequences  of  climate  anomalies  can  be  usei  as  a  means  of  crop 
forecasting  in  the  region  of  the  later  anomaly.  Various  techniques  are  being 
ievelopei  to  attempt  to  recognize  these  sequences  in  the  world  climate  pattern 
of  the  past  twenty  years.  When  special  patterns  are  observe!  they  are  tested 
against  longer  data  bases.  Some  of  the  iata  are  contained  in  time  series 
which  go  back  to  1  365.  A  paper  iescribing  the  relationship  of  crops  in 
Australia  and  Indonesia  to  the  climate  anomalies  of  the  tropical  Pacific  Ocean 
has  been  accepted  for  publication  by  the  journal  Olimate  Ohange  ^ l  9*  d 1  • 
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